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CHAPTER  1 


INTRODUCTION 


1.1  Background 

The  US  Air  Force  has  invested  heavily  in  the  development  of  advanced  surveillance  systems  and  technologies. 
Of  increasing  concern  is  the  threat  that  potential  adversaries  may  use  bistatic  technologies  to  take  advantage  of  our 
significant  investment  in  advanced  sensors.  With  relatively  inexpensive  receiver  systems,  an  adversary  could  host 
off  of  our  illuminations  to  operate  bistatically,  A  central  requirement  for  non-cooperative  bistatic  operation  is  the 
estimation  of  a  coherent  reference  signal.  This  estimate  is  used  to  correlate  against  received  waveforms  to  extract 
the  desired  signal.  As  illustrated  in  Figure  1.1,  a  coherent  reference  is  typically  obtained  by  measuring  a  direct  path 
signal  through  the  sidelobes  of  the  illuminator.  Conventional  methods  to  prevent  the  interception  of  the  direct  path 
signal  include  low  sidelobe  antennas,  physical  isolation,  and  the  use  of  spread  spectrum  waveforms.  These  methods 
will  become  inadequate  as  surveillance  sensors  migrate  to  space. 


Figure  1.1:  Non-cooperative  bistatic  receivers  require  a  coherent  reference  from  the  host  illuminator 


In  addition,  to  enhance  precision  engagement,  it  is  desirable  that  covert  data/voice  communication  be 
incorporated  simultaneously  with  the  surveillance  function.  This  embedded  communications  capability  should 
allow  for  the  message  signal  to  vary  with  angle  such  that  separate  information  can  be  directed  to  different  receivers 
located  at  various  angles.  Directionally  dependent  waveforms  would  be  ideal  for  space  based  sensing  of  RF  tags  on 
aircraft  (IFF),  ground  equipment  and  troops  and  would  be  especially  valuable  for  search  and  rescue  missions.  Of 
course,  the  radar  and  communications  waveforms  should  be  designed  so  that  they  can  be  radiated  at  the  same  time 
without  self-jamming. 

1.2  Proposed  Methodology 

As  shown  in  Figure  1 .2,  it  is  proposed  to  employ  an  (N+2)  element  antenna  array  where  N  of  the  elements  are 
used  by  the  host  radar  for  surveillance  purposes  while  two  of  the  elements  are  used  by  the  interferometer  for 
coherent  reference  denial  and  communications  purposes.  The  two  interferometric  elements  are  driven  separately 
with  an  independent  waveform  generation,  timing  and  control  circuit. 

Ideally,  the  interferometer  antenna  pattern  will  overlay  the  sidelobes  of  the  host  radar  main  antenna  pattern  with 
minimal  overlay  of  the  radar  main  beam.  The  interferometer  may  or  may  not  be  on  all  of  the  time.  However,  as  a 
minimum,  the  interferometer  will  be  on  while  the  host  radar  is  in  the  transmit  mode  with  the  objective  that  the 
interferometer  signal  will  mask  that  portion  of  the  host  radar  signal  emitted  through  the  radar  sidelobes.  In  this  way, 
a  coherent  reference  signal  is  denied  to  a  non-cooperative  bistatic  receiver. 
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Figure  1.2:  A  pair  of  interferometric  elements  driven  by  independent  waveform  generation,  timing  and  control  provides  a  means 

to  modulate  the  sidelobes  of  the  host  radar  main  array. 


To  increase  the  effectiveness  of  this  masking,  it  is  also  proposed  to  modulate  the  interferometer  antenna  pattern 
from  pulse  to  pulse  such  that  the  pattern  is  rotated  on  each  pulse.  This  is  illustrated  in  Figure  1 .3.  Due  to  the  nature 
of  the  interferometer  modulation,  the  proposed  method  of  denying  a  coherent  reference  is  more  effective  than 
masking  the  sidelobes  with  a  noise  signal. 

Also,  by  carefully  selecting  the  interferometer  modulation,  a  desired  message  signal  can  be  embedded  in  the 
sidelobes,  thereby  allowing  data/voice  communications  operations  simultaneously  with  the  surveillance  function. 
Using  a  multi-channel  phased  array  antenna  and  orthogonal  spatial-temporal  waveforms,  the  message  signal  can 
vary  with  angle.  In  this  manner,  separate  communications  signals  can  be  directed  to  different  receivers  located  at 
various  angles. 

13  Research  Objectives 

The  objectives  of  this  research  effort  are  summarized  below: 

1)  Develop  joint  spatial-temporal  antenna-based  signal  processing  techniques  and  waveforms  that  can  be  applied 
simultaneously  to  radar  and  communications  systems  without  self  jamming, 

2)  Design  multi-dimensional  waveforms  for  the  purposes  of  denying  a  coherent  reference  to  a  non-cooperative 
bistatic  radar  located  at  an  arbitrary  angle  while  simultaneously  permitting  radar  surveillance  and  multiple 
embedded  communications, 

3)  Design  diverse  waveforms  which  will  accommodate  multi-mission  operations  such  as  ground  and  airborne 
moving  target  indication,  tracking,  automatic  target  recognition,  and  foliage  and  ground  penetration, 

4)  Investigate  theoretical  measures  of  performance  for  sensor  and  communication  systems  incorporating  these  new 
and  novel  waveforms. 


Figure  13:  Pulse-to-pulse  phase  modulation  of  the  sidelobes  denies  a  coherent  reference  to  non-cooperative  receivers.  The 

sidelobe  modulation,  also,  permits  embedded  communications. 
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1.4  Report  Organization 

Bistatic  radar  fundamentals  are  reviewed  in  Chapter  2.  Various  issues  associated  with  the  successful  operation 
of  a  non-cooperative  bistatic  radar  are  discussed.  Use  of  the  interferometer  signal  to  mask  the  direct  signal  of  the 
host  radar  is  then  considered.  Performance  degradation  of  the  non-cooperative  bistatic  radar  is  shown  to  depend 
upon  the  delay-Doppler  ambiguity  functions  of  the  host  radar  and  interferometer  waveforms  and  their  cross 
ambiguity  function. 

Antenna  patterns  of  the  interferometer  and  host  radar  are  developed  in  Chapter  3.  It  is  noted  that  it  is  possible 
to  achieve  spatial  orthogonality  between  the  host  radar  and  interferometer  by  placing  a  broad  null  of  the 
interferometer  pattern  in  the  direction  of  the  main  beam  of  the  host  radar.  However,  it  is  shown  that  frequency 
hopping  is  required  to  steer  the  interferometer  pattern  while  maintaining  the  broad  null.  Various  waveforms  are 
considered  for  the  signals  driving  the  interferometer.  Conditions  required  by  the  wave  equation  are  investigated  for 
the  quasi-stationary  behavior  of  an  antenna  when  modulated  signals  are  radiated.  Response  of  the  interferometer 
element  pair  when  excited  by  two  different  modulated  signals  is  derived.  This  is  specialized  to  the  situation  where 
both  elements  are  excited  by  the  same  linear  frequency  modulated  pulse.  In  addition,  the  case  is  considered  where 
one  of  the  interferometer  signals  is  an  up  chirp  while  the  other  is  a  down  chirp. 

A  frequency-coded  waveform,  employing  Costas  sequences,  is  proposed  for  the  host  radar  signal  in  Chapter  4. 
The  concept  of  the  Costas  array  is  introduced  and  the  delay-Doppler  ambiguity  function  of  the  Costas  signal  is 
derived.  The  ambiguity  function  is  shown  to  have  a  thumbtack-like  shape  (i.e.,  a  single  narrow  central  peak 
surrounded  by  a  low-level  pedestal). 

Orthogonal  frequency  division  multiplexing  (OFDM)  is  chosen  as  the  modulation  scheme  for  the 
interferometer  in  Chapter  5.  After  a  tutorial  discussion  of  OFDM,  the  delay-Doppler  ambiguity  function  of  the 
OFDM  signal  is  derived. 

Performance  of  the  non-cooperative  bistatic  radar  is  considered  in  Chapter  6  assuming  that  the  host  radar  and 
interferometer  employ  Costas  and  OFDM  signals,  respectively.  This  requires  derivation  of  the  delay-Doppler  cross 
ambiguity  function  between  the  Costas  and  OFDM  signals.  Expressions  are  obtained  for  the  false  alarm  and 
detection  probabilities  of  the  non-cooperative  bistatic  radar. 

Finally,  a  summary  of  results  and  suggestions  for  future  work  are  provided  in  Chapter  7. 
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CHAPTER  2 


NON-COOPERATIVE  BISTATIC  RADAR  CONSIDERATIONS 


2.1  Bistatic  Radar  Fundamentals  [1] 

A  bistatic  radar  is  one  in  which  the  transmitter  and  receiver  are  physically  separated  and,  therefore,  have 
separate  antennas,  as  shown  in  Figure  2.1.  The  north-referenced  coordinate  system,  which  is  illustrated  in  Figure 


Figure  2.1 :  A  bistatic  receiver  employs  separate  antennas  for  transmitting  and  receiving. 


2.2,  is  frequently  used  to  describe  the  bistatic  radar  geometry.  The  bistatic  baseline  is  denoted  by  L  and  designates 
the  distance  separating  the  transmitter  and  receiver.  RT  and  Rr  denote  the  distance  from  the  transmitter  to  the  target 
and  the  distance  from  the  target  to  the  receiver,  respectively.  The  bistatic  angle,  (3  ,  is  measured  at  the  target  and  is 
the  angle  between  the  lines  whose  lengths  are  denoted  by  RT  and  Rr. 


Target 


Figure  2.2:  Bistatic  radar  geometry  is  described  using  the  north-referenced  coordinate  system. 
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The  direct-path  delay  between  the  transmitter  and  receiver  is 


tdp=± 

c 


(2.1) 


where  c  is  the  speed  of  light.  The  total  delay  related  to  the  detection  of  a  point  target  is 

-TP  _  Rt+Rr 


T  '  =  (2.2) 

The  bistatic  delay  is  defined  to  be  the  difference  between  the  total  delay  and  the  direct-path  delay  and  is  given  by 

rp  np  Rt  d"  Rr  —  L 

Tb  =  T  -  T  =  — - ^ - .  (2.3) 


Thus,  the  bistatic  delay  is  a  measure  of  the  bistatic  range  which  is  defined  to  be 


BR  —  R  +  R„  —  L. 

T  K 


(2.4) 


The  locus  of  constant  bistatic  range  is  an  ellipsoid  whose  foci  are  located  at  the  transmitter  and  receiver.  The  major 
semi-axis  length  is  given  by 


a 


L  +  BR  RT  +  Rr 


2  2 

while  the  minor  semi-axis  lengths  are  equal  and  given  by 

b  =  yla2-L2  /4 . 

For  the  special  case  of  a  monostatic  radar, 

L  =  tdp  =0 

and  the  ellipsoid  becomes  a  sphere  of  radius  Rr  =  RT.  The  monostatic  delay  related  to  a  point  target  is 


(2.5) 


(2.6) 

(2.7) 


Z  -  2Rt  -  2R* 
^  m  c  c 


(2.8) 


The  Doppler  frequency  observed  in  the  received  signal  is  proportional  to  the  time  derivative  of  the  bistatic 
range.  Let  the  target  velocity  vector  be  denoted  by  V  .  Assuming  a  stationary  transmitter  and  receiver,  (dL/dt)  =  0. 
The  time  derivative  of  the  bistatic  range  is  then  given  by 


d{BR) 


dt 


=  —2 1  €>  |  cos  8  cos(fi !  2) 


(2.9) 


where  8  is  the  angle  between  the  velocity  vector  and  the  bistatic  angle  bisector,  as  illustrated  in  Figure  2.3.  Hence, 
the  bistatic  Doppler  frequency  is 
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Figure  2.3:  The  bistatic  Doppler  frequency  depends  upon  the  angle  between  the  velocity  vector  and  the  bistatic  angle  bisector. 


r  1  d(BR)  2  v  ~  ~  im 

fB  =  -- - V- -  =  -J— [COs8cos(/3/2)  (2.10) 

at  Ao 

where  A.  0  is  the  wavelength  of  the  propagating  wave.  Observe  that  |  V  |  COS  8  is  the  projection  of  the  velocity 
vector  on  the  bistatic  angle  bisector.  It  represents  the  velocity  component  perpendicular  to  the  constant  range 
ellipsoid.  Along  the  bistatic  baseline,  /?  =  180°  and  the  target  Doppler  frequency  is  always  zero  independent  of  the 
velocity  vector  orientation. 

For  a  monostatic  radar,  (5  =  0°.  The  Doppler  frequency  becomes 

fm=lMcosS  (2.11) 

Ao 


where  8  is  now  the  angle  between  the  velocity  vector  and  the  radial  line  between  the  target  and  the  radar.  The 
monostatic  Doppler  frequency  geometry  is  shown  in  Figure  2.4. 

The  bistatic  triangle  of  Figure  2.2  is  readily  solved  to  locate  the  target  position  relative  to  the  transmitter 
and  receiver.  Three  pieces  of  information  are  needed  to  completely  and  unambiguously  specify  a  triangle,  one  of 

which  must  be  the  length  of  a  side.  Given  L,  then  any  two  of  BR,  0T  ,  and  0R  can  be  used  to  solve  the  triangle  to 
locate  the  target.  For  example,  using  BR  and  0R  , 


Figure  2.4:  The  monostatic  Doppler  frequency  depends  upon  the  angle  between  the  velocity  vector  and  the  radial  line  from  the 

target  to  the  radar. 


(BR)2  +  2L(BR) 

~  2[BR  +  L(\  +  sin  dR  )] 
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In  a  similar  manner,  using  BR  and  dT  , 


In  addition,  from  Figure  2.2, 


Using  the  law  of  sines,  it  follows  that 


(BR)2  +  2  L(BR) 

T  2[BR  +  L(\  -  sin0r)] 

p  =  eT-eR. 


sin  P  _  cos  6r  _  cos  0T 
L  ~  Rt  ~  Rr 


(2.13) 


(2.14) 


(2.15) 


As  a  result,  knowledge  of  any  three  of  the  four  parameters  0T  ,  0R  ,  L,  and  BR  is  sufficient  to  determine  the  fourth. 
Also,  it  can  be  shown  using  the  law  of  sines  that 


and 


(BR  +  L)  cos0r 
cos  9r  +  cos  9t 


_  (BR  +  L)cos9r 

Jvj’  “  , 

COS^J  +  COS0J 

Finally,  the  distance  along  the  baseline  can  be  expressed  as 

L  _  (BR  -F  X)sin(^r  -  0R) 
COS0R  +cos0T 


(2.16) 


(2.17) 


(2.18) 


The  discussion  in  this  section  assumed  that  the  bistatic  radar  transmitter  and  receiver  were  designed  to  be 
used  together  and  are  synchronized  for  maximum  performance.  Synchronization  issues  for  the  case  of  a  non- 
cooperative  bistatic  radar  receiver,  where  the  receiver  attempts  to  “hitchhike”  off  of  an  independent  source  of 
illumination,  is  discussed  in  Section  2.2.  This  mode  of  operation  is  termed  “non-cooperative”  because  the 
transmitter  is  not  designed  to  support  bistatic  radar  operation. 


2.2  Synchronization  Issues  for  a  Non-Cooperative  Bistatic  Radar  Receiver  [2,3] 

In  order  to  maximize  the  radar  coverage  and  detection  sensitivity  of  a  non-cooperative  bistatic  radar 
receiver,  the  receiver  must  either  know  or  be  able  to  estimate  the  following  transmitter  parameters: 

1)  Transmitter  location, 

2)  Antenna  polarization, 

3)  Carrier  frequency, 

4)  Antenna  beam  pattern  and  beamwidth, 

5)  Antenna  scan  characteristics, 

6)  Complex  envelope  waveform. 

Each  of  these  requirements  is  briefly  discussed  in  this  section. 
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2.2.1  Transmitter  Location 


As  pointed  out  in  Section  2.1,  solution  of  the  bistatic  triangle  requires  knowledge  of  the  transmitter 
position.  This  can  be  determined  in  a  variety  of  ways.  Although  non-cooperative,  the  transmitter’s  location  may  be 
known  a  priori.  For  instance,  the  position  of  airport  surveillance  radars  may  be  found  in  various  data  bases.  When 
not  known,  the  transmitter’s  location  can  be  estimated  using  emitter  location  techniques.  Such  techniques  typically 
involve  multiple  receiver  sites  or  a  moving  platform  and  are  generally  good  enough  so  as  not  to  significantly  affect 
target  location  accuracies. 

2.2.2  Antenna  Polarization 

The  antenna  polarization  may  be  known  a  priori  in  some  cases.  Otherwise,  the  polarization  can  be 
determined  by  measuring  the  direct  path  signal  using  antennas  matched  to  various  polarizations.  When  this  is  not 
possible,  a  45°  slant  polarized  antenna  may  be  used  as  a  compromise  by  the  bistatic  radar  receiver. 

2.2.3  Carrier  Frequency 

Waveform  power  can  be  maximized  by  tuning  the  bistatic  radar  receiver  to  the  known  or  estimated  carrier 
frequency  of  the  transmitted  waveform.  Techniques  for  determining  the  carrier  frequency  of  an  unknown  waveform 
are  well  known  and  have  been  implemented  in  electronic  support  measures  (ESM)  receivers.  In  fact,  frequency 
estimation  can  be  done  with  a  higher  accuracy  than  is  required  for  good  bistatic  radar  operation.  It  can,  also,  be 
done  very  rapidly.  For  example,  frequency  can  be  measured  on  a  pulse-to-pulse  basis  by  a  technique  known  as 
instantaneous  frequency  measurement  (IFM).  This  becomes  essential  when  hosting  off  radars  that  are  pulse-to-pulse 
frequency  agile. 

2.2.4  Antenna  Beam  Pattern  and  Beamwidth 

It  is  desirable  to  integrate  all  of  the  pulses  that  the  mainbeam  places  on  the  target.  Knowledge  of  the 
beamwidth  is  necessary  to  set  the  dwell  time  or,  equivalently,  the  coherent  processing  interval  (CPI).  The 
transmitted  pulse  amplitudes  are  modulated  by  the  scanning  of  the  transmitter  beam  for  long  dwells.  Reasonable 
knowledge  of  the  transmitting  antenna  beam  pattern  becomes  necessary  to  do  a  good  job  of  matched  filtering. 
Beamwidth  is  another  parameter  that  can  be  measured  by  ESM  systems.  This  is  generally  derived  by  monitoring  the 
envelope  of  the  transmitted  pulses  over  a  scan.  It  appears  that  an  accuracy  within  25%  is  sufficient  to  maintain  good 
bistatic  radar  receiver  performance. 

2.2.5  Antenna  Scan  Characteristics 

In  order  for  the  antenna  of  the  bistatic  radar  receiver  to  be  pointing  at  the  volume  of  space  simultaneously 
illuminated  by  the  transmitter,  it  is  necessary  to  know  the  transmitter  scan  characteristics.  This  is  referred  to  as 
pulse  chasing.  For  mechanically-scanned  radars,  rotating  at  a  constant  angular  rate,  the  scan  rate  of  the  transmitter  is 
estimated  by  tracking  the  times  when  the  mainbeam  flashes  past  the  receiver.  Knowledge  of  the  scan  time  enables 

estimates  to  be  made  of  the  azimuth  angle  toward  which  the  transmitter  antenna  is  pointing  at  any  time  instant. 

Generally,  it  is  not  possible  to  synchronize  a  bistatic  radar  receiver  to  the  elevation  scanning  of  a  transmitter  unless  a 
known  regular  elevation  scan  pattern  is  used.  To  overcome  this  lack  of  knowledge,  the  bistatic  radar  receiver 
typically  uses  a  receive  aperture  with  an  elevation  fan  beam  which  covers  all  elevations  of  interest. 

2.2.6  Complex  Envelope  Waveform 

The  complex  envelope  of  the  transmitted  waveform,  which  contains  the  amplitude  and  phase  modulation  of 
the  radar  pulse,  must  be  estimated  for  use  in  cross  correlation  with  the  radar  signal.  In  addition,  approximations  to 
the  pulsewidth  and  the  timing  of  pulses,  which  may  or  may  not  be  staggered,  are  needed  to  successfully  perform  the 
operations  of  Doppler  processing,  pulse  compression,  and  matched  filtering. 

As  pointed  out  in  Chapter  1,  one  of  the  objectives  of  this  effort  is  to  use  the  interferometer  signal  to  disrupt 
to  the  extent  possible  the  synchronization  of  the  non-cooperative  bistatic  radar  receiver  with  the  host  radar. 
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2.3.  Coherent  Reference  Denial 

In  this  report  it  is  assumed  that  the  non-cooperative  bistatic  radar  receiver  performs  detection  by  utilizing 
an  estimate  of  the  host  radar  transmitted  signal  for  correlation  with  the  received  data.  In  practice,  the  estimate 
should  be  continuously  updated  in  order  to  account  for  pulse-to-pulse  amplitude,  frequency,  and  phase  variations 
that  may  or  may  not  be  intentional. 

The  estimate  of  the  host  radar  transmitted  complex  envelope  is  typically  obtained  from  the  direct  path 
signal  observed  by  the  bistatic  radar  receiver.  In  general,  the  direct  path  signal  is  corrupted  by  clutter  and  receiver 
noise  which,  in  some  cases,  can  be  sufficient  to  cause  significant  signal  processing  losses  due  to  poor  estimates  of 
the  host  radar  complex  envelope.  However,  for  our  purposes,  it  will  be  assumed  that  the  interferometer  signal  is 
mainly  responsible  for  degradation  of  the  complex  envelope  estimate.  In  this  way,  there  is  no  reliance  solely  on  the 
clutter  and  receiver  noise  to  prevent  the  non-cooperative  bistatic  radar  receiver  from  hosting  off  of  the  monostatic 
radar.  Of  course,  to  the  extent  that  the  clutter  and  receiver  noise  help  to  mask  the  transmitted  radar  waveform  in  the 
direct  path,  the  more  difficult  it  will  be  for  the  non-cooperative  bistatic  radar  receiver  to  operate  successfully. 

With  reference  to  Figure  1.2,  let  the  transmitted  host  radar  and  interferometer  complex  envelopes  be 
denoted  by  sRAd  (t)  and  sIFM(t),  respectively.  By  definition,  the  energies  of  these  signals  are  given  by 


=  f 

Jo 

1  SRAD  (0  1  dt 

(2.19) 

=  f T* 

Jo 

1  1  dt 

(2.20) 

where  both  signals  are  assumed  to  be  of  the  same  duration,  Td.  To  account  explicitly  for  the  signal  energies,  it  is 
convenient  to  introduce  the  unit  energy  signals  fRAD(t)  and  (t)  such  that 


The  complex  envelopes  can  then  be  expressed  as 

srad(0  ~  ^ ^ rad /kad(0 
and 

SIFAf(  0  =  IFMflFM (0- 


(2.21) 


(2.22) 


(2.23) 


For  simplicity,  assume  that  the  direct  path  signal  between  the  monostatic  and  bistatic  radars  is  used  as  the 
coherent  reference  signal  for  the  correlator  in  the  non-cooperative  bistatic  radar  receiver.  With  reference  to  Figure 

2.2,  the  direct  path  consists  of  the  bistatic  baseline  which  is  of  length,  L.  Let  denote  the  antenna  gain  of  the 

monostatic  radar  in  the  direction  of  the  direct  path.  Consequently,  the  power  density  of  the  monostatic  radar  signal 
at  the  bistatic  radar  is 


WDP  — 
n  RAD  ~ 


(E^/TJGZ 


(2.24) 


Let  Gm  denote  the  antenna  gain  of  the  bistatic  radar  in  the  direction  of  the  direct  path.  Then  the  effective  aperture 
area  of  the  bistatic  radar  antenna  in  the  direction  of  the  direct  path  is 


/-<DP  12 


(2.25) 
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where  X  is  the  signal  wavelength.  It  follows  that  the  monostatic  radar  signal  power  intercepted  by  the  bistatic  radar 


pDP  _  WDP  ADP  _ 
rRAD  ~  r’  RAD  A BR  ~ 


77  /~iDP  s~yDP  j}2 

P-A  D  A  V\  vJj/D  VJnn  A 


AD  ^  MR  kjBR‘ 

m2L% 


(2.26) 


Similarly,  let  G^M  denote  the  antenna  gain  of  the  interferometer  in  the  direction  of  the  direct  path.  The 
interferometer  signal  power  intercepted  by  the  bistatic  radar  is  then  given  by 

P  piDP  (~*DP  <2. 

pDP  _  r'IFM(JIFMyJBR  "  n  27) 

Im  {Anfl}Td 

DP 

Assume  the  direct  path  signal  intercepted  by  the  bistatic  radar  receiver  is  delayed  by  an  amount,  T  ,  and  Doppler 
shifted  by  an  amount,  fDP.  If  nDP  (t)  denotes  the  complex  envelope  of  the  receiver  noise  in  the  direct  path,  the 
complex  envelope  of  the  direct  path  signal  may  be  modeled  as 

SD'«)  =  .  (2.28) 

+C,  /™  0-^')  «  V’P'o-r'+n”  (0 


where  k^D  and  k^M  are  complex  constants  such  that 


pDP  rp  _ 

rRAD  Jd  ~ 


77  /^iDPpiDP 
^RAD^MR^BR  A 

(An)2!? 


(2.29) 


p  pi  DP  fiDP  2  2 

I  l2=  PZD  =  im^2L2BR  •  (2-3°) 

With  reference  to  Figure  2.2,  the  signal  reflected  from  the  target  to  the  non-cooperative  bistatic  radar 
receiver  travels  from  the  transmitter  along  the  two  legs  of  the  triangle  whose  lengths  are  RT  and  Rr.  In  general,  this 
signal  consists  of  two  components,  one  from  the  monostatic  radar  and  one  from  the  interferometer.  The  two  legs  of 

the  triangle  are  referred  to  in  this  report  as  the  total  path.  Let  Gm,GIFM,GBR  denote  the  antenna  gains  of  the 
monostatic  radar,  interferometer,  and  bistatic  radar,  respectively,  in  the  target  direction  along  the  total  path.  Also,  let 
<JT  be  the  effective  target  radar  cross  section.  Since  the  power  density  of  the  monostatic  radar  signal  at  the  target  is 


pi jTP  _  (Brad  f_ D)^mr 


(2.31) 


the  power  density  of  the  monostatic  radar  signal  at  the  bistatic  radar  is 

Ur  TP  —  RAD  L Td)Gj 


W  = 

"  RAD 


(2.32) 


Hence,  the  target  signal  power  intercepted  by  the  bistatic  radar  due  to  the  monostatic  radar  is 


(2.33) 


p  piTP  piTP  _  j|2 

pTP  _  ^ RAD  ^MR  ^ BR  &T  A 

RAD~  (4nfR}R2RTd 

Similarly,  the  target  signal  power  intercepted  by  the  bistatic  radar  due  to  the  interferometer  is 

P  s-i  TP  /~i  TP  —  ^  2 

pTP  _  ^  IF M  U  IF M  ^  BR  °  T  A 

,FM  “  (4 x)'R}RX  '  <2M> 

Assinne  the  total  path  signal  received  by  the  bistatic  radar  receiver  is  delayed  by  an  amount,  tTP ,  and  Doppler 

TP  TP 

shifted  by  an  amount,  J  .If  n  ( t )  denotes  the  complex  envelope  of  the  receiver  noise  in  the  total  path,  the 
complex  envelope  of  the  total  path  signal  may  be  modeled  as 


i  TP  ?  TP 

where  kRAD  and  kIFU  are  complex  constants  such  that 


p  /~<TP  s~iTP  o2 

I  kTP  i2 _  pTP  T  _  ^RAD^m^BR  °T  A 

\  RAO  \  RAD  d  (4 JtfR^Rl 


Ti  /~*TP  s-iTP  _  o2 

UTP  |2 _  pTP  rn  _  £'JFMUIFMijBR  °TA 

KIFM  I  '  rIFM1d  ~ 


(An)  Rj.Rr 


The  bistatic  delay  is  defined  in  Equation  (2.3)  as 


_  _  ~TP  ~DP 
Tb=T  -T  , 


Similarly,  the  bistatic  Doppler  shift  is  defined  to  be 


fs  =fTP  ~fDP  ■  (2.38) 

f  TP  ) 

Ideally,  the  non-cooperative  bistatic  receiver  would  like  to  correlate  fRAD{t  —  TTF)eJ2n  with  the  total 

path  signal,  sTP  (t).  To  accomplish  this,  it  would  modify  the  direct  path  signal  by  a  delay,  T B ,  and  a  Doppler  shift, 
fB.  However,  in  practice,  T B  and  fB  are  unknown.  Hence,  the  bistatic  radar  receiver  would  utilize  estimates  of 
Tb  and  fB  .  Let  these  be  denoted  by 


^  ^  -  T 


(2.39) 


A  =fa-V- 


(2.40) 
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To  detect  a  possible  target,  the  bistatic  radar  receiver  would  then  correlate  SDP(t  —  fB)eJ~^B<'  T°  ]  with 
STP  (t).  This  yields  the  test  statistic  employed  for  the  detection  of  a  target. 


By  definition,  correlation  of  the  complex  waveforms  /  (t)  and  f2  (t)  is  given  by 


y  =  J/i  ’(0/2(0  dt 


(2.41) 


where  the  asterisk  denotes  the  operation  of  complex  conjugation.  Hence,  the  test  statistic  at  the  output  of  the 
correlator  in  the  non-cooperative  bistatic  radar  receiver  is 

/  =  J  [  sDP (t-fB) ]*e  -J2'?*-1*-**  sTP  (0  dt.  (2.42) 

Substituting  the  expressions  for  sDP  (t)  and  sTP  (t)  given  by  Equations  (2.28)  and  (2.35)  into  Equation  (2.42),  the 
detection  statistic  becomes 


/=  J I  (*£S,  )*  flc  (I  -  *Dr -t.)e 


e-,u A<"'  ■  +  <0 * /;„  (/  -  *°r -  f, ) 

-  j  U  f°F (t-rDP  -fB)  Dp 

e  e  -  j**  SbO-t  -tB) 

+  (nDP  ( t  -iB)  )*e~J  2*h  0  “  T°P  ~ te)  ) 

+  kppM  fIFM(t-rTP)eJ2*fTP<-'TF) 


+  nTP  (01  dt. 


(2.43) 


Carrying  out  the  product  in  the  integrand  of  / ,  it  is  seen  that  /  is  a  summation  composed  of  nine  integrals.  If  the  i,h 
integral  is  denoted  by  Ij,  the  detection  statistic  can  be  expressed  as 

/=£/,-  (2.44) 

1=1 


where  only  four  of  the  terms  do  not  involve  receiver  noise.  Focusing  on  these  four  terms,  it  is  convenient  to  define 
the  complex  constants 


k  _  (fcDP  \  *  lJP 
'h  —  k^radJ  /v rad 

*2=(0**SL 

*3  =  (0**£d  (2-45) 


k  _  (h-ur  \  *  UP 
"’A  \^1FM  >  "'IFM  • 


Referring  to  Equations  (2.3)  and  (2.39),  note  that 


t~tDP  -fB=t 


-TP  .  ^ 
-T  +  T 


(2.46) 
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Also,  with  reference  to  Equations  (2.38)  and  (2.40), 

/"+/,  =fP-V. 

Making  the  change  of  variable, 


G  -  t  -  T +  T , 


it  can  be  shown  that 


I\-kxe  ~jl*  {TP '  J  fMD  (G)  f^r,  (G-T)  e  ilttvadG. 

— oo 

By  definition,  the  delay  -  Doppler  ambiguity  function  of  the  host  radar  waveform  is 

X  (J>V)  =  —1—  ]  s'jm  (g)srad  ( G-T )  e  j2*vadG 
Erad  -t 

=  J  f*RAD  0 VWrAD  {<*-*)  e  iltCVadG 

— oo 

where  use  was  made  of  Equations  (2.22)  and  (2.23).  As  a  result,  Ii  may  be  expressed  as 

h=kxe  -]lnfTPx  x  (J,v)- 

In  a  similar  manner,  it  can  be  shown  that  I2  and  I3  are  given  by 

I2=k2e-j2*/TFr  J  (g-t)  ej2nvadG 


(2.49) 


(2.50) 


I3=k3  e  -J2nfTP’  J  fjFM  (a)  (G-T)  ej2Kvada.  (2.53) 

By  definition,  the  delay  -  Doppler  cross  ambiguity  function  between  the  interferometer  and  host  radar  waveform  is 


Z  (r,  v)  . . . —  j  s'jFM(g)Srad  (G-T)ei2*vada 

yEIFM  t,MD 

=  J  flFM  (e)fRAD  (<T-T)e  J2nVCdG. 


Note  that 


Z  (~T,-V)=  f  fIFM(G)fRAD(G +T)e-j2nvodG 


Z  *  (-T,  -  V)  =  J  f1FM  (g)  fxAo  ( G+T)e  J2Kva{dG. 


(2.56) 
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With  the  change  of  variable,  OC  =  O  +  T  ,  Z  *  (— T,  —  v)  becomes 

Z  *  (-r,  - v)  =  ]  f'MD  (a)fWM  (a  -r)  e  j2*v{a~T)doc 

(2.57) 

=  e  -JlnvT  ]  J' rad  (a)/™  (a -r)ej2n  vada. 

It  follows  that  I2  and  I3  can  be  expressed  as 

I2=k 2  e  ~j2n{fTP  ~v)t  t  (-T,  -  V)  (2.58) 

and 

I3=k3  e  ~j2nfT"r  Z(t,V).  (2.59) 

Finally,  it  can  be  shown  that 

It=k,e-lu/"’  ]/;„(<?)  (a -T)  eJlm  da.  (2.60) 

By  definition,  the  delay  -  Doppler  ambiguity  function  of  the  interferometer  waveform  is 

T  (/,  v)  =  ~  J  s'ifm(<j)  sIFM  (a  -t)  eJ2™a  da 

hlFM  -«  (2.61) 

=  j/™(0O  fiFM  (O’-'T)  ej2m<T  do. 

Consequently, 

I^kAej2KfTPx  r(T,  V).  (2.62) 

In  summary,  the  detection  statistic  used  by  the  non-cooperative  bistatic  radar  receiver  is 

l  =  e-J'-r”'[kt  x(ZtV)+k,e‘M’ t  (-T.-V)  (263) 

+  k3  Z  (t,  v  )  +  k4  r  (t,  t»)]  +  [terms  that  involve  receiver  noise]. 

The  desired  signal  component  of  the  detection  statistic  is 

ls=k j  e~j2nfTP  x  X  (T=  V).  (2.64) 

The  remaining  terms  represent  noise  and  interference.  In  general,  the  desired  signal  component  is  corrupted  by 
thermal  noise,  clutter,  multipath  and/or  other  propagation  effects.  However,  it  is  the  intention  of  this  effort  to  design 
the  interferometer  signal  such  that  it  effectively  masks  the  direct  path  signal  from  the  host  radar.  The  expression  for 
the  detection  statistic  in  Equation  (2.63)  clearly  shows  that  the  objective  of  coherent  reference  denial  depends  upon 
the  delay-Doppler  ambiguity  functions  of  the  host  radar  and  interferometer  waveforms  and  their  cross  ambiguity 
function. 
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CHAPTER  3 


INTERFEROMETER  AND  HOST  RADAR  ANTENNA  PATTERNS 


3.1  Interferometer  Antenna  Pattern 

In  addition  to  coherent  reference  denial,  the  interferometer  signal  is  to  be  used  for  communications  in  many 
different  directions.  Because  the  interferometer  will  serve  as  a  transmitter,  its  antenna  pattern  is  derived  in  this 
section  in  the  context  of  a  radiation  pattern.  However,  the  radiation  and  reception  patterns  of  an  antenna  are  known 
to  be  identical. 

Consider  a  point,  P0,  in  the  far  field  of  the  interferometer,  as  shown  in  Figure  3.1.  The  cartesian 
coordinates  of  P0  are  given  by  {%0 ,yo,zo)  while  its  spherical  coordinates  are  (ro ,6o,<po).  These  are  related  by 

^sin^eos^ 

y0  -  ro  s*n  Q0  sin  (f>0 .  (3.1) 

ze  =  ro  cos#a 


Figure  3.1 :  Both  Cartesian  and  spherical  coordinates  are  used  to  characterize  a  point  in  the  far  field  of  the  interferometer. 


Assume  each  antenna  in  the  interferometer  pair  consists  of  a  dipole  antenna  of  length,  /  ,  as  shown  in 
Figure  3.2.  Note  that  each  leg  is  of  length,  111.  Let  P  denote  an  arbitrary  point  in  the  far  field  of  the  antenna  with 
spherical  coordinates  {r,6,(j)).  Assuming  a  sinusoidal  excitation  at  frequency,  (00 ,  and  using  a  far  field 
approximation  [4],  the  electric  field  at  P  is  approximately  given  by  the  real  part  of 

E  =  —F(e)e-",-'e”’-'  (3.2) 

r 


where  K  is  a  constant, 


Po  -mo  -  27tfo 


(3.3) 


A0  is  the  signal  wavelength,  and 
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F(e)= 


COS 

■  cos  9 

-cos 

1 _ 

cos 

id  Q 
—  cos  9 

-cos 

— 

la  < 

_ 

L  2 

L  2  J 

_A> 

lAo] 

sin  9  sin  6 


(3.4) 


Figure  3.2:  Each  antenna  in  the  interferometer  pair  consists  of  a  dipole  antenna  of  length,  /  . 


Consistent  with  the  circular  symmetry  of  the  dipole,  the  field  is  seen  to  be  independent  of  (p .  The  most  frequently 
encountered  case  is  the  half-wave  dipole  for  which  the  total  length  of  the  dipole  antenna  is  /  =  /  2 .  Then 


F(e)= 


cos 

—  cosy 

2 

sin  6 


(3.5) 


Let  the  two  dipole  antennas  of  the  interferometer  pair  be  directed  in  the  z-direction  and  placed  on  the  x-axis 
at  and  X  =  +  d IFM  12,  respectively,  as  shown  in  Figure  3.3.  Thus,  diFM  is  the  spacing 

between  the  interferometer  elements.  Consider  a  point,  P,  in  the  near  field,  as  illustrated  in  Figure  3.4.  The  lengths 
of  lines  OP,  BP,  and  AP  are  denoted  by  r,  ri,  and  r2,  respectively.  The  angles  from  the  z-direction  to  lines  OP,  BP, 

and  AP  are  given  by  9, 9l ,  and92 ,  respectively.  Finally  y,  yl  ,  and  y2  ,  respectfully,  denote  the  angles  between 
the  x-axis  and  the  lines  OP,  BP,  and  AP.  As  point,  P,  moves  into  the  far  field  of  the  antennas,  the  following 
approximations  become  valid: 


2  2 

Figure  33:  The  dipole  antennas  are  placed  along  the  x-axis  at  A  and  B  and  are  directed  in  the  z-direction. 
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Figure  3.4:  When  P  is  in  the  near  field  of  the  interferometer,  the  geometry  is  relatively  complicated. 


r{~r2~r 

0{~92~  6  (3.6) 

7i  =  y2  ~r- 

The  far-fteld  geometry  is  shown  in  Figure  3.5. 

Assuming  both  antenna  elements  of  the  interferometer  pair  are  excited  by  the  same  sinusoidal  signal  at 
frequency,  6)g  ,  the  total  electric  field  at  point  P  when  P  is  in  the  near  field  is  given  by  the  real  part  of 


Figure  3.5:  The  geometty  simplifies  as  P  moves  to  the  far  field 
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However,  the  phase  angles, 


F{ox)=F(e2)~F{e). 


Por\  = 


Por 2  =  2tf 


are  very  sensitive  to  small  changes  in  X\  and  r2  on  the  order  of  a  wavelength  because  of  their  dependence  on  the 
electrical  distances  r,  /  Xg  and  r2/  X0.  From  Figure  3.5, 

rx  ~r  +Ar 

r2  ~  r  —  Ar .  (3.12) 

Using  Equations  (3.8),  (3.9),  and  (3.12)  in  Equation  (3.7),  the  total  electric  field  at  point,  P,  becomes 

E TOT  -  a  F(ey  Me-  \JmJ 

IFM  r  [ 

2K  r«  (k  )i  -iP°r  j<°ot  ..... 

= - F(#jcos  \p0  {Ar)\e  e  .  (3.13) 


From  Figure  3.5, 


A r  =  ^IFM  cos  y 
2 


where  7  is  the  angle  between  ,  the  unit  vector  in  the  r-direction  of  spherical  coordinates,  and  ^  ,  the  unit  vector 
in  the  x-direction  of  Cartesian  coordinates.  Hence, 

COS  y  =7  (3.15) 

r  x 

where  the  dot  indicates  a  dot  product  between  the  two  vectors.  Thus,  COS  Y  is  the  projection  of  7  on^  .  From 
Equation  (3.1),  it  follows  that 

cos  7  =  sin#  cos  0.  (3.16) 

Therefore,  utilizing  Equations  (3.14)  and  (3.16)  in  Equation  (3.13),  the  total  electric  field  at  point,  P,  is  expressed  as 

TOT  _  2 K  „J Pn  J  J  _  -jP0r  .  M* 

£ IFM 


F($)cos  ^-djpu  sin 0 cos ^  e  e' 
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2  K 


f  d ir>/  ^ 


=  —  F(0)cos  n\2®L 
r  L  v  to  J 

The  interferometer  array  factor  is  defined  to  be 


sin  9  cos  (j) 


_ jPor  ja>ot 

e  e 


(3.17) 


~  cos 


'  (d  \ 

n  -J^L  sin#cos0 

.  V  K  ) 


(3.18) 


Hence, 


rTOT  r?(Q\r  (Q  ,\  -  jfior  j®0t 

E 1FM  ~  F\@riFM  \@-><P) e  e 

r 


(3.19) 


Observe  that  the  resultant  field  involves  the  product  of  the  pattern  of  the  individual  (identical)  antenna  elements  and 
the  array  factor  which  is  a  function  of  the  electrical  distance,  dIFM  /  Aa  .  This  is  referred  to  as  the  principle  of 

pattern  multiplication.  If  antennas  other  than  dipoles  are  used  for  the  array  elements,  is  obtained  by  merely 
substituting  the  appropriate  F{0)  for  the  antenna  elements. 

The  above  analysis  makes  use  of  traditional  spherical  coordinates  in  which  the  azimuth  angle,  <p ,  is 
measured  from  the  x-axis.  However,  in  the  conventional  array  literature,  the  azimuth  angle,  a  ,  is  measured  from 
the  perpendicular  to  a  linear  array.  When  the  linear  array  is  placed  along  the  x-axis,  a  is  measured  from  the  y-axis, 
as  illustrated  in  Figure  3.6.  Note  that 

cos  </)  —  sin  cc .  (3.20) 


Figure  3.6:  When  a  linear  array  is  placed  on  the  x-axis,  the  azimuth  angle,  (X  ,  in  the  array  literature  is  typically  measured  from 

the  y-axis. 


In  terms  of  CC ,  the  total  electric  field  at  point,  P  is 

2  K 


fTOT  ~ 
^  IF  M  ~  ' 


(3.21) 


where 


i 
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Fifm  °0  ~  cos 


K 


\ 


' IFM 


v  K  j 


sin  6  sin  a 


(3.22) 


Radiation  patterns  are  conventionally  plotted  in  decibels  where  the  pattern  is  normalized  such  that  the 
maximum  value  is  OdB.  The  radiation  pattern  of  the  dipole  antenna  used  for  the  interferometer  elements  is 
considered  first.  In  decibels,  this  is  defined  to  be 


c(e)=ioiog]0 


F(e) 

He}m 


(3.23) 


where  max  denotes  the  maximum  value  of  F{6)  and  F(9)\s  given  by  Equation  (3.4).  Detailed  plots  of 
G(0)  are  presented  in  Figures  3.7,  3.8,  and  3.9  for  various  antenna  lengths  [5]. 

The  polar  plots  in  Figure  3.7  pertain  to  situations  where  l  <  Aa  .  It  is  seen  from  Figure  3.7  that  the 

radiation  pattern  for  this  case  includes  only  one  major  lobe  symmetrically  placed  around  the  antenna  axis.  Note  that 
the  beam  widths  of  the  main  lobes  decrease  with  an  increase  in  the  antenna  length.  Thus,  the  longer  the  antenna  is, 
the  more  directive  is  the  radiated  energy. 


o+\+o 

0° 


180° 


l«Xo 
l  -V  4 
l  -V  2 
/  =3V4 
/ 


Figure  3.7:  Radiation  patterns  are  plotted  for  the  dipole  antenna  with  /  <  . 


21 


Figure  3.8:  The  dipole  antenna  radiation  pattern  is  plotted  for  /  =  1 .25 Ag . 

Figure  3.8  shows  the  dipole  element  radiation  pattern  for  /  =  1 .25 Ag  Observe  that  there  are  now  three 

major  lobes  symmetrically  placed  around  the  antenna  axis.  The  lobe  with  maximum  gain  is  centered  at  9  —  90° . 
The  dipole  element  radiation  pattern  for  l  =  1 -5A0  is  presented  in  Figure  3.9.  As  in  Figure  3.8,  there  are  three 
major  lobes  symmetrically  placed  around  the  antenna  axis.  However,  now  there  are  two  lobes  of  equal  maximum 

intensity  located  at  6  =  40°  and  0  =  140°  .  With  a  still  further  increase  in  the  electrical  length  of  the  antenna,  a 
larger  number  of  lobes  begin  to  appear  in  the  element  radiation  pattern.  The  radiation  pattern  of  the  dipole  antenna 
is  seen  to  lose  its  directional  properties  for  antenna  lengths  larger  than  Ao  . 

The  radiation  pattern  of  the  interferometer  array  factor,  given  in  Equation  (3.22),  is  considered  next.  As 
indicated  in  Figure  1 .2,  it  is  first  assumed  that  the  interferometer  elements  are  placed  on  a  linear  grid  about  the  main 
radar  elements.  The  configuration  is  illustrated  in  Figure  3.10  where  the  spacing  between  antenna  elements  is 
denoted  by  d.  Independent  of  whether  the  N-element  main  radar  antenna  array  contains  an  odd  or  even  number  of 
elements,  note  that  the  spacing  between  the  interferometer  element  pair  is 


Figure  3.9:  The  dipole  antenna  radiation  pattern  is  shown  for  /  =  1 .5  Ac  . 


^ifm  ~  I 


d  _ 

(N  +  l)d] 

l  2  J 

Cl 

2  r\ 

=  (N  +  l)d. 


(3.24) 
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Figure  3.10:  The  antenna  configuration  consists  of  a  linear  grid  with  equal  spacing,  d.  The  main  radar  array  is  composed  of  N 

elements  where  N  can  be  either  an  (a)  odd  or  (b)  even  integer. 
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The  maxima  and  zeros  of  \Fifm  {0,  Cf)J  are  seen  to  be  a  function  of  both  6  and  OL .  To  visualize  the  azimuthal 

radiation  pattern  of  the  array  factor,  it  is  convenient  to  hold  6  fixed  and  vary  OL .  Alternatively,  to  visualize  the 
elevation  radiation  pattern,  it  is  convenient  to  hold  OL  fixed  and  vary  0  . 

The  azimuthal  radiation  pattern  is  considered  first.  For  this  purpose,  let  6  equal  the  fixed  angle,  9a .  The 
maxima  are  now  governed  by 


where 


Note  that 


f  d\ 

(iV  +  l);r  —  sin0osina  =  tfj;r 

v/U 


m  =  0,±1,±2, ... . 


sin  a  = 


m  ( 1 


{N  + 1)  v  ^  )  sin  90 


|sin  <  1, 

it  follows  that 

H<(V  +  l)^|-l|sine,|. 

Similarly,  with  6  —  0o,  the  zeros  of  the  array  factor  occur  for 


( d\  n 

(w  +  lk  —  sin^„sina  =  n  — 

UJ  2 


n  =  ±1,±3,±5,... . 


sin  a  = 


n  f  /O  1 


2(iV  +  l)UJsin^' 

Because  |sin  CX |  cannot  exceed  unity,  it  is  concluded  that 


hi  ^  2iN + lsin^l- 


Thus  far,  the  discussion  has  focused  on  the  values  of  OL  for  which  Fifm  { 0o,CL )  attains  either  its  maximum 
or  minimum  values.  The  question  arises  as  to  whether  or  not  the  azimuthal  radiation  pattern  contains  minor  lobes  in 
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addition  to  the  major  lobes.  To  answer  this  question,  the  derivative  of  Fim{90,a )  with  respect  to  Of  is 
considered.  With  reference  to  Equation  (3.25), 

— EM. &£  =  -(n  +  l)rf — 1  sin  $0  sin  +  sin  G0  sines'  cos  or.  (3.38) 

da  yk0)  I  v4J 


The  derivative  equals  zero  when  either 


(  d\ 

sin  (iV  +  l);r  —  sin  sin  or  =0 


(3.39) 


cos  or  =  0 . 


Note  that  Equation  (3.39)  is  satisfied  for 


(d\ 

(N  +  \)m  —  sin  Q0  sin  or  =  mn 


where 


m  —  0,±1,±2, ... . 


This  result  is  identical  to  that  of  Equations  (3.29)  and  (3.30)  and  yields  the  major  lobes  of  the  azimuthal  radiation 
pattern  for  which 


F,n(e„a)  =  ±  1. 


(3.43) 


Equation  (3.40)  is  satisfied  for 


,  n 
or  =  ±  — . 
2 


However,  it  should  be  remembered  that  the  derivative  equal  to  zero  at  a  point  only  guarantees  that  the  slope  of  the 
function  is  zero  at  that  point.  Consequently,  either  a  major  lobe,  minor  lobe,  or  inflection  point  is  possible  at 

a-±— . 

2 

Polar  plots  for  the  azimuthal  radiation  pattern  of  the  interferometer  array  factor  are  presented  in  Figure 
3.1 1,  where  N  is  allowed  to  range  from  N=2  to  N=7,  and 


2  2 


(3.45) 


For  this  special  case, 


dIFM={N  +  l)  -f 


(3.46) 
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Fifm  —,a  =  cos  (7V  +  l)— sinof  . 
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Figure  3.11 :  The  azimuthal  radiation  pattern  of  the  interferometer  array  factor  is  plotted  for  d  —  A0  /  2  ,  #o  =  Jt  /  2 ,  and 
d™  equal  to  (a)  3Ao  /  2  ,  (b)  2Ag  ,  (c)  5/l0  /  2  ,  (d)  3/l0 ,  (e)  7Aa  /  2  ,  and  (f)  4Ao . 

Because  6o  —  •f ,  the  plots  in  Figure  3.11  show  the  horizontal  radiation  pattern  in  the  x-y  plane  at  z=0.  With 
reference  to  Equation  (3.31)  and  (3.33),  the  maxima  occur  for 


sin  Cf  = 


(JV  +  l) 


m  =  0, ±  1, ±  2, •  •  •,  ±  floor\ 


N  + 1 


(3.49) 


and  floor  (•) denotes  the  largest  integer  less  than  or  equal  to  its  argument.  Similarly,  from  Equations  (3.36)  and 
(3.37),  zeros  appear  when 


sin  or  = 


(N  + 1) 


(3.50) 


where 


n  =  ±  1,  ±3,-",  ±ofloor{N  +  l) 

and  ofloor  (•) denotes  the  largest  odd  integer  less  than  or  equal  to  its  argument.  Because 


1  =  C°S  (jV  +  l)^sinor  , 

V.2  J  L  2 


(3.52) 
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the  radiation  pattern  assumes  either  positive  or  negative  real  values.  In  Figure  3.11  the  lobes  corresponding  to 
positive  values  of  the  cosine  function  are  denoted  by  “+”  while  those  corresponding  to  negative  values  of  the  cosine 
function  are  denoted  by 

In  general,  the  azimuthal  radiation  pattern  of  the  interferometer  array  factor  can  be  shown  to  have  the 
following  properties  when  the  conditions  of  Equation  (3.45)  hold: 

1 .  For  all  integer  values  of  N,  either  odd  or  even, 

a.  symmetry  exists  about  both  the  x-  and  y-  axes, 

b.  the  number  of  major  lobes  equals  (2  N  +2), 

c.  there  are  no  minor  lobes, 

d.  major  lobes  always  occur  at  OC  =  0°  and  1 80° 

2.  When  N  is  an  odd  integer, 

a.  the  polarity  of  the  major  lobes  strictly  alternate  in  sign  from  one  lobe  to  the  next, 

b.  major  lobes  always  occur  at  (X  =  ±90°  . 

3.  When  N  is  an  even  integer, 

a.  focusing  entirely  above  or  below  the  x-axis,  the  polarity  of  the  major  lobes  strictly  alternate  in 
sign  from  one  lobe  to  the  next, 

b.  the  first  major  lobes  positioned  on  either  side  of  the  x-axis  have  the  same  polarity, 

c.  nulls  always  occur  at  CC  —  ±90°  . 

These  properties  are  clearly  demonstrated  by  the  plots  in  Figure  3.11. 

The  nulls  appearing  at  (X  —  ±90°  when  N  is  an  even  integer  are  of  particular  interest.  As  pointed  out  by 

the  discussion  leading  up  to  Equation  (3.44),  the  derivative  of  the  array  factor  is  zero  for  a  =  ±90°  .  Therefore, 
the  nulls  at  these  angles  have  zero  slope.  This  is  consistent  with  the  property  that  the  first  major  lobes  positioned  on 
either  side  of  the  x-axis  have  the  same  polarity.  This  is  illustrated  in  Figure  3.12  where  a  cartesian  plot  of 

FmM  {it  12,  Ct),  corresponding  to  N=2,  is  presented  with  a  varying  from  0°  to  180°  .  Observe  that  the  null  at 
(X  =  90° ,  where  the  polarity  of  FWM  ijt  /  2,  a)  is  unchanged  and  the  slope  is  zero,  is  noticeably  broader  than  the 

nulls  at  a  =  19.5°  and  160.5° ,  where  the  polarity  of  F1fm  {jt  /  2,  a)  does  change.  To  avoid  self-jamming  of  the 
radar  and  communications  waveforms,  it  may  be  desirable  to  steer  the  interferometer  pattern  such  that  the  main 
beam  of  the  host  radar  is  centered  in  this  broad  null.  Steering  of  the  interferometer  and  placement  of  the  broad  null 
is  investigated  in  Section  3.3. 


Alpha  in  Degrees 

Figure  3.12:  The  Cartesian  plot  of  Fifu  (jt  /  2,  (X) ,  corresponding  to  N=2,  illustrates  the  broad  null  that  appears  at 

a  =  90°. 
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Additional  insight  into  the  radiation  pattern  of  the  interferometer  array  factor  is  gained  by  considering  the 
case  for  which 

J  ^ a  n  7T 

d  =  -t’eo=T-  (3-53) 

z  o 


|  =  |  COS  (-/v  +  l^siim  . 


Polar  plots  for  this  case  are  presented  in  Figure  3.13  where,  as  with  Figure  3.1 1,  the  number  of  antenna  elements  in 
the  main  radar  are  allowed  to  vary  from  N=2  to  N=7.  Because  60  —  It  /  6 ,  the  plots  of  Figure  3.13  illustrate  the 

radiation  pattern  as  seen  looking  into  a  cone  emerging  from  the  origin  whose  axis  is  oriented  in  the  z-direction  and 
whose  sides  are  at  an  angle  of  30°  from  the  z-axis. 


At  a  specified  angle,  Ct0 ,  the  vertical  pattern,  \Fjfm  (&,  OCo  )|  is  helpful  in  understanding  the  lobing  structure 

of  the  azimuthal  pattern,  |  F1fm  {0o  ,  a) | .  For  example,  consider  0Co  =  ±90°  .  As  before,  let  d  —  Xot  1 .  From 
Equation  (3.24),  the  vertical  radiation  pattern  is  given  by 


^™^±y)|=IC0S  (Ar+1)f sin o  . 


(3.55) 
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Figure  3.13:  The  Azimuthal  radiation  pattern  of  the  interferometer  array  factor  is  plotted  for 
d  =  Xa  /  2 , 0o  =  It  1 6 ,  and  dffM  equal  to  (a)  3X0  /  2 ,  (b)  2Xg  ,  (c)  5\  /  2 ,  (d)  3/l0 ,  (e)  1 X0  /  2  ,  and  (f)  4XQ . 
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This  vertical  pattern  is  plotted  in  Figure  3.14  as  N  varies  from  N  =  2  to  N  =  7.  First,  focus  on  the  patterns  at  the 
coordinates  90  —  7t / 2 ,«  =  ±7t / 2  .  By  comparing  the  plots  in  Figure  3.11  with  those  in  Figure  3.14,  it  is 
important  to  note  that  the  vertical  pattern  experiences  a  broad  null  for  the  very  same  values  of  N  that  yields  a  broad 
null  in  the  horizontal  pattern.  Next,  focus  on  the  patterns  at  the  coordinates  90  =  7t  1 6,CC0  =  ±7T  /  2  .  The  major 

lobes  at  a  =  ±K  /  2  in  Figures  3.13  (b)  and  (f)  are  seen  to  coincide  with  the  major  lobes  at  9  =  30°  in  Figures  3.14 
(b)  and  (f)  for  N  =  3  and  7,  respectively.  Similarly,  the  nulls  at  a  =  ± 7T  /  2  in  Figure  3.13  (d)  are  seen  to  coincide 

with  the  null  at  9  -  30°  in  Figure  3.14  (d)  for  N  =  5.  Finally,  the  values  at  a  =  ±K  12  in  Figures  3.13  (a),  (c), 
and  (e)  are  seen  to  be  due  to  “cuts”  through  the  major  lobes  located  at  9  =  41.8°, 0  =  23.6° and  0  =  34.9°, 
respectively,  in  Figures  3.14  (a),  (c),  and  (e). 

3.2  Main  Radar  Antenna  Pattern 

As  shown  in  Figure  3.10,  the  main  radar  antenna  is  assumed  to  be  composed  of  N  dipole  antenna  elements 
equispaced  by  a  distance,  d,  along  the  x-axis.  Assuming  each  of  the  N  elements  is  excited  by  the  same  sinusoidal 
source  at  frequency,  CO  0,  the  total  electric  field  at  a  point  P  in  the  far  field  with  coordinates,  (r,  9,  a) ,  can  be 
shown  to  be  given  by  [5] 


(d)N=5 


(e)N=6 


(f)N=7 


Figure  3.14:  The  vertical  radiation  pattern  of  the  interferometer  array  factor  is  plotted  for  d  —  Aa  /  2,  <X0  —  +7T  /  2,  and  dIFM 
equal  to  (a)  3A0  /  2 ,  (b)  2Ag  ,  (c)  5Ag  /  2 ,  (d)  3Ag  ,  (e)  7  Ag  /  2  ,  and  (f)  4  Ag  . 


eZ  •*—F(e)FUD(e,a)e 


JP‘reJ‘°°' 


(3.56) 


where,  as  in  Section  3.1,  K  is  a  constant,  /?0  is  defined  in  Equation  (3.3),  and  FRAD(9,a) is  the  main  radar  array 
factor  with 
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/  x  _  1  sin  [Md/K)  sin#  sin  a] 

r RAD\P’a)  —  T7  :  r  /  ,  ,  2  \  .  2T" •  T'  (3.57) 

A'  sin[;r(z///l0  )sm#sinorJ 

Observe  that  the  maximum  value  of  |-F^0  (#,Or)j  is  unity.  Of  course,  with  N=2,  the  above  array  factor  should 
reduce  to  that  for  the  interferometer.  When  N=2, 

1  MMi  )sin0sintfj 


Frad{@’&)  ~ 


2  sink-f sin0sincd 


Making  use  of  the  trigonometric  identity. 


sin  2x  =  2  sin  x  cos  x, 

Equation  (3.58)  simplifies  to 

Frad  (#,#)  =  cos [n(d  /  \ )  sin  6  sin  a] . 
This  is  identical  to  Equation  (3.22)  where  diFM  appears  in  place  of  d. 

Of  particular  interest  is  the  case  for  which 


d  =  2/2. 


Then  the  main  radar  array  factor  becomes 


F  (f)  1  sin[N(;r/2)sinff  since] 

1010  ’  N  sin[(^  /  2 )  sin  #  sin  a] 


For  this  choice  of  d,  it  can  be  shown  that  the  array  factor  given  by  Equation  (3.62)  does  not  contain  any  grating 

lobes.  The  array  factor  has  a  maximum  value  of  unity  which  occurs  for  (X  =  0° ,  180°.  The  zeros  of  the  array  factor 
occur  when  its  numerator  is  zero  but  the  denominator  is  not  simultaneously  zero.  Thus,  zeros  result  when 

N— sin#  sin  a  =  nn:  (3.63) 


n  =  ±1,±2,±3,---. 

Since  CC  =  0°  corresponds  to  a  maximum,  note  that  n  =  0  is  excluded.  At  a  null 


sin  or  = 


TV  sin# 


Because  sin  Or|  cannot  exceed  unity,  it  follows  that 


uVsin# 


Consequently,  at  a  null,  n  can  be  any  integer  which  is  restricted  to  the  range. 


0  <  In  < 


Nsin6 

2 


The  sidelobes  of  the  main  radar  array  factor  are  now  investigated.  Assume  N  is  large  enough  such  that  the 
denominator  of  Equation  (3.62)  varies  slowly  relative  to  fluctuations  in  the  numerator.  The  sidelobe  peaks  then 
correspond  to  maxima  of  the  numerator.  These  result  when 

7t  WiK 

N — sin  #  sin  a  = -  (3.68) 

2  2 

where  m  is  an  odd  integer.  Equivalently,  sidelobe  peaks  occur  when 


sma  = 


TV  sin# 

where  m  is  an  odd  integer.  Note  that  m  =  ±1  in  Equation  (3.69)  yields 


a  =  sin  1  ± 


L  Wsin#J 

However,  according  to  Equation  (3.65),  the  first  nulls  on  either  side  of  the  main  beam  at  (X  =  0°  occur  for 


•  -iL  2 

or  =  sin  ± -  . 

Ns\n6 


(3.71) 


iin  1  + — - —  <  sin-1  ± 
NsintfJ 


2 

N  sin# 


(3.72) 


the  sidelobe  peaks  corresponding  to  m  =  ±1  fall  within  the  main  beam  centered  at  OT  =  0° .  Hence,  m  -  ±1  are 
excluded  from  the  allowable  values  of  m.  Because  |sin  Qr|  cannot  exceed  unity,  note  from  Equation  (3.69)  that 


m\  <  |iV  sin  #| . 


(3.73) 


In  summary,  it  is  concluded  that  the  sidelobe  peaks  occur  for  integer  values  of  m  such  that  Equation  (3.69)  is 
satisfied  where 

l<|w?|<  |iVsin  #|.  (3.74) 

The  width  of  the  main  beam  at  OT  =  0°  is  of  interest.  Let  the  beamwidth  be  measured  between  the 
first  nulls  on  either  side  of  the  main  beam.  Denote  the  angles  of  these  nulls  by  (X  .. .  From  Equation  (3.65), 


sink •-„„«)= ± 


2 

N  sin# 


(3.75) 


For  N  sufficiently  large  relative  to  sin  # , 


a\s,-null  <<L 


(3.76) 
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Then 


It  follows  that 


The  beamwidth  is  then  given  by 


a.s,  «± - . 

l-nul1  N  sin  9 


BW  =2  a 

DYVV -nulls  Z 


A^Jsin  0\ 


(3.79) 


In  some  cases,  the  beamwidth  is  measured  between  the  3-dB  points  of  the  main  beam.  Let  the  angles  at 
which  the  3-dB  points  occur  be  denoted  by  Cly_dB  .  By  definition,  these  are  the  angles  for  which  the  magnitude  of 

the  array  factor  is  equal  to  1  /  yfl  =  0.707.  It  follows  that 


1  sin[Af  f  sin  9  sin  (a3_dB )] 
N  sin[f  sintfsin^,^)] 


=  0.707. 


For  suitable  choices  of  9  and  N 


a3-dB  «  S*n  #1^1. 


sin(a3d9 )  ~  a3_d 


sin  ~  sin  9  sin(a3_dB )  =  a3_dB  ^  sin  9. 


With  these  approximations,  Equation  (3.80)  simplifies  to 


Mir 

sin#  =  0.707 a3_dB  — sin  9 . 


(3.83) 


(3.84) 


This  is  a  transcendental  equation  which  can  be  solved  either  graphically  or  numerically.  The  horizontal  radiation 
pattern  is  obtained  when  9  —K 1 2  .  Equation  (3.84)  becomes 


•  "  N7C~\ 
an  «3 -dB  —  “  °- 


707a, 


The  solution  to  this  equation  is 


«3  -dB  =  ± 


0.886 

N 


(3.85) 


(3.86) 


as  can  be  verified  through  substitution.  It  follows  that 
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BW} _dB  -  2|«3 


1.772 


for  the  main  beam  of  the  horizontal  radiation  pattern. 


The  following  properties  of  the  main  radar  array  factor  hold  as  N  increases: 

1.  The  main  lobe  narrows, 

2.  The  number  of  sidelobes  increase, 

3.  The  peak  magnitudes  of  the  sidelobes  decrease. 

By  way  of  example,  the  horizontal  radiation  pattern  is  plotted  in  Figure  3.15  for 

=  i  0  =  =  6. 

2  2 


Analytically,  this  is  given  by 


(n  ^  1  sin  [3^  sin  or] 

^rad\  \  ~  T  •  r_  ■  T 


6  sin[fsinor] 


The  first  nulls  on  either  side  of  the  main  beam  at  OC  =  0°  occur  for 


Figure  3.15:  The  azimuthal  radiation  pattern  of  the  main  radar  array  factor  is  plotted  for  a  =  ,  U  —  — ,  and  N=6. 


Consequently, 


as,  ,=±19.5°. 

f-null 


BWSI  =39° 

r -nulls 


The  approximate  expression  for  this  beamwidth,  as  given  by  Equation  (3.79),  yields 

4  4  2 


V -nulls 


N  6  3 


radians  =  38.2  . 


These  two  results  are  seen  to  be  in  close  agreement.  In  addition,  the  approximate  expression  in  Equation  (3.87)  for 
the  beamwidth,  as  measured  between  the  3-dB  points,  yields 
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As  mentioned  in  Section  3.1,  it  may  be  desirable  to  steer  the  interferometer  pattern  such  that  the  broad  null 
of  the  interferometer  is  centered  on  the  main  beam  of  the  radar.  With  reference  to  Figures  3. 1 1  (e)  and  3.15,  observe 

that  the  broad  null  of  the  interferometer  is  positioned  at  a  =  90°  while  the  main  beam  of  the  radar  is  at  a  =  0° . 
Issues  involved  with  steering  the  interferometer  and  positioning  the  broad  null  are  discussed  in  Section  3.3. 

3.3  Steering  the  Interferometer  Along  With  Positioning  the  Broad  Null 

As  shown  in  Figure  3.11,  a  broad  null  appears  in  the  horizontal  radiation  pattern  of  the  interferometer  array  factor 
when  d  =  V2,0O  =x!  2,  and  N  is  an  even  integer.  To  minimize  interference  in  target  returns  received  by  the 

main  radar,  which  could  consist  of  both  the  interferometer  and  radar  signals,  it  is  desirable  that  the  broad  null  of  the 
interferometer  be  centered  on  the  main  beam  of  the  radar.  In  addition,  a  steering  capability  should  be  incorporated 
into  the  interferometer  so  that  separate  communications  signals  can  be  directed  to  different  receivers  located  at 
various  angles.  These  considerations  are  treated  in  this  section. 

3.3.1  Steering  the  Interferometer 

The  interferometer  pair  can  be  steered  by  inserting  an  appropriate  phase  shift  into  each  channel.  In  particular,  with 
reference  to  Figure  3.4,  let  the  phase  shift  associated  with  the  element  at  position  B  be  y/l  and  the  phase  shift 
associated  with  the  element  at  position  A  be  \f/2  where 

(3-95) 


The  total  electric  field  at  point,  P,  is  given  by  the  real  part  of 


E]°TM=\f-A9  e^jeM-'  (3, 


where  GJn  is  the  frequency  of  the  sinusoidal  excitation.  Let  P  be  in  the  far  field  such  that  the  approximations  of 
Equations  (3.6),  (3.8),  (3.9),  and  (3.12)  are  valid.  It  follows  that 


■F{6)e~ip°rejVl 


~jP  c(Ar)  +  ejP  o(Ar)e  -■/'(AyOlgM-' 


However,  the  bracket  in  Equation  (3.97)  can  be  expressed  as 
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From  Equations  (3.14),  (3.16),  and  (3.20) 

A r  =  d,FM  sin  9  sin  a . 

2 

Also,  making  use  of  Equations  (3.3)  and  (3.24),  Equation  (3.97)  becomes 

where  the  interferometer  array  factor  is  now  given  by 

^  d  ^  A  \j/ 

F IFM(9, a, by/)  =  cos  (iV  +  l);r  —  sin  9  sin  a - —  . 

L  \K)  2  _ 


(3.100) 


(3.101) 


It  is  seen  that  the  array  factor  is  independent  of  y/x .  Only  the  phase  difference  between  the  phase  shifts  applied  to 
the  two  interferometer  elements  plays  a  role  in  steering  the  array  factor. 


The  maximum  value  of  | Fifm  ( 9 ,  OC,  A  ^)j  is  unity  and  occurs  when 

f  <0  -  a  ■  Ay/ 

[N  +  \)n  —  sin  6  sin  a - —  =  mi r 

.  2 


where 


m  =  0,±1,±2,- 


The  minimum  value  of  Km  (9,  a.  Ay/]  is  zero  and  occurs  when 


(iV  +  l);r  —  sm6sina-^-  =  n— 

UJ  2  2 


(3.102) 


(3.103) 


n  —  ±1,±3,±5,*". 


(3.104) 


The  maxima  and  zeros  of  \Fifm(9,(X,Ay\  are  seen  to  be  a  function  of  the  phase  difference,  Ay/ .  Hence,  the 
interferometer  array  factor  can  be  steered  through  an  appropriate  choice  for  Ay/ . 


It  is  convenient  to  measure  A y/  'm  units  of  Jt  radians.  Consequently,  let 

A  y/  =  kn 


(3.105) 


where  k  is  an  appropriately  chosen  constant.  To  consider  the  azimuthal  radiation  pattern,  let  9  equal  the  fixed 
angle,  90  .  The  maxima  are  now  governed  by 
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sm  a  = 


(N  +  Oy-sin  60 

Ao 


(3.106) 


where  m  is  an  integer  such  that 


-+|<(W  +  l{^-l|Sin^|. 

1  \Ao) 


Similarly,  the  locations  of  the  nulls  are  determined  from 

{(«  +  &) 

/v  —  _ _ -  -A-  .1 _ _ _ _  L 


sma  = 


(Ar  +  1)(t)sin^ 


where  n  is  an  odd  integer  such  that 


i(w  +  £)  <{N  +  l^j-  |sin0o|. 


(3.107) 


(3.108) 


(3.109) 


As  was  done  in  Section  3.1,  it  can  be  shown  that  minor  lobes  cannot  exist  except  at  OC  =  ±7T  /  2  where  either  a 
major  lobe,  minor  lobe,  or  inflection  point  is  possible. 


To  gain  insight  into  the  consequences  of  steering  the  interferometer  array,  let 


d=^,e 0=~. 
2  2 


The  array  factor  is  now  given  by 


^ if M  ~  COS 


■1  .  k\ 

— sm  a  —  \n  . 

2J  J 


Equations  (3.106)  and  (3.107),  governing  the  maxima,  simplify  to 


sin  a  = 


2m +  k 
N  + 1 


where  m  is  an  integer  such  that 


k  ^  N  + 1 

m+—  < - 

2  2 


while  Equations  (3.108)  and  (3.109),  governing  the  nulls,  reduce  to 


smor  = 


n  +  k 
N  + 1 


(3.110) 


(3.111) 


(3.112) 


(3.113) 


(3.114) 


where  n  is  an  odd  integer  such  that 
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|h  +  £|<N  +  1.  (3.115) 

For  the  unsteered  interferometer,  it  is  seen  from  Figure  3.11  that  a  major  lobe  occurs  at  a  =  0°  for  all 
integer  values  of  N.  Consider  steering  this  lobe  to  a  =  10°  when  N  =  2  and  N  =  7.  For  the  unsteered  array,  the 
major  lobe  at  a  =  0°  results  when  m  =  0.  Hence,  let  m  =  0  and  Of  =  10°  in  Equation  (3.1 12).  Solution  for  k  yields 


when  N  =  2  and 


k  =  0.522 


k  =  1.392 


(3.116) 

(3.117) 


when  N  =  7.  Polar  plots  of  are  presented  in  Figure  3.16  for  these  two 


cases. 


Comparisons  of  Figure  3.1 1(a)  with  Figure  3.16  (a)  and  Figure  3.11  (f)  with  Figure  3.16  (b)  reveal  that  the 

major  lobe  at  GT  =  0°  has,  in  fact,  been  steered  to  <2  =  10°.  Note  that  symmetry  about  the  y-axis  no  longer  exists 
in  the  steered  arrays  although  symmetry  about  the  x-axis  is  maintained.  For  the  case  in  which  N  =  2,  the  broad  null 
at  OC  =  90°  in  Figure  3.11  (a)  has  disappeared  and  the  two  major  lobes  at  CC  =  41.8°  and  OC  =  138.2°  have 

merged  together  into  a  single  broad  lobe  centered  at  OC  =  90°  in  Figure  3.16  (a).  Also,  at  Of  =  —90° ,  note  that 
there  is  now  a  minor  lobe  in  the  steered  array  whereas  there  is  a  null  in  Figure  3.1 1(a).  For  the  case  in  which  N  =  7, 
the  major  lobe  at  Of  =  90°  in  Figure  3.11  (f)  has  disappeared  and  the  two  major  lobes  at 

Of  =  48.6°  and  OC  —  131.4°  have  merged  together  into  a  single  broad  lobe  centered  at  Of  =  90°  in  Figure  3.16  (b). 
Also,  at  Of  =  -90° ,  the  major  lobe  in  Figure  3.1 1  (f)  has,  in  effect,  split  into  two  major  lobes  at  or  =-55.7°  and 


«  =  -124.3°  as  well  as  introduced  a  minor  lobe  of  opposite  polarity.  We  see  that  severe  distortion  can  occur  in 
the  array  factor  radiation  pattern  when  the  interferometer  is  steered. 


0  r  180 


(a)  N-2  k *0.52 2 

Figure  3.16:  The  azimuthal  radiation  pattern  of  the  interferometer  array  factor  is  plotted  when  the  major  lobe  at  CC  =  0°  in 
Figure  3. 11  is  steered  to  OC  =  10°  for  d  =  2,6 0  =  7T  /  2  ,  and  dIFM  equal  to  (a)  3  /  2  and  (b)  4/L0 


3.3.2  Positioning  the  Broad  Null  by  Steering  the  Interferometer 

A  major  objective  for  steering  the  interferometer  is  to  center  its  broad  null  on  the  main  beam  of  the  radar. 
As  pointed  out  in  Section  3.2,  the  main  beam  of  the  radar  is  located  at  a  =  0° .  This  is  illustrated  in  Figure  3.15. 
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Consequently,  using  the  results  of  Section  3.3.1,  we  attempt  to  locate  a  broad  null  at  (X  =  0°  by  steering  the 
interferometer. 

We  first  consider  the  case  for  N  =  2,  d  =Ao  / 2  ,  and  90  =7t  12 .  As  seen  in  Figure  3.1 1  (a),  a  broad  null 

occurs  at  (X  —  —90°  .  This  resulted  for  ft  =  — 3 ,  in  Equation  (3.50).  Hence,  let  «  =  — 3 ,  (X  =  0°  and  N  =  2  in 
Equation  (3.1 14)  in  order  to  rotate  the  broad  null  at  (X  =  —90°  to  (X  —  0° .  Solution  for  k  yields 

k  =  3.  (3.118) 

,  as  obtained  using  Equation  (3.111),  is  shown  in  Figure  3.17  (a).  As  expected, 

there  is  a  null  at  (X  =  0°  .  However,  it  is  not  a  broad  null !  This  is  confirmed  by  noting  that  the  polarities  of  the 

major  lobes  at  (X  =  ±19.5°  are  of  opposite  signs.  Interestingly,  the  major  lobes  in  Figure  3.17  (a)  are  positioned  at 
the  null  in  Figure  3.1 1  (a)  and  the  nulls  in  Figure  3.17  (a)  are  positioned  at  the  major  lobes  in  figure  3.1 1  (a). 


A  polar  plot  of 


7t 


F 'ifm  ~  ,  a,3n 


<a)N-2  k-3  (b)  N*3  k-1 


Figure  3.17:  The  azimuthal  radiation  pattern  of  the  interferometer  array  factor  is  plotted  when  nulls  in  Figure  3.1 1  are  steered  to 
CX  —  0°  for  d  =  /  2 , 0o  =  K  /  2  and  dIHM  equal  to  (a)  3Xo  /  2  and  (b)  2Ag  . 


Having  failed  with  N  even,  an  attempt  is  now  made  to  locate  a  broad  null  at  (X  =  0°  for  N  odd.  Consider 
N  =  3,  d  =  \  12 ,  and  90  =X  I 2 .  As  seen  in  Figure  3.11  (b),  a  narrow  null  occurs  at  CX  =  —14.5°  when 

ft  =  —1  in  Equation  (3.50).  Hence,  let  ft  =  —1,  CX  =  0°  and  N  =  3  in  Equation  (3.1 14)  in  order  to  rotate  the  null  at 
CX  =  —1 4.5°  to  CX  =  0°  .  Solution  for  k  yields 

k  =  1.  (3.119) 

A  polar  plot  of\F!FM  {x  /  2,  CX,X^ ,  as  obtained  using  Equation  (3.11 1),  is  shown  in  Figure  3.17  (b).  As  expected, 

there  is  a  null  at  CX  =  0° .  However,  because  of  the  alternation  in  polarity  of  the  major  lobes  located  at 
CX  =  ±14.5°,  the  null  at  Cf  =  0°  is  not  a  broad  null!  Interestingly  enough,  there  are  now  broad  nulls  at 

CX  =  ±90  As  was  the  situation  with  N  =  2,  the  major  lobes  in  Figure  3.17  (b)  are  positioned  at  the  nulls  in  Figure 
3.1 1  (b)  and  the  nulls  in  Figure  3.17  (b)  are  positioned  at  the  major  lobes  in  Figure  3.1 1  (b). 

Although  it  was  possible  to  steer  nulls  to  CX  =  0°  for  N=2  and  N=3,  these  were  not  broad  nulls.  This  raises 
the  question,  “Is  it  possible  to  place  a  broad  null  at  CX  =  0°  ?”  From  Equation  (3.114)  the  nulls  of 
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occur  when 


IFM 


(n  \ 

—  ,a,k7C 

U  J 


sin  a  = 


n  +  k 
N  + 1 


where  n  is  an  odd  integer.  The  solution  for  k  with  (X  =  0°  is 


(3.114) 


Thus,  k  is  an  odd  integer.  Note  that 


COS 


x  +  - 


k  =  —n  . 

kK 


=  ±sinx 


for  k  an  odd  integer.  It  follows  that  the  array  factor  in  Equation  (3.111)  becomes 


IFM 


7C 


,a,kn  ~  ±sin 


(N  +  \)tc  . 

- - —sin  a 


(3.120) 


(3.121) 


(3.122) 


when  a  null  is  steered  to  a  —  0°  .  To  have  a  broad  null  at  (X  —  0° ,  it  is  necessary  that  the  major  lobes  located  on 
either  side  of  (X  =  0°  have  the  same  polarity.  However,  this  is  not  possible  because 


IFM 


( Jl  (  ji 

— ,  i x ,  ku  =  —Fjp M  —  —cx,  kn 

V2  )  \2  j 


(3.123) 


It  is  concluded  that  it  is  not  possible  to  steer  a  broad  null  to  (X  =  0° .  In  addition,  the  fact  that  the  unsteered 
interferometer  has  an  array  factor  given  by 


Fifm  ^y,«,oj- 


cos 


fN  +  1 


V 


brsm  a 


J 


(3.124) 


while  the  interferometer  steered  to  have  a  null  at  (X  =  0°has  an  array  factor  given  by  equation  (3.122)  explains  why 
the  maxima  of  one  occur  at  the  nulls  of  the  other. 


3.3.3  Placement  of  Interferometer  Along  the  y-Axis 

To  obtain  a  broad  null  at  (X-  0°,  it  is  necessary  to  place  the  interferometer  on  the  y-axis  while  the  linear 
array  of  the  main  radar  remains  along  the  x-axis.  As  before,  let  the  two  dipole  antennas  of  the  interferometer  pair  be 
directed  in  the  z-direction.  However,  as  shown  in  Figure  3.18,  the  elements  are  now  placed  on  the  y-axis  at 
y  —~dlFM  /  2  and  y  =  +  dIFM  /  2 ,  respectively.  The  corresponding  geometry  for  evaluation  of  the  electric  field 
at  point,  P,  is  illustrated  in  Figure  3.19. 
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The  analysis  proceeds  exactly  in  the  same  manner  as  in  Section  3.1  where  the  interferometer  was  placed 
along  the  x-axis.  The  only  difference  is  that,  in  place  of  J ,  the  angle,  8 ,  is  involved  where  8  is  the  angle  between 

Cl r  and  Cl  ,  the  unit  vector  in  the  y-direction.  Consequently,  cos  5  is  the  projection  of 


Clr  on  Qy  and  is  given  by 


Cos  8  =  Clr  .  Gy  =  sin  0  sin  (j)  =  sin  0  cos  (X  . 


(3.125) 


►  y 


Figure  3.18:  The  dipole  antennas  are  placed  along  the  y-axis  at  A  and  B  and  are  directed  in  the  z-direction. 


P 


Figure  3.19:  The  geometry  is  shown  for  evaluation  of  the  electric  field  at  point,  P. 


Corresponding  to  Equation  (3.14),  we  have 

A  r  =  djFM_  CQS  g  _  djFM_  sin  e  cos  ft  (3.126) 

2  2 
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(3.127) 


In  the  far  field  it  follows  that  the  total  electric  field  at  point,  P,  is  given  by 

EJZ  - 

r 

where  the  interferometer  array  factor  is  given  by 

d 

F  ,FM  (6,  a )  =  cos[7T  (  )  sin  0  cos  a  ].  (3.128) 

AO 

This  is  identical  to  the  expression  for  the  array  factor  given  in  Equation  (3.22)  except  that  sin  OC  has  been  replaced  by 
cos  Cl! . 


Because  the  interferometer  elements  no  longer  straddle  the  main  radar  elements  on  the  x-axis,  the  spacing 
between  the  interferometer  elements  is  no  longer  given  by  Equation  (3.24).  Instead  it  is  convenient  to  measure  the 
interferometer  spacing  in  units  of  half  wavelength.  Therefore,  the  interferometer  spacing  is  defined  to  be 


dlFM 


(3.129) 


where  ks  is  a  constant.  When  d,  the  spacing  between  the  main  radar  antenna  elements,  is  a  half  wavelength,  ks  is 
analogous  to  (N+l)  of  Equation  (3.24).  Substitution  of  Equation  (3.129)  in  Equation  (3.128)  yields 


k  K 

F  ( #,  a )  =  cos[  -£—  sin  #  cos  a  ]. 


(3.130) 


The  maxima  of  \Fifm{6,  OC) |  occur  for 


2m 

Cos  a  =  - - 

ks  sin# 


(3.131) 


where  m  is  an  integer  such  that 


\m\< 


k .  sin# 


Similarly,  the  nulls  of  \F1fm{6,  a) |  are  governed  by 


(3.132) 


n 

Cos  OC  = - 

k  sin  # 


(3.133) 


where  n  is  an  odd  integer  such  that 

|n|  <  ks  sin  # .  (3.134) 

Polar  plots  of  the  azimuthal  radiation  pattern  of  the  interferometer  array  factor,  where  the  elements  are  on 
the  y-axis,  are  presented  in  Figure  3.20  for  0  =  nil  and  k  s  =  3,  5,  and  7.  Note  that  these  values  of  k  s  correspond  to 
diFM  =  3AC  /  2, 5Aa  /  2  and  7  A0  /  2 ,  respectively.  As  expected,  because  k  s  is  an  odd  integer  for  each  of  the  plots, 

a  broad  null  is  seen  to  exist  at  a  =  0  0  in  each  case.  Comparison  of  Figures  3.11  (a),  (c),  and  (e)  with  Figures  3.20 
(a),  (b),  and  (c),  respectively,  reveals  that  the  horizontal  radiation  pattern  of  the  interferometer  positioned  on  the  y- 
axis  is  identical  to  that  of  the  interferometer  positioned  on  the  x-axis  except  that  one  is  rotated  90  0  with  respect  to 
the  other. 
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Figure  3.20:  The  azimuthal  radiation  pattern  of  the  interferometer  array  factor,  where  the  elements  are  on  the  y-axis,  is  plotted 
for  0  =  K  /  2  and  d  IFm  equal  to  (a)  3Ag  /  2 ,  (b)  5  /  2  and  (c)  1  /  2  . 


3.3.4  Evaluation  of  Broad  Null  at  a  =0° 

In  Section  3.3.2  the  interferometer  was  steered  to  produce  a  null  at  a  =  0°  ,  as  shown  in  Figure  3.17  (a). 
However,  this  was  not  a  broad  null  .  In  Section  3.3.3,  a  broad  null  was  located  at  a  =  0°  by  positioning  the 
interferometer  elements  on  the  y-axis,  as  shown  in  Figure  3.20  (a).  In  both  cases  the  spacing  between  the 
interferometer  elements  is 


dlFM  - 


-  34, 


(3.135) 


It  is  of  interest  to  compare  the  width  of  the  steered  null  in  Figure  3.17  (a)  to  that  of  the  broad  null  in  Figure  3.20  (a). 
From  Equation  (3.122)  the  horizontal  radiation  pattern  plotted  in  Figure  3.17  (a)  is  given  by 

| FiFM&’a-2n:)  =  sin  sin  a  (3.136) 

whereas  the  horizontal  radiation  pattern  plotted  in  Figure  3.20  (a),  as  obtained  from  Equation  (3.130),  is 

(3.137) 


3  K 

^IFM 

COS 

—  cos  a 

L  2  J 

Converting  to  dB,  Equations  (3.136)  and  (3.137)  become 


201og 


10 


^IFM  (  2  ’  a’37U) 


=  20  log 


10 


sin 


In  . 

— sin  a 
2 


(3.138) 


and 


20  log 


10 


=  20  log 


cos 


3  n 


-cos  a 


(3.139) 


Evaluations  of  Equations  (3.138)  and  (3.139)  about  the  nulls  at  CX  =0°  are  tabulated  in  Table  3.1.  With  reference  to 
the  table,  note  that  the  width  of  the  steered  null,  when  measured  between  the  -30  dB  points,  is  approximately  0 . 8° 
while  that  of  the  broad  null  is  in  excess  of  1 2°.  Hence,  at  the  -30  dB  points,  the  broad  null  is  more  than  1 5  times 
wider  than  the  steered  null.  At  the  -20  dB  points,  the  width  of  the  steered  null  is  approximately  2.5°  while  that  of 
the  broad  null  is  in  excess  of  22°.  Thus,  at  the  -20  dB  points,  the  broad  null  is  more  than  8.8  times  wider  than  the 
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steered  null.  It  is  concluded  that  the  broad  null  is  considerably  more  effective  than  the  steered  null  in  establishing 
spatial  orthogonality  between  the  interferometer  and  radar  signals. 

For  the  purpose  of  a  later  discussion,  the  broad  nulls  of  Figure  3.20  (b)  and  (c)  are  also  examined.  The 
horizontal  radiation  pattern  plotted  in  Figure  3.20  (b)  is 


p 

J  IFM 

(it  \ 
-Z’a 

12  J 

~5it 

cos 

— cos  a 

L  2  J 

(3.140) 


while  that  plotted  in  Figure  3.20  (c)  is 


(it  N 

'in 

Fl™\2,0C; 

~ 

cos 

— cos  a 

L  2  J 

(3.141) 


Converting  to  dB,  as  in  Equations  (3.138)  and  (3.139),  the  results  are  tabulated  in  Table  3.2.  For  the  broad  null  in 
Figure  3.20  (b)  the  width  is  approximately  10°  when  measured  between  the  -30dB  points  and  is  approximately  18° 
when  measured  between  the  -20  dB  points.  For  the  broad  null  in  Figure  3.20  (c)  the  width  is  in  excess  of  8°  when 
measured  between  the  -30  dB  points  and  is  in  excess  of  14°  when  measured  between  the  -20  dB  points.  In  general, 
the  width  of  the  broad  null  is  narrowed  as  diFM  increases  by  selecting  larger  values  of  ks. 


Table  3.1 

3X, 

Comparison  of  Steered  Null  to  Broad  Null  at  a  =  0°  for  dIFM  -  — ~ 

Zi 

Steered  Null  (N=2)  Figure  3.17  (a) 

Broad  Null  (ks  =  3)  Figure  3.20  (a) 

a 

(Degrees) 

( 7t 

Fjfm  —’a->3ltj 
(dB) 

a 

(Degrees) 

FlFM^2,a, 

(dB) 

0 

-oo 

0 

-OO 

+0.1 

-41.7 

±1 

-62.9 

±0.2 

-35.7 

±2 

-50.8 

±0.3 

-32.2 

±3 

-43.8 

±0.4 

-29.7 

±4 

-38.8 

±0.5 

-27.7 

±5 

-34.9 

±0.6 

-26.1 

±6 

-31.7 

±0.7 

-24.8 

±7 

-29.0 

±0.8 

-23.6 

±8 

-26.7 

±0.9 

-22.6 

±9 

-24.7 

±1.0 

-21.7 

±10 

-22.9 

±1.1 

-20.9 

±11 

-21.3 

±1.2 

-20.1 

±12 

-19.8 
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Table  3.2 

5/1  7/1 

Evaluation  of  Broad  Null  at  a  =  0°  for  d!FM  =  — —  and  — — 

2  2 

a 

(Degrees) 

Broad  Null  (ks  =  5) 
Figure  3.20  (b) 

(dB) 

Broad  Null  (ks  =  7) 
Figure  3.20  (c) 

Fl™{l,a) 

(dB) 

0 

-OO 

-OO 

±1 

-58.5 

-55.5 

±2 

-46.4 

-43.4 

±3 

-39.4 

-36.4 

+4 

-34.3 

-31.4 

±5 

-30.5 

-27.6 

+6 

-27.3 

-24.4 

±7 

-24.6 

-21.7 

±8 

-22.3 

-19.4 

±9 

-20.3 

-17.4 

±10 

-18.4 

-15.5 

+11 

-16.8 

-14.0 

±12 

-15.4 

-12.5 

The  behavior  of  the  vertical  pattern  at  the  broad  null  is,  also,  of  interest.  Letting  OC  =  0°  in  Equation 
(3.130),  the  expression  for  the  interferometer  array  factor  becomes 


F ifm(  ®>0)  —  cos 


(3.142) 


This  yields  the  vertical  pattern  at  a  =  0°.  These  are  plotted  in  Figure  3.21  for  k5  =  3,  5,  and  7.  In  all  three  cases, 
observe  that  broad  nulls  exist  at  6  =  ±90° .  In  fact,  examination  of  Equation  (3.142)  reveals  that  broad  nulls 

at  6  =  ±90°  always  occur  when  ks  is  an  odd  integer.  Therefore,  assuming  that  the  main  radar  antenna  elements  are 
positioned  on  the  x-axis  while  the  interferometer  elements  are  placed  along  the  y-axis,  broad  nulls  occur  broadside 
to  the  main  radar  antenna  in  both  the  horizontal  and  vertical  planes  whenever  kE  is  an  odd  integer.  Also,  note  that 
the  vertical  patterns  in  Figures  3.21  (a),  (b),  and  (c),  where  the  interferometer  elements  are  along  the  y-axis  and  a  = 
0,  are  identical  to  those  in  Figures  3.14  (a),  (c),  and  (e),  where  the  interferometer  elements  are  along  the  x-axis  and 
a -±k  12. 
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(3.143) 


Converting  the  magnitude  of  Equation  (3.39)  to  dB  results  in 


2Ologlo|F/™(0,O)|}=2Ologlo 


Values  for  Equation  (3.143)  about  6  =  90°  are  tabulated  in  Table  3.3  for  ks  =  3,  5  and  7.  By  comparing  the  entries 
in  Table  3.3  with  the  corresponding  entries  in  Tables  3.1  and  3.2,  it  is  concluded  that  the  broad  nulls  in  the 
horizontal  and  vertical  planes  have  identical  behaviors. 


Figure  3.21 :  The  vertical  radiation  pattern  of  the  interferometer  array  factor,  where  the  elements  are  on  the  y-axis,  is  plotted  for 
a  =  0  and  dIFM  equal  to  (a)  3Aa  /  2  ,  (b)  5 A0  /  2  and  (c)  1  Aa  /  2  . 


Table  3.3 

34  54  74 

Evaluation  of  Broad  Null  at  0  =  90°  for  diFM  =  ,  and 

2  2  2 

a 

Broad  Null  (ks  =  3) 

Broad  Null  (ks  =  3) 

Broad  Null  (ks  =  3) 

Figure  3.21  (a) 

Figure  3.21  (b) 

Figure  3.21  (c) 

8 

|  FIFM 

\Fifu  i0’0] 

(Degrees) 

(dB) 

(dB) 

(dB) 

90 

—OO 

-OO 

-OO 

90+1 

-62.9 

-58.5 

-55.5 

90+2 

-50.8 

-46.4 

-43.4 

90±3 

-43.8 

-39.4 

-36.4 

90±4 

-38.8 

-34.3 

-31.4 

90+5 

-34.9 

-30.5 

-27.6 

90±6 

-31.7 

-27.3 

-24.4 

90±7 

-29.0 

-24.6 

-21.7 

90±8 

-26.7 

-22.3 

-19.4 

90+9 

-24.7 

-20.3 

-17.4 

90±10 

-22.9 

-18.4 

-15.5 

90+11 

-21.3 

-16.8 

-14.0 

90112 

-19.8 

-15.4 

-12.5 

3.3.5  Comparison  of  the  Interferometer  and  Radar  Field  Strengths  Broadside  to  the  Radar 

In  this  section  the  electric  field  at  a  point  in  the  far  field  due  to  the  interferometer  is  compared  to  the 
electric  field  at  the  same  point  due  to  the  radar.  Attention  is  focused  on  those  points  in  the  main  beam  of  the  radar 
which  also  fall  within  the  broad  null  of  the  interferometer.  For  this  comparison,  identical  excitations  are  assumed 
for  the  interferometer  and  radar  array. 

With  reference  to  Equation  (3.21),  the  total  electric  field  due  to  the  interferometer  is 

E  (3.21) 

where  the  array  factor  is  given  by  Equation  (3.130).  Similarly,  from  Equation  (3.56),  the  total  electric  field  due  to 
the  radar  is 


E  (3.56) 

where  the  array  factor  is  given  by  Equation  (3.57).  It  follows  that  the  ratio  of  the  field  magnitudes  is 


7tot 

J1FM 


7tot 

J  RAD 


_  2 

“jvj/we.oOl 


(3.144) 


Conversion  of  the  ratio  to  dB  results  in 


20  log 


10  « 


7tot 

J  IF M 


\pT°T 
^ RAD 


=  2Ologlo{ivJ0,a)|} 


(3.145) 


N 


-20 log  10  f^D(0,a)|}-2Ologlo  - 


Let  the  point  at  which  the  fields  are  measured  be  in  the  x-y  plane  where  the  elevation  angle  is  d  =7t  1 2  . 
Assume  the  main  radar  has  the  azimuthal  radiation  pattern  shown  in  Fig  3.15  for  which  N=6  and  d=Xo/2.  Equation 
(3. 145)  then  becomes 


20  log 


I  pTOT 
^  IF M 


7TOT 
' RAD 


=  20  log 


10 


cos 


kn 


-cos  a 


-20  log 


10 


1  sin  [3tt  sin  or] 


6  sin[fsinar] 


-9.5  (dB). 


(3.146) 


Note  that  the  first  term  on  the  right  side  of  Equation  (3.146),  corresponding  to  the  interferometer  array  factor,  is 
tabulated  in  Tables  3.1  and  3.2  for  ks  =  3,5,  and  7.  These  are  values  of  ks  for  which  a  broad  null  appears  at  a  =  0. 
Values  for  the  second  term,  corresponding  to  the  radar  array  factor,  are  tabulated  in  Table  3.4.  From  Equation 
(3.92),  the  3-dB  bandwidth  of  the  main  beam  was  shown  to  be  approximately  16.9  °.  This  agrees  closely  with  the 
values  given  in  Table  3.4  for  a  =  ±  8  0  and  a=  ±  9°.  Utilizing  Tables  3.1,  3.2,  and  3.4,  results  are  tabulated  in 
Table  3.5  for  interferometer  spacings  corresponding  to  ks  =  3,  5,  and  7. 
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Recall  that  a  variety  of  interferometer  element  spacings  are  possible  because  the  interferometer  elements 
are  placed  along  the  y-axis  while  the  main  radar  elements  are  positioned  along  the  x-axis.  Consequently,  the  main 
beam  of  the  radar  can  be  made  as  narrow  as  desired  to  accommodate  the  width  of  the  broad  null  of  the 
interferometer.  Alternatively,  the  broad  null  of  the  interferometer  can  be  made  as  wide  as  needed  to  accommodate 
the  beamwidth  of  the  radar  main  beam.  The  half-power  points  of  the  radar  coincide  with  the  3-dB  points  of  its  main 
beam.  For  this  example,  the  3-dB  points  are  approximately  equal  to  ±  8.5°.  From  Table  3.5  it  is  seen  that  the 
interferometer  field  magnitude  is  at  least  31.0,  26.6,  and  23.7  dB  less  than  the  radar  field  magnitude,  respectively, 
for  the  cases  corresponding  to  k  s  =  3,  5,  and  7.  Recall  that  these  results  were  obtained  assuming  identical 
excitations  for  the  radar  and  interferometer.  Since  the  interferometer  excitation  is  likely  to  be  considerably  smaller 
than  the  radar  excitation,  placement  of  the  broad  null  of  the  interferometer  at  the  center  of  the  main  beam  of  the 
radar  is  likely  to  be  an  effective  technique  for  preventing  the  interferometer  signal  from  interfering  with  the  desired 
radar  target  returns. 


Table  3.4 

Evaluation  of  Radar  Main  Beam 

a 

(Degrees) 

Main  Beam  (N=6) 
Figure  3.15 

FRAD{l,a) 

(dB) 

0 

0 

±1 

© 

i 

±2 

-0.18 

±3 

-0.34 

±4 

-0.59 

+5 

-1.01 

±6 

-1.39 

±7 

-1.94 

±8 

-2.58 

±9 

-3.2 

±10 

-4.2 

±11 

-5.2 

±12 

-6.4 
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Table  3.5 


Ratio  of  Interferometer  to  Radar  Electric  Field  Magnitudes  in  the  Main  Beam  of  the  Radar  for  d]FM  = 

^,^-and^L 


a 

(Degrees) 

■ 

■ 

§§ 

1 

■ 

■ 

pTOT 

Z'RAD 

(dB) 

pTOT 
^  RAD 

(dB) 

pTOT 

C'RAD 

(dB) 

0 

-  OO 

-OO 

-OO 

±1 

-72.4 

-68.0 

-65.0 

±2 

-60.1 

-55.7 

-52.7 

±3 

-53.0 

-48.6 

-45.6 

±4 

-47.7 

-43.2 

-40.3 

±5 

-43.4 

-39.0 

-36.1 

±6 

-39.8 

-35.4 

-32.5 

±7 

-36.6 

-32.2 

-29.3 

±8 

-33.6 

-29.2 

-26.3 

±9 

-31.0 

-26.6 

-23.7 

±10 

-28.2 

-23.7 

-20.8 

±11 

-25.6 

-21.1 

-18.3 

±12 

-22.9 

-18.5 

-15.6 

3.4  Use  of  Frequency  Modulated  Signals  to  Steer  the  Interferometer 


When  the  interferometer  elements  are  on  the  y-axis,  the  interferometer  array  factor  is 


1FM 


(0,a)  = 


cos 


A 


(d 


IFM 

Xo 


sin  6  cos  a 


Recalling  that 


po  = 


27rf0  _  In 
c  Xo  ’ 


the  array  factor  can,  also,  be  expressed  as 


(3.128) 


(3.3) 
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F  ifm  (#,a)  =  cos  ^-dWM  sin#  cos  a  . 


(3.147) 


The  relation  in  Equation  (3.3)  can  be  interpreted  as  having  arisen  from  the  wave  equation. 


In  the  far  field  the  propagating  wave  can  be  approximated  as  a  spherical  wave  whose  description  depends 
only  upon  the  radial  distance,  r,  and  sinusoidal  frequency  of  excitation,  f0.  Assuming  spherical  symmetry,  the  wave 
equation  is  given  by  [6] 

d2Mr.Q]  =  J_  d 2  Mr,  03 

dr2  c2  dt2 

where  c  is  the  wave  propagation  velocity.  One  solution  to  this  equation  is  of  the  form 


K  -jPor  „  >2^/o' 


,s(r,f)  = — e  ’  °r  e 


(3.149) 


This  can  be  verified  by  direct  substitution  of  Equation  (3.149)  into  Equation  (3.148).  Note  that 

=.  &K  e-”-’ ew-' 
dr2 


(3.150) 


JL  32[ry(r,Q]  =  _  f  2^Y  Kc~jfi.r  e 
c2  3 12  i,  c  J 

Equality  between  Equations  (3.150)  and  (3.151)  exists  provided 

po=^. 

c 


(3.151) 


pO  can,  also,  be  written  as 


c=  tofo 


po  = 


(3.152) 


Examination  of  Equation  (3.144)  suggests  that  the  array  factor  can  be  varied  by  changing  pO  which  is  a 
function  of  the  sinusoidal  frequency  of  excitation.  Therefore,  it  should  be  possible  to  make  the  array  factor  vary  in 
time  by  exciting  the  interferometer  with  a  frequency  modulated  signal.  This  observation  is  examined  in  Section 
3.4.1. 
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3.4.1  Quasi-Stationary  Considerations 


(3.161) 


?{t)= 4 


r#W  +  2  tf(r) 

at 


12 


+  - 


d2P{t)  2xd2g^ 

dt 2  t*2 


Now  assume  that  the  modulation,  g(t),  varies  slowly  relative  to  the  carrier  frequency,  f0,  such  that 

\dg{t)\ 


dt 


«fo- 


Making  use  of  Equation  (3.158),  observe  that 


(3.162) 


and 


dt  c  dt 


t£M=2*r-^ 

dt  c  dt 3 


(3.163) 


(3.164) 


On  the  other  hand, 


It  is  likely  that 


27rf{t)  =  2jz 


fo  + 


dgit) 


dt 


A  r 

tr~~«  fo 


(3.165) 


(3.166) 


where  tr  denotes  the  time  required  by  the  propagating  wave  to  traverse  the  distance,  r.  Finally,  assuming 

\d2g{t)\ 


dt 2 


«fo. 


(3.167) 


it  is  concluded  that  Equation  (3.161)  is  approximately  given  by 

f  2>r(<)  V 


V 


C  ) 


(3.168) 


Therefore,  Equation  (3.158)  is  approximately  valid  and  the  propagating  wave  can  be  approximately  modeled  by 
Equation  (3.157). 


We  say  that  the  interferometer  behaves  in  a  quasi-stationary  manner  when  Equations  (3.149)  and  (3.3), 
which  are  strictly  valid  only  for  a  sinusoidal  constant  frequency  excitation,  can,  to  a  very  good  approximation,  be 
generalized  to  Equations  (3.157)  and  (3.158)  for  frequency  modulated  excitations.  However,  as  pointed  out  by 
Equation  (3.162),  the  modulation  must  vary  slowly  relative  to  the  carrier  frequency.  It  follows  that  the  array  factor, 
given  by  Equation  (3.147),  can  be  generalized  to 
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(3.169) 


^IFM  iPi^X,  t)  ~  COS 


Pit) 


d ip m  sin  6  cos  a 


for  frequency  modulated  excitations  in  this  class.  Introducing  the  time  variant  wavelength,  defined  by 

(3.170) 


!(,)=  c 

the  array  factor  can  also  be  expressed  as 

FlFM(0,a,t)=  cosl  K\ 


M 

( djFM  \  ■ 


sin  0  cos  Of 


L  \^{t)j 

The  time-variant  array  factor  is  investigated  in  the  following  sections  for  selected  modulations. 


(3.171) 


3.4.2  The  Time-Variant  Array  Factor  for  a  Chirp  Signal  Applied  to  the  Interferometer  Elements 

In  this  section  the  same  chirp  signal  is  applied  to  the  two  interferometer  elements.  In  Section  3.4.3  an  up 
chirp  is  applied  to  one  element  while  a  down  chirp  is  applied  to  the  other.  Both  cases  are  handled  by  the  general 
analysis  given  below. 


With  reference  to  Figure  3.18,  let  the  signals  applied  to  the  interferometer  Elements  at  A  and  B  be  given  by 


sA(t)  = 


Vm  cos  (<k?2  (o)] 


T  T 

- <  t  <  — 

2  2 


0  ,  otherwise 


(3.172) 


and 


sB(t)  = 


FmCOS  [M/o'-s  i(0+g2  ('))] 


0 


-L<,<L 

2  2 


otherwise 


(3.173) 


where  gj  (t)  and  g2  (t)  are  suitably  chosen  phase  modulations  such  that  the  quasi-stationary  assumptions  of  Section 
3.4.1  are  valid.  In  terms  of  complex  signals  sA  (t)  and  sB  (t)  can  be  expressed  as 


^(0  = 


Re  {V  ej2"f°‘  ei2ng'{,)  ei2ng2{l'> ) 


L<t<L 

2  2 


0 


otherwise 


(3.174) 


and 
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sB(t)  = 


_T<  kT 

Re  (Vm  ej2nLt  ejlng'w  ej2ng2(t) }  ,  2~f  ~  2 


,  otherwise 


By  definition,  the  instantaneous  frequency  of  sA  (t)  and  sB(  t),  respectively,  is 


(3.175) 


f(f)=L+M>+Ml 

JaX>  °  dt  dt 


(3.176) 


(3.177) 


Using  the  notation  of  Figure  3.19,  and  assuming  quasi-stationary  behavior,  the  total  electric  field 
at  point,  P,  is  given  by  Re  {is  ™TM  }  where 


c  c 


(3.178) 


gTOT  _  K  pig  \e-jpB(‘h ej2tfo‘ e-j2’%i(<)ej2”:g2(t)  p(g  \e-jPA‘)r2 eJ2*fo‘ eJ2,tfo< ej2m{‘)  p  J7^ 
^  IF  M  r  _  2 


With  P  in  the  far  field  the  analysis  proceeds  in  the  same  manner  as  in  Section  3.1 .  Then 


c  \_  dt  dt 


(3.180) 


2tr[,  dgx(t)  dg2(t) 


fo~ 

c 


(r+Ar) 

dt  dt 


(3.181) 


where  A r  is  given  by  Equation  (3.126).  In  the  far  field,  therefore,  Equation  (3.179)  becomes 


F{6)\  e 


T,,  ,dgi«)  ,  ,r  <*«.(')  r1  rfg'(0  .1 

•|.(/*+-S_)Ar  ~  rJ  e-j2*glO)  +eJ  c  La  <*  ’  <*  J  eJ2*gA‘) 


- j-W ^ 

r,  C  L  ■* 


dt  dt  j  eJ'2trg2(t)  gjlnfot 


(3.182) 


JK  -M(foJsA0)r-MlAr) 

=—F(0)FIFM(0,a,t)  e  ^  J  e32KgAt)  em°' 

r 
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where  the  interferometer  array  factor  is  given  by 


FjFM{0,a,t)=  cos{^-  +2 *£,(*)}• 

Finally,  from  Equations  (3.176)  and  (3.177),  note  that 

f  ,  dg2{t)  _  fA{t)+  fB{t) 


Jo  dt 


(3.183) 


(3.184) 


dgx{t) 

dt  2 

As  a  result,  Equations  (3.182)  and  (3.183)  can,  also,  be  expressed  as 

eZ  «—  F(e)Fm(e,a,t) 


(3.185) 


(3.186) 


g-J^K/A ^)+fB{t))r-(fA(t)-fB {*))£*] ^2ng2(t)  jin fQ t 


^ifm  A{t)+ f B(t))Ar  {fA{t )  fs  (O)^]"*"  In  gj  (O}-  (3.187) 

Having  completed  the  above  analysis,  we  now  consider  the  case  for  which  the  same  chirp  signal  is  applied 
to  the  two  interferometer  elements.  Let 


,W=o 


(3.188) 


Then,  with  reference  to  Equations  (3.1 72)  and  (3.1 73),  it  follows  that 

_!<  <L 

sA(t)  =sB(t)  =  -  V’-c0s^27t(<fot+d  V)]»  2  2 

0  ,  otherwise 

The  corresponding  instantaneous  frequency  is 

f  A  (*)  =  fB  fr)  =  fo+ jUt  =  f(t) . 


Observe  that 


/,(')+/»W  =  2(/.+/«) 


(3.189) 


(3.190) 


(3.191) 


(3.192) 
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(3.193) 


Consequently,  when  the  same  chirp  signal  is  applied  to  the  two  interferometer  elements,  the  expression  in  Equation 
(3 . 1 87)  for  the  interferometer  array  factor  becomes 


Fifm (d, a, t)  =  cos  —  (fo  +  JLu)Ar  =cos—f(t)Ar  . 

_  c  J  L  c 

Making  use  of  Equations  (3.1 58)  and  (3.126),  the  array  factor  can,  also,  be  expressed  as 

FWM  (0>  «»*)  =  cos  — ®  d1FM  sin  6  cos  a  . 

Observe  that  this  is  consistent  with  Equation  (3.169). 


(3.194) 


(3.195) 


As  was  done  in  Equation  (3.129),  it  is  convenient  to  measure  the  spacing  of  the  interferometer  elements  in 
units  of  half  wavelength.  Then 

4U™ = 4?  -  -  ^[f. +4TW1 


[l+AOi] 


(3.196) 


where  A^(/)  is  the  instantaneous  frequency  deviation  given  by 

4 =  =&■ 

In  addition,  it  is  convenient  to  define 

Jo 

which  can  be  interpreted  as  the  instantaneous  deviation  from  ks .  Note  that 

Af(l)_Ak,{,) 

f.  K 

Hence,  Equation  (3.196)  can,  also,  be  written  as 

2  im  2  [  ks 

It  follows  that  alternative  expressions  for  the  interferometer  array  factor  are 


Fim{e,a,t)  =  cos--^  Sind 

^  Jo 


cos  a 


(3.197) 


(3.198) 


(3.199) 


(3.200) 


(3.201) 
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Now  assume  a  1%  change  in  the  carrier  frequency  at  time,  t2,  such  that 

A f{t2)  Ak(t2) 

J  v  11  ~ — *-^=0.01. 


The  corresponding  array  factor  is 


„  (it  f  2\2\k 

J  =  C0S — - — cosor  . 


(3.205) 


(3.206) 


Its  magnitude  is  plotted  in  Figure  3.22(b).  Instead  of  the  broad  null  at  (X  =  0° ,  there  is  now  a  minor  lobe  peak 
whose  magnitude  is  0.324.  This  is  only  9.8  dB  down  from  the  major  lobe  peak  value  of  unity  and  raises  concern 
that  significant  energy  could  be  radiated  in  the  direction  of  the  radar  main  beam.  Even  so,  as  is  seen  from  Table  3.6, 
there  is  very  little  steering  of  the  major  lobe  positions. 
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Using  Equation  (3.210),  the  magnitude  of  the  minor  lobe  peak  at  (X  =  0°  is  approximately  given  by 

Ak(t)-  =0.21-  =  0.33.  (3.211) 

W2  2 
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This  value  is  very  close  to  the  peak  value  of  0.324  observed  in  Figure  3.22  (b). 

At  some  time  instant,  t3,  suppose  it  is  desired  that  the  minor  lobe  peak  at  O'  =  0°be  30dB  down  from  the 
major  lobe  peak  of  unity.  From  Equation  (3.210)  it  would  be  required  that 

201ogIO|MJ(*3)f  |  =  -30dB .  (3.212) 

Solution  of  this  equation  results  in 

A*,(/3)  =  0.02. 

As  before,  let  ks  =  2 1 .  Then  Equation  (3 . 1 99)  gives 


A/fc)  M.fe)  0.02  0M1 

f,  K  21 


(3.213) 


This  corresponds  to  a  0.1%  change  of  the  carrier  frequency.  The  corresponding  array  factor  is 


1FM 


(n  \ 

—,a,t3  =cos| 

V2  J 


21.02br 


■cos  a 


(3.213) 


Its  magnitude  is  plotted  in  Figure  3.22  (  c  ).  At  CX  =  0° ,  the  value  of  the  interferometer  pattern  is 


IFM 


=  0.03298 


(3.214) 


which  is  29.6  dB  down  from  the  major  lobe  peak  value  of  unity.  Consequently,  a  0.1%  change,  or  less,  of  the  carrier 

frequency  does  a  very  good  job  of  approximating  the  broad  null  at  (X  =  0° .  However,  as  is  shown  in  Table  3.7, 
there  is  negligible  steering  of  the  major  lobe  positions. 


If  the  broad  null  at  CC  —  0°  is  to  be  maintained,  it  is  concluded  that  a  linear  frequency  modulated  pulse,  or 
chirp  signal,  is  not  effective  for  steering  the  interferometer.  In  the  next  section  we  investigate  the  interferometer 
response  when  an  up  chirp  is  applied  to  one  of  its  elements  while  a  down  chirp  is  applied  to  the  other  element. 


Table  3.7 

Comparison  of  Major  Lobe  Positions  in  First  Quadrant  for  a  0.1%  Change  in  Carrier  Frequency 

4 m 

fo 

Major  Lobe  Positions 

a 

(Degrees) 

0 

17.8 

31 

40.4 

48.2 

55.2 

61.6 

67.6 

73.4 

79 

84.5 

90 

0.001 

17.9 

31.1 

40.4 

48.2 

55.2 

61.6 

67.6 

73.4 

79 

84.5 

90 
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3.4.3  The  Time-Variant  Array  Factor  for  Two  Different  Chirp  Signals  Applied  to  the  Interferometer 
Elements 

With  reference  to  Equations  (3.172)  and  (3.173),  let 

gx{t)  =  llt-  (3.215) 

and 

g2(/)  =  0.  (3.216) 


Therefore,  the  signals  applied  to  the  interferometer  elements  at  A  and  B  are 

_t_<  <t_ 

Vm  cos[2 z(f0t+JiT)\>  2  ~  *  ~  2 

0  ,  otherwise 

_  T_  <  <T_ 

sB(t)  =  .  ^cos[2>r(/./-/i^)]f  2  ~  1  ~  2 

0  ,  otherwise 


(3.217) 


(3.218) 


The  corresponding  instantaneous  frequencies  are 

fA{t)=fo+Vt 


Note  that 


/.(')+/,(<)  =  V. 


(3.219) 

(3.220) 

(3.221) 


while 

fA{t)-fB(t)  =  20. 

Thus,  with  an  up  chirp  applied  to  one  of  the  interferometer  elements  and  a  down  chirp  applied  to  the  other,  the 
interferometer  array  factor  of  Equation  (3.187)  becomes 


F1FM  (#>  a.  U  r)  =  cos|~  I fobr-tM*  2Jtfi  y| . 


(3.222) 


Observe  that  the  array  factor  now  depends  upon  r  as  well  as  9,  (X,  and  t.  The  dependence  upon  0  and  a  is  made 
obvious  by  utilizing  Equation  (3.126)  which  yields 
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2jt  ~  f  d  t2 1 

Fifm  (^.  a>l> r)  ~  cosl —  “ — ^-sin^cos  a-fltr  +2  ji /x — >.  (3.223) 

[c[_  2  J  2  j 

It  is  of  interest  to  evaluate  the  array  factor  at  a  fixed  point  in  space  for  which  the  coordinates  are  60,CC0 , 
and  r0 .  Making  use  of  Equation  (3.129),  note  that 

2a-  fodjpu  qq  cosa°  _  /A  sjn^  cos a  =-^sin 6  cosor  =  (j)^  (3.224) 

c  2  c  2  2 

where  (j)AF  is  a  constant.  From  Equations  (3.219)  and  (3.220)  the  maximum  frequency  deviation  of  the  up  and  down 
chirp  signals  is 

,  ,(T\  ,(  T\  J  T\  JT\  ,  nT  ,  pT  ^ 

f°=f{irf\-irf{-i)-f{2)=f° *2  -f° +i =tiT  <3-225) 


Hence, 


It  follows  that 


2  k  2  n  fD 

Mro  ~  T  ro  —  iF 

c  c  T 


(3.226) 


(3.227) 


Consequently,  at  the  point  in  space  for  which  the  coordinates  are  Ba ,  OLo ,  and  r  ,  the  interferometer  array  factor 
can  be  written  as 


r  t 2 

fifm (&o ,a0,t,r0)  =  cos  2 jtfAFt -Inn  —  - <Paf 


(3.228) 


where  fAF  and  <p ^  are  array  factor  constants  dependent  upon  position.  Thus,  at  any  point  in  the  far  field,  the  array 
factor  is  a  linear  frequency  modulated  pulse  whose  carrier  frequency  and  phase  are  determined  by  position. 

From  Equation  (3. 1 82)  the  total  complex  electric  field  at  a  point  in  the  far  field  is  given  by 


(3.229) 


Therefore,  the  total  real  electric  field  is 


Re{E™  }~—F(e)Fm(0,a,t,r) 


(3.230) 


cos  2 n(L+t*\-IS 


As  before,  consider  a  point  whose  coordinates  are  Q0,(X0,  and  r0.  Following  the  previous  development,  note  that 


60 


flA r 
c 


— —  ^ 1FM  sin  6.  cos  an 
c  T  2 


_  1  Id  ksK  sing>  cos  Or.  (3.231) 

Vo  T  4  00 

=  cosa0  =  fCA 

J  o 

is  a  carrier  offset  frequency.  In  addition,  define  the  carrier  phase  to  be 

(/)ca=2k^-  .  (3.232) 

c 

The  total  real  electric  field  at  the  point  whose  coordinates  are  60,Ct0,  and  ro  is  then  given  by 

Re{«™  }«  —  F(».  )F*«  («.  r.) 

r 

o 

COS  [24/;+/Jf-(U  (3.233) 

where  fCA  and  (f>CA  are  carrier  constants  dependent  upon  position. 

The  substitution  of  Equation  (3.228)  into  Equation  (3.233)  results  in 

Rete}“— ^ F(0„)cos 

ro 

cos  [2/r(/o+/c<>-0j.  (3.234) 

This  is  recognized  as  a  double  side  band  suppressed  carrier  amplitude  modulated  signal  where  the  amplitude 
modulation  is  a  linear  frequency  modulated  waveform  whose  instantaneous  frequency  decreases  during  the  chirp.  It 
is  interesting  to  note  that  the  carrier  is  a  constant  frequency  sinusoid  even  though  an  up  chirp  and  down  chirp  were 
applied  at  the  interferometer  elements.  At  t  =  0  the  electric  field  is  given  by 

Re{ej ™  }  |,=0  I  -  —  F{eo  )cos(-  <t>AF  )cos(-  (f>CA  ) .  (3.235) 

ro 

This  can  be  a  small  value  at  various  points  in  the  far  field.  A  typical  sketch  of  the  electric  field  is  shown  in  Figure 
3.23,  assuming  fAF  to  be  much  smaller  than  (f0  +  fCA). 

Although  fAF ,  (f)AF ,  fCA  and  <pCA  change  with  position,  the  basic  waveform  for  the  electric  field  is 

unchanged.  Consequently,  steering  the  interferometer  with  up  and  down  chirps  applied  to  the  interferometer 
elements  has  little  meaning  even  though  the  array  factor  is  a  function  of  both  time  and  the  observation  point  in  the 
far  field.  In  addition,  the  waveform  sketched  in  Figure  3.23,  is  not  a  good  choice  for  masking  the  direct  path  signal 
from  the  main  radar. 


27tfAFt-2KH  —  -<t>AF 
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It  is  shown  in  Section  3.4.4  that  a  frequency-hopped  waveform  is  capable  of  maintaining  the  broad  null 
while  steering  the  interferometer  and  degrading  performance  of  the  non-cooperative  bistatic  radar. 


0.5  1 

t (secs) 


Figure  3.23:  When  up  and  down  chirps  are  applied  to  the  two  interferometer  elements,  the  electric  field  at  a  fixed  point  in  the  far 
field  is  a  suppressed  carrier  amplitude  modulated  signal  whose  modulation  is  a  linear  frequency  modulated  waveform,  as 

sketched  in  the  figure. 


3.4.4  The  Time-Variant  Array  Factor  for  a  Frequency-Hopped  Signal  Applied  to  the  Interferometer 
Elements 


It  is  intended  that  the  interferometer  signal  be  used  to  1)  mask  over  many  directions  the  direct  path  signal  of 
the  main  radar  and  2)  communicate  separate  information  to  various  receivers  located  at  different  angles.  In  addition, 
the  interferometer  signal  should  not  interfere  with  the  mission  of  the  main  radar.  In  order  to  maintain  spatial 

orthogonality  between  the  interferometer  and  main  radar,  it  is  desirable  that  the  broad  null  at  a  =  0"  be  preserved 
while  the  interferometer  is  steered  in  different  directions.  These  objectives  can  be  met,  as  shown  in  this  section,  by 
using  a  frequency-hopped  signal  to  drive  the  two  elements  of  the  interferometer. 


Let  both  elements  of  the  interferometer,  assumed  to  be  positioned  on  the  y-axis,  be  driven  by 


sifm (0  —  0  iT.) 

1=0 


where 


t>i(t)  = 


Vm  cos  \lK  (f0  +  foi )  t  ] ,  0  <  t  <  Ts 


0  , 


otherwise 


Consequently,  during  the  i  —  interval,  given  by 

iTt<t<(i  +  \)Ts, 

the  interferometer  signal  is  a  constant  frequency  sinusoid  with  frequency,  {fa  +  foj). 


(3.236) 


(3.237) 


(3.238) 


Following  the  quasi-stationary  development  in  Section  3.4.2,  the  interferometer  array  factor  can  be 
expressed  as 


Fmt  (0,  a,  t)  =  cos  dIFM  sin  6  cos  a 


(3.195) 
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where 


2 nf(t) 


During  the  i-  interval 


/}(>) 


=  2f(4 ±2"),  ,T,  <  ( <  (i  + 1)7; 


As  in  Section  3.3.3,  let  the  spacing  between  the  interferometer  elements  be 


d  —  k  — 
uifm  -  Ks  2 


Then 


P{t)d1FM  _  ks7[X0  r 

2  2c  1/0  ‘/o,J  2 


|  |  _j_  J  oi 


fol 

fo 


where  use  was  made  of  the  fact  that 


Now  define 


Note  that 


C  =  Kfo  ■ 


Ak,  =  . 

‘  fo 


Ll  =  ^l 

fo  K 


Hence,  Equation  (3.240)  can  be  written  as 

fi{t)d if M  _  k/C 


1  + 


Ak; 


=  {ks+Akl)^. 


Substitution  of  Equation  (3.243)  into  Equation  (3.195)  yields 


P ~'lFM  (^»  ~ 


(ks  +  Ak')—  sin#  cos  # 


When  6  -It  /  2 ,  the  horizontal  radiation  pattern  becomes 


1FM 


( n  \ 
-,a,t 
V2 


=  cos 


(K  +  Akt  )^-  cos  a 


(3.158) 

(3.239) 

(3.129) 

(3.240) 

(3.152) 

(3.241) 

(3.242) 

(3.243) 

(3.244) 

(3.245) 
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A  broad  null  exists  at  OC  =  0°  provided  (ks  +  Akt )  is  an  odd  integer.  To  maintain  a  broad  null  at  f0  (i.e.,  when 
Aki  —  0),  let  ks,  also,  be  an  odd  integer.  The  broad  null  is  then  preserved  by  requiring  that  Akt  be  plus  or  minus 
an  even  integer.  For  specified  values  of  A ki ,  k$  and  fa ,  fgi  is  determined  from  Equation  (3.242)  according  to 


f  = 

J  Ol 


(3.246) 


In  this  way,  a  broad  null  always  exists  at  OC  =  0°  even  though  the  frequency  changes  from  one  interval  to  the  next. 

The  interferometer  pattern  of  Equation  (3.245)  can  be  steered  in  various  directions  through  the  appropriate 
choice  of  Akt  and,  therefore,  foi .  In  general,  when  (ks  +  Ak^is  an  odd  integer,  there  are  (ks  +  A kt  —  l)  major 

lobes  completely  contained  in  the  upper  half  of  the  interferometer  pattern  centered  at  OC  =  0° , 
(ks  +A kt  —  l)  major  lobes  completely  contained  in  the  lower  half  of  the  interferometer  pattern  centered  at 

OC  =  ±1 80° ,  and  two  major  lobes  centered  at  (X  =  ±90° .  Thus,  the  number  of  major  lobes  changes  with  a  change 
in  Akt .  This  can  be  utilized  to  direct  the  interferometer  signal  towards  a  specified  direction  during  the  i-  interval. 

To  illustrate  the  steering  of  the  interferometer  by  means  of  a  frequency-hopped  signal,  while  maintaining  the 
broad  null  at  OC  =  0° ,  radiation  patterns  are  presented  in  Figure  3.24  for  ks  =  21  and  A kl  =  0  and  ±  2  .  Note  that 

this  involves  the  smallest  allowable  values  for  A&(.  | .  Of  course,  greater  differences  between  the  patterns  will  result 

for  larger  values  of  M- 


As  expected,  a  broad  null  exists  at  OC  =  0°  for  all  three  cases  presented  in  Figure  3.24.  Observe  that  there  are 
18,  20,  and  22  major  lobes,  respectively,  in  the  upper  half  of  the  patterns  shown  in  Figure  3.24  (a),  (b),  and  (c). 
Also,  each  pattern  contains  two  major  lobes  centered  at  OC  =  ±90° . 


Positions  of  the  major  lobes,  indicated  by  x’s,  and  nulls,  indicated  by  o’s  ,  are  shown  in  Figure  3.25  for  the  three 
patterns  plotted  in  Figure  3.24.  The  polarities  of  the  major  lobes  are  indicated  by  +  and  -  signs  placed  over  the  x’s. 


Comparison  of  the  major  lobe  positions  in  Figure  3.24  (a)  with  those  in  Figure  3.24  (b)  reveals  that  the 
major  lobes  in  the  two  figures  are  separated  by  a  minimum  of  1°  for  20°  <  OC  <  80°  .  A  similar  conclusion  arises 
when  the  major  lobe  positions  in  Figure  3.24  (b)  are  compared  with  those  in  Figure  3.24  (c).  Making  the  same 
comparisons  with  regard  to  the  nulls,  it  is  seen  that  they  are,  also,  separated  by  a  minimum  of  1°  for 
20°  <  OC  <  80°  .  Consequently,  significant  separations  exist  even  though  {ks  +  Aki )  differs  by  the  minimum 

allowable  value. 
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Figure  3.24:  The  azimuthal  radiation  pattern  of  the  interferometer  array  factor  is  plotted  for  a  frequency-hopped  excitation  with 

9  =  7T !2,ks  =  21,  and  A  k(  equal  to  (a) -2,  (b)  0,  and  (c)  +2. 


Figure  3.25:  Positions  of  the  major  lobes,  indicated  by  x’s,  and  nulls,  indicated  by  o’s,  are  shown  for  the  three  patterns  plotted  in 

Figure  3.24. 

However,  even  though  the  difference  in  (ks  +  Akt)  is  twice  as  large  when  comparing  the  pattern  of  Figure 

3.24  (a)  with  that  of  Figure  3.24  (c),  observe  for  54°  <  a  <  66°  that  some  separations  between  major  lobe 
positions  and  null  positions  are  less  than  1°.  Consequently,  making  the  difference  in  (ks  +  Akj )  larger  does  not 

guarantee  large  separations  relative  to  neighboring  major  lobe  and/or  neighboring  null  positions.  It  should  be  noted 
that  some  of  the  closely  spaced  major  lobes,  when  comparing  Figure  3.24  (a)  to  3.24  (c),  are  of  opposite  polarity. 
The  reversal  in  sign  could  help  to  thwart  the  non-cooperative  radar  in  its  attempt  to  obtain  a  high  quality  reference 
from  the  direct  path  signal. 

Nevertheless,  this  example  indicates  that  it  should  be  possible  to  select  suitable  values  for  ks  and  A kt  such 
that  the  interferometer  is  steered  in  desired  directions  while  maintaining  the  broad  null  at  a  =  0° . 
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CHAPTER  4 

A  RADAR  WAVEFORM  EMPLOYING  FREQUENCY  DIVERSITY 


4.1  Frequency-Hopped  Waveforms  [7],  [8] 

Frequency-hopped  waveforms  are  a  subclass  of  frequency-coded  signals.  In  particular,  consider  the  radar 
pulse  waveform  of  duration  NT  which  consists  of  N  contiguous  different  frequency  bursts,  each  of  duration  T. 
Assume  the  hopping  frequencies  are  appropriately  selected  from  a  set  of  N  frequencies  equally  spaced  according  to 
(fc+fn)  where 


f 


n 


n  =  0,1,  ...,(N  -1)  . 


(4.1) 


The  spacing  is  given  by 

A/  =  (/c  +  f„+i )  -  (fc  +  /„ )  =  ~~~  n  =  0,1, (N  -1)  .  (4.2) 

A  linear  stepped-firequency  pulse  results  when  the  frequencies  are  chosen  consecutively  such  that  the  n* 
frequency  burst  has  the  frequency  (fc+f„),  where  fn  is  given  by  Equation  (4.1).  On  the  other  hand,  the  frequencies 
can,  also,  be  scrambled.  For  this  purpose,  consider  an  arbitrary  ordering  of  the  N  integers  contained  between  0  and 
(N-l).  Note  that  there  are  N!  such  orderings.  Let  a  particular  sequence  of  ordered  integers  be  represented  by  the  set 

{ea}  =  {e0,ex,e2,...,e^_  ,>  .  (4.3) 

Now  let  the  subscript,  n,  indicate  the  ordering  of  the  frequency  bursts  contained  within  the  radar  pulse.  The 
frequency  of  the  n'h  burst  is  then  chosen  to  be 

f.  +  f,  =  f.+  y  ;  »=<u  ■-  <N-»  ■  <44> 

The  resulting  waveform  is  referred  to  as  a  scrambled-ffequency  pulse.  Because  {0n}  contains  each  integer  between 
0  and  (N-l)  only  once,  the  frequencies  of  the  N  bursts  within  the  radar  pulse  are  each  different. 

Analytically,  the  frequency-hopped  pulse  can  be  expressed  as 

s(t)=  t  s„(t-nT)  (4.5) 

n=0 


where 


Sn(t) 


A  cos[2;r(fc  +fn)t] , 
0 


0  <  t  <T 
otherwise 


(4.6) 


The  envelope  of  this  waveform  is  illustrated  in  Figure  4.1  for  both  a  linear  stepped-frequency  pulse  and  a  scrambled- 
ffequency  pulse.  When  the  frequencies  are  chosen  consecutively,  the  linear  stepped-frequency  pulse  is  a  discrete 
approximation  to  the  linear  frequency  modulated  pulse  or  chirp  waveform.  The  bandwidth  of  the  frequency-hopped 
pulse  is  NAf,  where  Af  is  given  by  Equation  (4.2),  and  is  assumed  to  be  much  less  than  the  carrier  frequency,  fc 
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Figure  4.1 :  The  frequencies  of  a  frequency-hopped  pulse  may  be  either  (a)  consecutive  or  (b)  scrambled. 


Many  factors  favor  the  use  of  frequency-hopped  waveforms  for  the  pulsed  radar  signal.  These  are  briefly 
discussed  below: 

1)  By  appropriate  ordering  of  the  frequencies,  it  is  possible  to  eliminate  the  range-Doppler  ambiguity  of  the 
linear  frequency  modulated  or  chirp  waveform.  It  is  shown  in  the  following  sections  that  a  thumbtack 
approximation  to  the  ambiguity  function  is  possible. 

2)  For  a  specified  value  of  N  there  are  multiple  orderings  of  the  frequencies  that  result  in  a  thumbtack  shaped 
ambiguity  function.  Coherent  reference  denial  is  enhanced  by  changing  the  frequency  scrambling  from  one 
radar  pulse  to  the  next.  This  prevents  the  non-cooperative  bistatic  radar  from  relying  on  the  integration  of  a 
single  radar  waveform  estimate  over  the  coherent  processing  interval.  In  addition,  the  non-cooperative 
radar  is  forced  to  repeatedly  estimate  a  new  radar  waveform  from  each  pulse  received  in  the  direct  path. 

3)  Another  benefit  from  modifying  the  frequency  ordering  from  one  radar  pulse  to  the  next  is  that  the  main 
radar  is  likely  to  be  less  susceptible  to  interfering  FM  signals  and  repeater  jammers. 

4)  When  the  radar  pulse  consists  of  N  frequency  bursts,  the  receiver  can  be  channelized  with  N  channels 
where  each  channel  is  tuned  to  a  unique  frequency  and  has  a  bandwidth  equal  to  Af.  This  bandwidth  is  1/N 
times  the  total  signal  processing  bandwidth.  Therefore,  extremely  high  time-bandwidth  products  can  be 
achieved  utilizing  receiver  components  that  have  instantaneous  narrow-band  characteristics. 

5)  Channelization  allows  the  capability  for  selective  limiting.  In  a  conventional  FM  receiver  using  a  single 
limiter,  strong  CW  interference  anywhere  in  the  total  signal  bandwidth  is  likely  to  capture  the  FM  receiver. 
As  a  result,  the  target  return  signals  are  suppressed  and  unlikely  to  be  detected  with  a  preset  threshold 
crossing.  In  a  channelized  system,  CW  interference  suppresses  target  returns  in  only  one  of  the  channels. 
Hence,  the  target  signal  power  is  reduced  only  by  a  factor  of  [(N-l)/N]2.  If  N  is  greater  than  10,  CW 
interference  can  be  20  to  30  dB  above  the  per-channel  target  signal  power  without  preventing  target 
detection  in  a  fixed-threshold  system. 

6)  Modem  technology,  employing  phase-locked  voltage  controlled  oscillators  and  digital  frequency 
synthesizers,  allows  for  the  generation  of  sophisticated  stepped  frequency  waveforms.  In  fact,  specialized 
stepped-frequency  systems  have  been  developed  that  achieve  up  to  10  GHz  of  bandwidth  and  very  large 
time-bandwidth  products. 

4.2  The  Permutation  Matrix  (7] 

Frequency-hopped  waveforms  are  conveniently  represented  by  an  NxN  matrix,  denoted  by  A,  for  which  the  N 
rows  correspond  to  the  N  frequencies  and  the  N  columns  correspond  to  the  N  time  intervals  during  which  a  burst 
occurs.  The  entry  in  the  ith  row  and  j*  column,  designated  by  ay,  is  set  equal  to  unity  if  and  only  if  frequency,  fc+f|, 
is  transmitted  in  the  time  interval  t,  where  i,  j  =0,  1,  ....  (N-l).  Otherwise,  ay  is  set  equal  to  zero.  The  matrix.  A,  is 
known  as  the  permutation  matrix.  Because  the  elements  of  A  take  one  of  two  possible  values,  A  is,  also,  called  a 
binary  matrix. 
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The  binary  matrix  for  a  quantized  linear  FM  signal  is  shown  in  Figure  4.2.  For  this  case,  N=7.  A  dot  in  the 
ij‘h  block  indicates  ay=l  whereas  a  blank  indicates  a,j  =  0.  Note  that  there  is  only  one  dot  per  column  and  per  row. 
This  is  consistent  with  the  requirement  at  any  given  time  slot  that  only  one  frequency  is  transmitted  and  each 
frequency  is  transmitted  only  once. 
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Figure  4.2:  Permutation  matrix  for  a  stepped  frequency  linear  FM  waveform  has  nonzero  entries  only  along  one  of  the 

diagonals. 


An  N*N  matrix  contains  a  total  of  Q=N2  matrix  elements.  Thus,  when  the  elements  are  binary,  the  number 
of  different  matrices  possible  is  2Q.  However,  of  these  only  N!  matrices  can  be  obtained  by  permuting  the  N 
integers  contained  between  0  and  (N-l).  It  is  for  this  reason  that  A  is  referred  to  as  a  permutation  matrix.  Some  of 
these  are  better  than  others  for  designing  the  signal  pattern  of  a  radar  waveform. 

When  a  frequency-hopped  signal  is  reflected  from  a  target  and  received  by  an  observer,  it  is  shifted  in  time 
depending  upon  the  range  to  the  target  and  in  frequency  depending  upon  the  radial  velocity  of  the  target.  The 
amounts  of  these  shifts  are  determined  by  comparing  shifts  in  both  time  and  frequency  of  a  replica  of  the  transmitted 
waveform  with  the  actual  received  signal  and  noting  for  which  combination  of  time  and  frequency  shifts  the 
coincidence  is  greatest.  This  may  be  thought  of  as  counting  the  number  of  coincidences  between  l’s  in  the 
permutation  matrix  A  with  1  ’s  in  a  shifted  version  of  A  for  which  all  entries  have  been  shifted  r  units  to  the  right  (r 
is  negative  if  there  is  a  shift  to  the  left)  and  s  units  upward  (s  is  negative  if  there  is  a  shift  downward).  Let  the 
shifted  permutation  matrix  be  denoted  by  A(r,s). 

The  number  of  such  coincidences,  C(r,s),  is  the  unnormalized  autocorrelation  between  A  and  A(r,s).  C(r,s) 
clearly  satisfies  the  following  conditions: 

C(0,0)  =  N  (4.7) 

0<C(r,s)<N  ;  r*0  and  s?*0  (4.8) 

In  addition,  because  time  and  frequency  are  restricted  to  the  intervals  t0  <  t  <  t(N.i)  and  fc  <  f  <  fc  +  fcN.,),  it  follows 
that 


C(r,s)  =  0  ;  |rj>  N  and/or  js|  >  N  .  (4.9) 

The  two-dimensional  autocorrelation  function,  C(r,s),  is  called  the  ambiguity  function  and  can  be  thought 
of  as  the  total  coincidence  between  the  actual  returning  signal  and  the  transmitted  signal  shifted  by  r  units  in  time 
and  s  units  in  frequency.  It  is  useful  to  think  of  the  permutation  matrix,  A,  as  a  two-dimensional  template  of  N2  cells 
which  is  opaque  at  the  (N2-N)  cells  where  ajj=0  and  transparent  at  the  N  cells  where  a^l .  The  total  signal  energy 
emerging  through  the  N  transparent  windows  is  summed  to  give  the  value  of  C(r,s)  when  the  template  is  shifted  r 
units  on  the  time  axis  and  s  units  on  the  frequency  axis. 
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For  example,  with  regard  to  the  permutation  matrix  shown  in  Figure  4.2,  C(0,0)  =  7,  C(1 ,1)  =  6,  C(2,2)  =  5, 
C(3,3)  =  4,  C(4,4)  =  3,  C(5,5)  =  2,  C(6,6)  =  1  and  C(r,s)  =  0  for  all  other  choices  of  r  and  s.  This  yields  the  well 
known  diagonal  ridge  associated  with  the  linear  stepped  FM  signal.  In  the  real  world  the  returning  signal  is  always 
noisy.  For  large  values  of  N  and  a  very  noisy  environment,  this  signal  pattern  would  be  likely  to  produce  spurious 
targets  in  addition  to  the  ambiguities  caused  by  the  ridge. 

For  radar  applications  it  is  desirable  to  have  an  ambiguity  function  that  approaches  the  ideal  thumbtack.  In 
this  context  a  good  frequency  hopping  pattern  should  result  in  small  values  for  C(r,s)  when  (r,s)  *  (0,0).  Such 
patterns  have  been  proposed  by  John  Costas  [9]  and  are  discussed  in  the  next  section. 


4.3  Costas  Arrays  [9],  [10] 

John  P.  Costas  attempted  to  find  those  NxN  permutation  matrices  for  which 

max  C(r,s)  =  l  . 

(r,s)  *  (0,0) 


(4.10) 


Such  matrices  are  known  as  Costas  arrays.  In  addition  to  radar,  the  frequency  hop  patterns  associated  with  Costas 
arrays  may  be  useful  in  other  applications,  such  as  spread  spectrum  communication  systems,  where  the  objective 
may  be  to  achieve  either  jamming  resistance,  low  probability  of  intercept,  or  frequency  diversity  for  a  selectively 
fading  channel. 

The  problem  of  constructing  Costas  arrays  can  be  described  as:  Place  N  l’s  in  an  otherwise  null  N  by  N 
matrix  such  that  each  row  contains  a  single  1  as  does  each  column.  Make  the  placement  such  that  for  all  possible 
(r,s)  shift  combinations  of  the  resulting  permutation  matrix  relative  to  itself,  at  most  one  pair  of  1  ’s  will  coincide. 

Let  M(N)  denote  the  number  of  Costas  arrays  that  can  be  generated  from  the  N!  (NxN)  matrices  obtained 
by  permuting  the  N  integers  contained  between  0  and  N-l.  The  density  of  Costas  arrays  is  defined  to  be  M(N)/N!. 
These  values  are  tabulated  in  Table  4.1  for  3  <  N  <12.  Although  the  number  of  Costas  arrays  increases  with  N,  the 
density  decreases  significantly.  Thus,  a  random  search  by  means  of  a  computer  program  can  be  rather  lengthy. 

An  example  of  a  Costas  array  for  N=7  is  shown  in  Figure  4.3.  This  array  corresponds  to  the  particular 
sequence  of  ordered  integers  given  by 

{0n}={3,6,O,5,4,l,2}  (4-11) 

The  corresponding  ambiguity  matrix  is  shown  in  Figure  4.4  where  the  time  shift,  r,  and  the  frequency  shift,  s,  both 
extend  over  the  integers  in  the  range  from  -6  to  +6.  A  block  of  the  ambiguity  matrix  is  left  blank  when  C(r,s)=0. 
Note  that  C(0,0)=7  and  that  the  nonzero  values  of  C(r,s)  are  equal  to  unity. 


TABLE  4.1 

NUMBER  AND  DENSITY  OF  COSTAS  ARRAYS 


N 

3 

m 

5 

6 

Hi 

8 

9 

10 

11 

12 

M(N) 

4 

12 

40 

116 

200 

444 

760 

2,160 

4,368 

7,852 

M(N)|N! 

0.67 

0.5 

0.33 

0.16 

0.04 

0.011 

■ 
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Figure  43:  The  permutation  matrix  shown  is  one  of  200  possible  Costas  arrays  for  N=7. 


The  ambiguity  matrix  is  readily  constructed  from  the  difference  triangle  associated  with  the  sequence,  {0n} . 
By  definition,  the  element  in  the  i*  row  and  jth  column  of  the  difference  triangle  is  determined  by 


d>,j  “  ^(i+j-D  fyj-b 


(4.12) 


where  (i+j)  <  N. 


For  example,  consider  the  sequence  in  Equation  (4.1 1)  which  generates  the  Costas  array  in  Figure  4.3.  The 
difference  triangle  is  formed  by  successively  considering  i=l,2,3,4,5,6  and,  for  each  i,  letting  j  assume  all  values 
such  that  (i+j)  <  7.  Hence,  for  i=l  and  j=l, 2, 3, 4,5 ,6,  we  obtain 
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Figure  4:4:  The  ambiguity  matrix  is  shown  for  the  Costas  array  of  Figure  4.3. 
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dj  j  =  9l—0o  =6  —  3  =  3 
dj  2  =  62-9x  =  0  -  6  =  -6 
d13  =  03 -$2  =5-0  =  5 
di4  =  6>4  - =  4-5  =  -1 
dj5  =  05-04  =1-4  = -3 
d.,6  =O6-05  =2-1=1 


Similarly,  for  i=2  andj=l,2,3,4,5,  there  results 
d2l=92  -  90  =  0  -  3  =  -3 
d2 4  =  5  -  6  =  -1 

d2,3  =^4-^2  =4-0  =  4  (4.14) 

d2  4  =$5  —  #3  =1  —  5  =  —4 
d2,5=06~04=2-4  =  -2 


The  complete  difference  triangle  is  shown  in  Figure  4.5. 


It  can  be  shown  that  {0„}  yields  a  Costas  array  if  and  only  if  the  difference  triangle  formed  from  {0n}  has 
no  repeated  terms  in  any  row.  The  ambiguity  matrix  is  constructed  by  recognizing  that  d,j=k  indicates  that  a  1  is  to 
be  placed  in  the  block  of  the  ambiguity  matrix  corresponding  to  r=i  and  s=k.  For  example,  d,>2  =  -5  implies  that  a  1 
belongs  in  the  matrix  corresponding  to  r=4  and  s=  -5.  The  matrix  is  completed  by  recognizing  that 

C(r,s)  =  C(-r,-s)  .  (4.15) 

The  number  of  entries  in  the  difference  triangle  is  (N-l)N/2.  Utilizing  Equation  (4.15),  it  follows  that  there  are  (N- 
1)N  1  ’s  in  the  ambiguity  matrix. 


{  6n  }  =  { 3,  6,0,5,  4, 1,2} 
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Figure  4.5:  This  difference  triangle  can  be  used  to  construct  the  ambiguity  matrix  presented  in  Figure  4.4. 


The  ambiguity  matrix  provides  a  first-order  approximation  to  the  actual  ambiguity  function  because  it  does 
not  account  for  fractional  shifts.  A  detailed  evaluation  of  the  actual  ambiguity  function  for  the  frequency-hopped 
signal  is  likely  to  reveal  many  more  minor  peaks  and  a  few  peaks  whose  values  are  larger  than  unity. 
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The  frequency-hopped  signal  generated  from  a  Costas  array  is  referred  to  as  a  Costas  signal.  A  three- 
dimensional  view  of  the  ambiguity  function  for  the  Costas  signal  generated  from  the  Costas  array  of  Figure  4.3  is 
shown  in  Figure  4.6.  The  cut  along  the  time  shift  axis  for  zero  frequency  shift  is  shown  in  Figure  4.7.  Observe  that 
the  peak  value  of  the  ambiguity  function  has  been  normalized  to  unity. 


Figure  4.6:  The  ambiguity  function  of  the  Costas  signal  generated  from  the  Costas  array  of  Figure  4.3  approximates  the  shape  of 

a  thumbtack. 


The  plot  in  Figure  4.6  shows  a  few  peaks  reaching  a  value  of  0.29=2/7.  However,  most  peaks  reach  the  expected 
level  of  0.143=1/7.  The  ambiguity  matrix  does  not  predict  any  minor  peaks  along  the  time  shift  axis  for  zero 
frequency  shift.  Nevertheless,  many  minor  peaks,  most  of  which  are  less  than  or  equal  to  0.1,  are  revealed  in  Figure 
4.7. 


Figure  4.7:  The  ambiguity  matrix  does  not  predict  the  minor  peaks  along  the  time  shift  axis  for  zero  frequency  shift. 


In  summary,  the  Costas  array  generates  a  Costas  signal  whose  ambiguity  function  approaches  the  shape  of 
the  ideal  thumbtack.  The  ratio  between  the  main  lobe  and  the  pedestal  is  approximately  N.  However,  the  pedestal  is 
not  smooth.  The  ambiguity  function  has  a  few  peaks  higher  than  predicted  by  the  ambiguity  matrix  and  many  dips 
occur  at  various  locations  within  the  pedestal.  The  pedestal  stretches  from  the  origin  as  far  as  ±  NT  in  time  shift 
and,  essentially,  as  far  as  ±  N/T  in  frequency  shift.  The  number  of  major  sidelobe  peaks  predicted  from  the 
ambiguity  matrix  is  equal  to  (N-l)N.  There  are  approximately  three  times  as  many  blanks  than  1  ’s  in  the  ambiguity 
matrix. 
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4.4  Ambiguity  Function  of  the  Costas  Signal  [9] 

Without  loss  of  generality,  frequency-hopped  waveforms  can  be  analyzed  in  terms  of  baseband  signals  by 
considering  their  complex  envelopes.  The  complex  envelope  of  the  frequency-hopped  signal,  described  by 
Equations  (4.5)  and  (4.6)  and  normalized  to  have  unit  energy,  is  given  by 


where 


f  RAD  (0  — 


1 

Z  Pn(t-nT) 

n=0 

A  VNT 

Aej2*f"’  , 

0<t<T 

0  , 

otherwise 

(4.16) 


(4.17) 


By  definition,  the  energy  of  fRAo(t)  is 

NT  NT 

f  |W0f*=l  W‘)  CdM* 

0  0 

1  NT  N-l  N-l 

=  -i—r  I  I  P„(t-nT)  Z  Pm (t  —  mt) dt 

A  NT  0  n=°  m=0 

.  NT  rN_i  "I 

=  A*NT  l  {.?o|p"(t_nT|  }dt  <4',8) 

+  7^1  {l  Ip.(t-nT)  Pm(t-mT)|dt  . 

A  NT  0  [n=0  m=0  J 

n  #  m 


Note  that 

(,-„T,r=F- nTS,s(n+1)T 

[  0  ,  otherwise 

Also,  because  the  frequency  bursts  do  not  overlap  when  n*m, 

pn(t-nT)p*nl(t-mT)  =  0,  n*m. 


(4.19) 


(4.20) 


It  follows  that 


NT 

[  |fRAD(t)|  dt  =  1  .  (4.21) 


Therefore,  fRAD(t)  has  been  normalized  to  have  unit  energy. 


The  delay-Doppler  ambiguity  function  was  defined  in  Equation  (2.50)  to  be 

NT 

%(f,v)=  I  fiLu^c)  fRAD(<7'~T)ej2,ru<7  d<7  (2.50) 

0 
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where  the  infinite  limits  of  integration  were  truncated  because  of  the  finite  extent  of  fRAi>(t).  The  normalized 
ambiguity  function  is  known  to  have  the  following  properties: 


|^(T,U)|<l=|j(0,0)|  , 
\Z(~T,  -V)\  =  \Z(T,V)\  , 


I  /  \%(t,v)\z  drdv  =1 


(4.24) 


Because  of  the  odd  symmetry,  indicated  by  Equation  (4.23),  it  is  necessary  to  consider  only  nonnegative  values  of  t 
andu. 

Substitution  of  Equation  (4.16)  into  the  expression  for  the  ambiguity  function  results  in 

1  I  I  p'(cr-nT)  pm(<r-T-mT)  e**9a  da  .  (4.25) 

A  NT  n=°  m=° 

Separating  the  terms  in  the  double  summation  into  those  for  which  n=m  and  those  for  which  n*m,  we  obtain 


X(?,v)  =  I  X  p'n  (O-  -  nT )  pn  (a  —  T  —  nT )  e]2KVa  da 

A  N1  n=0 

+  ~r^r  J  X  X  P*n  (a  -  nT)  p  (a  -  t — mT)  eJ2lcvff  da 
A  NT  n=0  m= 0 

n*m 


Attention  is  first  focused  on  those  terms  in  Equation  (4.26)  for  which  n=m.  Make  the  change  of  variables 


p  =  a  -  nT  . 


(4.27) 


The  terms  for  which  n=m  then  become 


1—  £  ej2)r,'nT  J  p*n(/7)pn(p-f)ej2,r^  dp 
A  NT  n=0 

=  -L  S  ei""r 

AN  n=0 


(4.28) 
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where  the  autocorrelation  function  of  pn(t)  was  defined  to  be 


0n„(^)=^  f  P l(p)  Pn(p-'f)  ei2nVP  dp 


(4.29) 


Now  focus  on  those  terms  for  which  n*m.  As  before,  make  use  of  the  change  of  variables  in  Equation 
(4.27).  Those  terms  for  which  n*m  then  become 

—L-IIe'1""  1  p;(/»P„[p-r-(m-n)T]e1!'"’  dp 

A  NT  n=0  m=0 

n*m  (4.30) 

=  —7—  I’  I  0nm  [r  +  (m  -n)T,t>]  ej2*"  "T 

AN  n=° 

n*m 


where  the  cross-correlation  function  between  pn(t)  and  pm(t)  was  defined  to  be 


=  1  p'n(p)  P m(P~r)ei2’nJp  dp 

I  -°° 


(4.31) 


as 


Combining  the  results  obtained  for  n=m  and  n^^m,  the  ambiguity  function  of  Equation  (4.26)  is  expressed 


Z(T'V)=^ 


t  <pnm [r  +  (m-n)T,p] 

m-0 
n^m 


I  I 

n=0 


>  e 


j2m>nT 


(4.32) 


The  magnitude  of  the  ambiguity  function  represents  the  magnitude  of  the  coherent  processor  response  to  a  pulse 
train  arriving  with  delay,  x,  and  Doppler  shift,  u.  (This  assumes  that  the  Doppler  shift  is  approximately  unchanged 
over  the  different  frequency  hops  within  the  pulse.)  As  pointed  out  in  Section  4.3,  Costas  developed  his  approach  of 
frequency-hopped  waveforms  in  an  attempt  to  shape  the  ambiguity  function  in  the  form  of  a  thumbtack.  The 
success  of  this  approach  depends  upon  the  characteristics  of  (PnnCt,u)  and  (j>n,m(V)). 

An  analytical  expression  is  now  obtained  for  tpnn(t,v).  The  case  for  which  0  <  x  <  T  is  considered  first. 
Note  that  p n(p)  is  zero  outside  the  interval  [0,T]  while  p„(p-x)  is  zero  outside  the  interval  [x,  x  +T],  Since  x  >0, 
their  product  is  zero  outside  the  interval  [x,T],  Substitution  of  Equation  (4.17)  into  Equation  (4.29)  then  yields,  for  0 
<  x  <  T,  that 
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Finally,  since  both  p„(t)  and  pra(t)  are  zero  outside  the  interval  [0,T],  p*  (/?)  and  pm(p-T)  do  not  overlap  for  |t|  >  T. 
It  follows  that  cpnni(T,u)  is  zero  for  |x|  >  T.  Combining  the  above  results,  the  cross-correlation  function  between  p„(t) 
and  pm(t)  is  given  by 
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(4.39) 


^  2  C£  M)  (T  M)]  ^  -j[gg,„  (T+T)+2!T  f„  T)  U.I  <  j 

T  «nm  (T  -  |t|)  ’  11 

0  ,  |t|>T 

As  required,  note  that  <Pnm(t,t))  reduces  to  cpnn(t,u)  for  n=m. 

With  reference  to  Equation  (4.32),  it  is  now  shown  that  the  central  peak  of  the  ambiguity  function  is 
associated  with  the  summation  involving  (pnn(r,u)  while  the  sidelobe  responses  in  the  pedestal  of  the  ambiguity 
function  are,  in  the  main,  associated  with  the  summation  involving  cpnm[T+(m-n)T,u]. 

First,  consider  the  t  =  0  axis.  Along  this  axis  the  ambiguity  function  is  given  by 


Z(M)  =  — T77  2 
A~N  n=0 


N- 1 


(/>nn  (M)  +  Z  <t>nm  K™  “  *)T M  e 

m-  0 


jljrvnT 


For  n^m,  where  n  and  m  are  integers,  observe  that 

|(m  -  n)|  T  >  T  . 

From  Equation  (4.39),  it  is  concluded  that 

0-mKm-n)  T,i>]  =  0  ,  n*m  . 
Consequently,  Equation  (4.40)  simplifies  to 

X( 0,v)  =  “  I  K  (0,u)  ej2*u  nT  . 


AN  "=0 


However,  from  Equation  (4.35), 
<t>  m(0,v)  =  A 


2  sin(  n  v  T)  j!ruT 


K  V  T 


(4.40) 


(4.41) 


(4.42) 


(4.43) 


(4.44) 


It  follows  that 


Z(0,v) 


sin(  tv  v  T )  ^j2jn>nT 

71  V  NT  n=0 

g  jT  u  nt  sin  (k  v  NT  ) 
n  v  NT 


where  use  was  made  of  the  result  that 


N- 1 


I 

n=0 


jilt v NT  _  •. 

jlJtVnT  __  ^ 

~  ei2*vT  -  1 
eJKvNT  sin(  nv  NT) 

ejKVT  sin(^ur) 


(4.45) 


(4.46) 
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Observe  that 


*(0,0)  =  1  (4.47) 

and  that  the  first  zeros  of  %(0,u)  occur  for 

V  =  ±  —  .  (4.48) 

NT 

As  a  result,  the  Doppler  resolution  of  the  N-pulse  frequency-hopped  signal  is  identical  to  that  of  a  CW  pulse  of 
duration  NT. 


Costas  [9]  introduces  the  normalized  frequency  variable 

y  =  i)T  .  (4.49) 

Since  the  frequency  channel  spacing  is  1/T,  as  determined  by  Equation  (4.2),  y  is  a  measure  of  frequency  in  units  of 
frequency-hop  spacing.  The  ambiguity  function  along  the  x  =  0  axis  can  now  be  expressed  as 


*(0,  y)  =  ej*Ny 


sin(;r  Ny) 
n  N  y 


for  which  the  first  zeros  are  given  by 


(4.50) 


(4.51) 


The  ambiguity  function  along  the  u=0  axis  is  considered  next.  From  Equation  (4.32) 


*(T,0)  = 


1 


N-l 

I 

A"N  n=0 


<t>m  (r>°)  +  t  <t>nm  [r+(m-n)T,0] 

m=0 

n*m 


(4.52) 


The  summation  involving  the  cross-correlation  terms  can  no  longer  be  ignored.  However,  we  first  concentrate  on 
the  contribution  by  the  autocorrelation  terms  alone.  For  this  purpose,  define 


*'(*,  0) 


1 

a2n 


I  (r,0)  . 

n=0 


From  Equation  (4.35) 


g-j2irfnr 


M  -T 

H>t 


With  reference  to  Equation  (4.4),  it  follows  that 


(4.53) 


(4.54) 
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m  <t 


(4.55) 


x'(?>  0)  = 


(t-M) 

NT 


N-l 


I 


Because  the  sequence  of  integers  in  the  set,  {0n},  consists  of  all  integers  between  0  and  (N-l)  and  each  integer 
appears  only  once,  the  ordering  of  terms  in  the  summation  of  Equation  (4.55)  is  inconsequential.  Therefore,  the 
summation  can  be  written  as 


It  follows  that 


Z  e  T  =  Z  e  T 

n=0  n=0 
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-j2  7T-—T 


N 
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\  T  J 


%'(?$)  = 


.f  N'l 

•  ( N-0 
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l  T  J 

N  sin 


K  —  T 

V  T  j 


Observe  from  Equation  (4.54)  that  <pnr,(r,0)=0  for  |t|  >  T.  Hence, 

J'(T,0)  =  0  ,  \t\  >T  . 


(4.56) 


,  r<T  . 


(4-57) 


(4.58) 


It  is  seen  that  the  x-axis  response  beyond  |x|  =  T  is  entirely  due  to  the  summation  over  the  cross-correlation  terms  in 
Equation  (4.52). 

Costas  [9]  introduces  the  normalized  delay  variable 


x 


(4.59) 


Since  each  subpulse  is  of  width,  T,  x  is  a  measure  of  delay  in  units  of  subpulse  width.  Note  that  |  x  |  <T  implies  |x| 
<1 .  For  this  case  /'( x,0)  becomes 


rw>)  =  (i-|*|)  e"1'1""1’* 


sin(;rNx) 
N  sin(^x)  ’ 


(4.60) 


Observe  that 


r(o,o)  =  i 


(4.61) 


Since  Equation  (4.22)  is,  also,  valid,  it  can  be  deduced  from  Equation  (4.61)  that  the  cross-correlation  terms 
do  not  contribute  to  the  ambiguity  function  at  the  origin.  The  first  zeros  of  x'(  x,0)  to  either  side  of  the  origin  occur 
at 


or,  equivalently,  at 


N 


(4.62) 


(4.63) 


This  is  a  measure  of  the  delay  resolution  for  the  Costas  signal. 

It  is  of  interest  to  evaluate  the  Doppler  and  delay  resolutions  in  terms  of  the  time-bandwidth  product  of  the 
frequency-hopped  pulse.  The  duration  of  this  pulse  is 

Td  =  NT  .  (4.64) 

Because  each  subpulse  is  of  duration  T,  the  subpulse  bandwidth  is  approximately  1/T.  In  addition,  the  frequency 
hops  span  a  frequency  range  of  extent  equal  to  NAf  =  N/T.  Thus,  the  bandwidth  of  the  frequency-hopped  pulse  is 
approximately 


As  a  result,  the  time-bandwidth  product  of  the  Costas  signal  is  approximately 


Note  that 


while 


Td(BW)  =  (NT)  (y)  =  N1  . 

(BW)  1  _  1 
N2  ”  Td  ”  NT 


Td  =  1  =  T 
N2  BW  N 


(4.65) 


(4.66) 


(4.67) 


(4.68) 


Recall  that  the  Doppler  resolution,  from  Equation  (4.48),  is  1/(NT)  and  the  delay  resolution,  from  Equation  (4.63),  is 
T/N.  Therefore,  the  compression  ratio  for  the  Costas  signal  is  N2  in  both  the  Doppler  and  delay  domains. 

The  cross-correlation  terms  in  Equation  (4.52)  are  now  examined.  For  t  =  0  note  that 

9nm[t+(tn-n)T,0]  becomes  <Pnm[(m-n)T,0J.  Because  n^m,  )(m-n)T|  >  T.  It  follows  from  Equation  (4.39)  that. 
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0nJ(m-n)T>O]  =  O>  n*m 


(4.69) 


This  is  consistent  with  the  observation,  previously  pointed  out,  that  the  cross-correlation  terms  do  not  contribute  to 
the  central  peak  of  the  ambiguity  function  at  t  =  u  =  0. 

With  reference  to  Equation  (4.39),  note  that  (Pnm[?+(m-n)T,0]  is  not  identically  zero  only  over  the  range  of 
x  governed  by 


T+(m-n)T|  <  T 


(4.70) 


This  inequality  can  be  rewritten  as 


-T  <  r  +(m-n)T  <  T 


(4.71) 


or,  equivalently,  as 


-  T  -  (m-n)T  <  r  <  T  -  (m-n)T  .  (4.72) 


Hence,  along  the  t-axis,  (pnm[x+(m-n)T,0]  spans  an  interval  of  length  2T  centered  at  -(m-n)T.  Since  n^rri  in  the 
summation  over  m  in  Equation  (4.52),  the  interval  given  by  -  T  <  t  <  T  is  excluded.  Nevertheless,  for  a  fixed  value 
of  n  and  values  of  m  that  differ  by  unity,  the  corresponding  cross-correlation  terms  do  have  intervals  that  overlap. 
The  overlapping  of  these  terms  produces  reinforcement  and  cancellation  effects  which  are  quite  complex  because  of 
the  amplitude  and  phase  relationships  involved.  It  is  difficult  to  predict  the  exact  peak  values  of  the  sidelobes. 
However,  they  are  expected  to  be  on  the  order  of  1/N,  as  suggested  by  the  relevant  Costas  array. 

From  Equations  (4.4)  and  (4.37)  note  that 


or  =  f„  -f„ 


-v  = 


9„  -9„ 


-  v 


(4.73) 


For  |r+(m-n)T|<T  ,  use  of  Equation  (4.73)  in  Equation  (4.39)  yields 


|0mn[7+(m-n)T,  t>]|  =  A2 


(T-|f+(m-n)T|) 

T 


•  r  /  9 n  —  9  m 

-v)  (T  — |T+(m  —  n)T|)] 

n(0°~e*  - 

v)  (T-|r+(m-n)T|) 

(  T 

(4.74) 


Consider  the  case  for  which 


t  -  (n  —  m)T 

Then  the  magnitude  in  Equation  (4.74)  becomes 


(4.75) 
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where  dL,ra+i  is  the  difference  given  by  Equation  (4.12).  Consequently,  the  frequencies  at  which  the  cross¬ 
correlation  terms  peak,  for  a  specified  value  of  L,  can  be  obtained  from  the  difference  triangle  of  Section  4.3.  These 
peak  positions  are,  also,  reflected  in  the  corresponding  Costas  array.  Although  the  actual  peak  values  of  the 
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sidelobes  in  the  pedestal  are  difficult  to  predict,  they  are  expected  to  be  small  relative  to  the  central  peak.  This  is 
illustrated  in  the  following  examples. 

The  examples  are  based  on  an  N=10  and  an  N=30  Costas  array.  For  N=10,  the  firing  order  is  given  by  the 
sequence  of  ordered  integers  in  the  set 

{en}  =  {1,3,  7,  4,  9,  8,  6,  2,  5,0}.  (4.84) 

For  N=30,  the  firing  order  is  determined  from  the  set 


'2,  8,  26, 18,  25, 15, 16, 19,  28,  24, 12,  7,  23,  9,  29,  1 
27,  21,  3, 11,  4, 14, 13, 10, 1,  5, 17,  22,  6,  20,  0  j 


(4.85) 


They  are  referred  to  as  Welch-codes  because  they  can  be  constructed  from  a  theorem  proposed  by  L.  R.  Welch  [1 1], 

Various  views  of  the  ambiguity  surface  for  the  10-pulse  Welch  code  are  shown  in  Figure  4.8  [9]. 
jjjr(x,  y)|  is  plotted  where  x  and  y  are  the  normalized  variables  of  Equations  (4.59)  and  (4.49),  respectively.  In 

Figure  4.8(a)  observe  that  the  density  of  peaks  increases  towards  the  origin,  as  predicted  by  the  difference  triangle 
and  the  ambiguity  matrix  for  the  Costas  array.  A  different  perspective  of  the  ambiguity  surface  is  shown  in  Figure 
4.8(b).  Note  that  sidelobe  peaks  far  from  the  origin  remain  below  the  1/N  level.  However,  those  near  the  origin 
tend  to  have  their  values  above  the  1/N  level  due  to  the  overlapping  of  the  cross-correlation  terms.  Nevertheless, 
most  peaks  remain  below  the  2/N  level. 


(a) 


(b) 


Figure  4.8:  Various  views  of  the  ambiguity  surface  for  the  10-pulse  Welch  code  are  shown:  (a)  an  overview,  (b)  a  different 
perspective,  (c)  Doppler  response  for  x=0,  (d)  delay  response  for  y=0,  (e)  pedestal  region  near  the  highest  sidelobe  peak. 


The  Doppler  frequency  response  for  x=0  is  shown  in  Figure  4.8(c).  As  predicted  by  Equation  (4.50),  a 
|sin(/T  N  y)  /  7t  N  y  |  behavior  is  seen  with  the  zero  of  the  main  lobe  occurring  at  y=0.1=l/N.  The  delay  response 

for  y=0  is  shown  in  Figure  4.8(d).  From  x=0  to  x=l  the  |sin(;rNx)  /  N  sin(7T  x)|  behavior,  as  noted  in 
Equation  (4.60),  is  seen.  Beyond  x=l,  the  q>nn(x,0)  terms  vanish  and  only  the  9nm(x,0)  cross-correlation  terms 
contribute  to  the  ambiguity  function.  As  predicted  by  Equations  (4.75)  and  (4.80),  observe  that  |/Jf(x,0)|  =  0  for 
all  integer  values  of  x. 


The  largest  sidelobe  peak,  due  to  overlapping  of  cross-correlation  terms,  was  found  to  have  a  value  of  0.21. 
Because  N=10,  a  single  cross-correlation  term  would  result  in  a  sidelobe  peak  of  height  0.1.  Thus,  the  largest 
sidelobe  peak  has  a  value  which  is  2.1  times  larger  than  would  be  expected  from  a  single  cross-correlation  term. 
The  general  neighborhood  of  the  largest  sidelobe  peak  is  shown  in  Figure  4.8(e).  The  phasing  and  overlap  of  the 
cross-correlation  terms  are  seen  to  produce  isolated  peaks  approximately  twice  the  1/N  value  predicted  by  the  Costas 
array. 


Various  views  of  the  ambiguity  surface  for  the  30-pulse  Welch  code  are  shown  in  Figure  4.9  [9],  As  was 
done  with  Figure  4.8,  |)}f(x,  y)|  is  plotted  where  x  and  y  are  the  normalized  variables  of  Equations  (4.59)  and 

(4.49),  respectively.  Pedestal  spiking  above  1/N  is  evident  near  the  main  peak  of  Figure  4.9(a).  However,  away 
from  the  main  peak,  most  of  the  sidelobes  have  peak  values  approximately  equal  to  1/N.  The  largest  pedestal  peaks 
are  on  the  order  of  2/N.  Clearly,  the  pedestal  peaks  can  be  reduced  by  increasing  the  value  of  N. 


Figure  4.9(b)  shows  a  plot  of  |)jf(x,  y)|  where  values  are  evaluated  only  at  integer  values  of  x  and  y.  As 

would  be  expected  from  Equations  (4.75)  and  (4.80),  the  pedestal  region  of  the  ambiguity  surface  is  either  zero  or 
1/N  depending  upon  whether  or  not  there  is  a  cross-correlation  peak.  The  peaks  of  the  individual  cross-correlation 
terms  are  clearly  displayed  by  this  artifice.  However,  note  that  triangles  indicate  isolated  peaks  while  the  flat-topped 
sections  indicate  two  or  more  frequency-adjacent  peaks.  The  density  of  peaks  is  seen  to  increase  near  the  main  lobe, 
as  predicted  by  the  ambiguity  matrix  for  the  corresponding  Costas  array. 
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Figure  4.9:  Various  views  of  the  ambiguity  surface  for  the  30-pulse  Welch  code  are  shown:  (a)  an  overview,  (b)  plot  using  only 
integer  values  for  x  and  y,  (c)  Doppler  response  at  x=0,  (d)  delay  response  at  y=0,  (e)  complete  delay  axis  response,  (f)  Doppler 
response  at  x=25,  (g)  Doppler  response  at  x=5,  (h)  Doppler  response  at  x=1.4,  (i)  pedestal  region  near  the  highest  sidelobe  peak. 


The  Doppler  response  at  x=0  is  shown  in  Figure  4.9  (c)  for  values  of  y  between  0  and  1.  As  was  the  case 
with  Figure  4.8(c),  the  behavior  is  of  the  form  |sin(/T  N  y)  1 7t  N  y  | .  As  expected,  the  first  zero  of  the  main  lobe 

occurs  at  y=0.0333=l/30.  The  delay  axis  response  at  y=0  is  shown  in  Figure  4.9(d)  for  values  of  x  between  0  and  1 . 
As  with  Figure  4.8(d),  the  behavior  is  dictated  by  Equation  (4.60).  Consequently,  the  first  zero  appears  in  the 
vicinity  of  x=l/30.  The  complete  positive  delay  axis  for  y=0  is  shown  in  Figure  4.9(e).  The  main  peak  is  hard  to 
see  because  of  the  scale  chosen  for  the  x-axis.  However,  this  choice  of  scale  clearly  demonstrates  that  the  sidelobes 
are  very  well  behaved  with  relatively  few  excursions  above  the  1/N  value. 

Figure  4.9(f)  is  a  frequency  cut  taken  at  a  delay  of  x=25.  With  reference  to  the  set  of  ordered  integers  in 
Equation  (4.85),  the  25th  row  of  the  difference  triangle  is  given  by 

d25,j  =  ^25+j-l  0j-i  (4.86) 

where  j  assumes  all  positive  integer  values  such  that 

25  +  j  <  30  .  (4.87) 

As  a  result, 

d  25,1  =025-00=17-2=15 
d25,2  =026-01  =  22-8=14 

d25,3  =  021  -02  =  6  -  26=  -20  (4.88) 

d25,4  =028-03  =  20-18=2 
d25,5  =  029  -04  =  0  -  25=  -25  . 

Recall  from  the  discussion  of  the  difference  triangle  in  Section  4.3  that  dy=k  implies  a  peak  at  x=i  and  y=k. 
Therefore,  for  the  frequency  cut  at  x=25,  peaks  are  predicted  at  y=  -25,  -20,  2,  14,  and  15.  These  sidelobe  peaks  are 
clearly  seen  in  Figure  4.9(f).  Observe  that  the  isolated  cross-correlation  terms  have  their  peak  values  equal  to 
1/N=0.0333.  However,  the  two  peaks  at  y=14,15  experience  mutual  interference  effects  that  result  in  a  wider  peak 
whose  magnitude  exceeds  1/N  but  is  less  than  2/N. 

A  frequency  cut  for  x=5  is  shown  in  Figure  4.9(g).  Because  this  is  closer  to  the  main  peak  at  the  origin, 
many  more  sidelobes  appear  than  was  the  case  for  x=25.  Near  the  center  of  this  plot  the  mutual  interference  effects 
due  to  groups  of  frequency-adjacent  cross-correlation  peaks  are  clearly  evident.  The  largest  peak  values  are 
approximately  equal  to  2/N. 

The  largest  peak  in  the  pedestal  region  was  found  to  be  at  x=l  .4  and  had  a  value  equal  to  0.078  which  is 
2.34  times  the  1/N  value  of  0.0333.  Figure  4.9(h)  shows  a  frequency  cut  at  this  delay  with  the  largest  peak  occurring 
near  y=3.  The  ambiguity  surface  for  the  region  near  the  pedestal  highest  peak  is  shown  in  Figure  4.9(i).  The  rapid 
undulations  of  this  surface  are  to  be  compared  with  the  comparable  surface  in  Figure  4.8(e).  Note,  also,  the  lower 
level  of  the  surface  for  N=30  as  compared  to  that  for  N=10.  Increasing  N  from  10  to  30  has  increased  both  the  range 
and  Doppler  resolution  by  a  factor  of  3  and  has  lowered  the  pedestal  region  by  about  the  same  factor. 
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An  advantage  of  frequency-hopped  waveforms  is  their  flexibility  with  regard  to  ordering  of  the  frequency 
steps.  For  some  radar  applications  it  may  be  more  desirable  to  have  a  lower  pedestal  than  the  thumbtack  shape 
approximated  by  the  Costas  signal.  Figure  4.10  shows  the  ambiguity  surface  of  a  10-pulse  frequency  stepped  linear 
FM  signal  having  spacings  equal  to  1/T.  The  classic  FM  ridge  is  clearly  in  evidence  and  contains  most  of  the 
volume  under  the  surface.  Comparison  of  this  surface  with  that  shown  in  Figure  4.8(a),  which  was,  also,  for  a  10- 
pulse  burst,  reveals  the  much  lower  pedestal  of  the  linear  FM  waveform.  Although  this  is  an  advantage,  there  is 
potential  for  considerable  ambiguity  due  to  the  ridge  and  nearby  sidelobe  structure. 


-3.33  -1.37  0.00  1.67 

NORMALIZED  FREQUENCY 


Figure  4.10:  The  ambiguity  surface  is  shown  for  a  10-pulse  frequency  stepped  linear  FM  signal  with  spacing  equal  to  1/T. 
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CHAPTER  5 


AN  INTERFEROMETER  WAVEFORM  EMPLOYING  ORTHOGONAL 
FREQUENCY  DIVISION  MULTIPLEXING 


5.1  Orthogonal  Frequency  Division  Multiplexing  [12] 


5.1.1  Overview  of  OFDM 

Orthogonal  frequency  division  multiplexing  (OFDM)  is  a  modulation  technique  that  has  been  proposed 
recently  for  use  in  broad-band  mobile  multimedia  applications.  Because  frequency-selective  and  time-variant 
channels  are  typically  encountered  in  these  applications,  channel  equalization  must  be  incorporated  to  achieve  an 
acceptable  level  of  performance.  The  quality  of  service  is  usually  measured  in  terms  of  data  rate  and  bit  error  rate. 

Fading  channels  can  be  frequency  selective  when  multipath  propagation  is  experienced  and  can  become 
time-variant  when  the  transmitter  and/or  receiver  are  in  motion.  The  computational  requirements  for  an  equalizer  to 
be  used  to  compensate  for  a  frequency-selective  channel  increases  in  complexity  with  the  data  rate.  Furthermore, 
the  system  resources  (or  overhead)  needed  for  the  estimation  of  time-variant  channel  parameters  can  be 
considerable. 

To  overcome  these  drawbacks  OFDM  has  been  considered  for  high  data-rate  applications.  In  effect,  OFDM 
divides  the  radio  channel  into  many  narrowband  subchannels  which  appear  to  be  frequency  nonselective.  Thus,  the 
task  of  channel  equalization  is  simplified  to  estimating  a  single  complex  factor  (or  transfer  function)  for  each 
subcarrier.  Such  a  channel  estimation  can  be  realized  by  inserting  pilot  symbols  with  known  modulation  into  the 
transmitted  signal.  Based  on  these  pilot  symbols  the  receiver  can  measure  the  channel  transfer  functions  for  each 
subcarrier  using  interpolation  techniques.  In  this  case,  each  subcarrier  can  be  demodulated  coherently. 

Alternatively,  differential  modulation  can  be  performed  in  a  manner  similar  to  that  used  with  M-ary 
differential  phase  shift  keying  (M-DPSK).  With  OFDM  this  has  been  extended  for  multi-level  modulation  to  a 
combined  differential  amplitude  and  phase  modulation  in  order  to  achieve  a  higher  bandwidth  efficiency  and  still 
rather  good  performance.  This  modulation  technique  is  termed  differential  amplitude  and  phase  shift  keying 
(DAPSK).  The  performance  of  OFDM  systems  with  differential  modulation  compares  quite  well  with  systems 
using  nondifferential  modulation  and  coherent  demodulation.  Furthermore,  the  computational  complexity  in  the 
demodulation  process  is  quite  low  for  differential  modulation. 

Channel  coding  is  an  important  ingredient  for  OFDM  systems.  Transmission  over  a  frequency-selective 
channel  implies  that  some  of  the  subcarriers  are  strongly  attenuated  and  errors  occur  even  with  a  high  average  signal 
power.  In  this  flat  fading  situation  an  efficient  channel  coding  leads  to  a  very  high  coding  gain.  For  this  reason, 
OFDM  systems  will  always  make  use  of  channel  coding. 

Also,  apart  from  simplifying  the  channel  equalization  with  respect  to  single  carrier  systems,  OFDM  offers 
more  possibilities  to  adapt  the  transmission  system  to  the  communications  channel  by  varying  the  modulation  of  the 
subcarriers.  However,  an  optimal  adaptation  is  limited  to  quasistationary  channels. 

As  compared  to  single  carrier  systems,  the  most  important  advantage  of  the  OFDM  transmission  technique 
is  obtained  in  frequency  selective  channels.  The  signal  processing  in  the  receiver  is  rather  simple  in  this  case 
because,  after  transmission  over  the  communications  channel,  the  orthogonality  of  the  OFDM  subcarriers  is 
maintained  and  the  channel  interference  effect  is  reduced  to  a  multiplication  of  each  subcarrier  by  a  complex  factor. 
Therefore,  equalizing  the  signal  is  very  simple,  whereas  equalization  may  not  be  feasible  in  the  case  of  single  carrier 
transmission  covering  the  same  bandwidth. 
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5.1.2  OFDM  Signal  Characterization 

As  was  done  in  Section  4.4,  the  discussion  is  carried  out  at  baseband  through  the  use  of  complex  envelopes. 
OFDM  was  chosen  for  the  interferometer  signal  because  it  can  be  made  to  closely  resemble  the  Costas  signal 
selected  in  Chapter  4  for  the  radar  waveform.  This  should  greatly  enhance  the  ability  of  the  interferometer  signal  to 
mask  the  radar  signal  transmitted  through  the  sidelobes  of  the  main  radar. 

In  general,  an  OFDM  signal  consists  of  M  subcarriers,  each  of  duration  Ts  The  baseband  frequencies  of  the 
subcarriers  are  given  by 


m 


f  =  ^;m  =  0,l,  ...,  (M  —  1)  . 


(5.1) 


Consequently,  their  spacing  is 


A/  =  /m+1  ~fm  = 


m  + 1 


m 


1 


=  — ;  m=  0,  1,  (M - 1)  . 


(5.2) 


When  Ts  =T  and  M=N,  the  OFDM  signal  frequencies  are  identical  to  those  used  with  the  Costas  signal. 
The  mth  subcarrier  signal  is  described  analytically  by 


gjt) 


B  eJlKfm'  ,  0  <  t  <  Ts 

0  ,  otherwise 


(5.3) 


Each  subcarrier  can  be  modulated  independently  with  the  complex  modulation  symbol,  Sk>m,  where  k  refers  to  time 
and  m  refers  to  the  subcarrier  frequency.  Thus,  the  kth  OFDM  block  extending  over  the  interval  [k  Ts,  (k+l)TJ  is 
characterized  by 

M-l 

st(t)=  X  Sk>„  g„(t-kT,)  .  (5.4) 

m=0 

The  total  interferometer  signal,  consisting  of  K  OFDM  blocks  is  given  by 

K-l  K-l  M-l 

s,FM(t)  =  Z  sk(t)  =  I  X  Sk  gm(t -kTs)  .  (5.5) 

k=0  k— U  m— U 


The  energy  of  the  interferometer  signal  is 


kts  2  kts 

®ifm  =  ^  |sifm(^)|  dt  =  J  sIFM(t)  sIFM(t)  dt 
0  0 


(5.6) 


Because  the  OFDM  blocks  do  not  overlap  when  k*£, 

gm(t-kTs)g*n(t-£Ts)  =  0,  .  (5.7) 

Consequently,  the  expression  for  EffM  simplifies  to 
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5.13  OFDM  Signal  Recovery  Employing  a  Guard  Interval 

As  shown  in  Equation  (5.9),  all  subcarriers  are  mutually  orthogonal  over  the  time  intervals  given  by  [kTs, 
(k+l)Ts]  where  k=0,  1,  (K-l).  This  orthogonality  enables  the  employment  of  correlation  to  recover  the 

modulation  symbols,  Sk,m.  Specifically,  for  the  k*  OFDM  block,  note  that 
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K-l  K-l  M-l 

sifm(0  =  S,  8k(t)  =  sk>m  gm(t-kT)  .  (5.18) 

The  guard  interval  is  removed  at  the  receiver  and  only  the  time  intervals  [kT,  (k+l)T-Tg],  where  k=0,  1,  .... 
(K-l),  are  actually  processed.  As  a  result,  the  guard  interval  in  the  kth  block  is  a  pure  system  overhead.  Observe  that 
each  processing  interval  is  of  length 

(k  +  1)T  -  Tg  -kT  =  T  - Tg  =  Ts  .  (5.19) 

Consequently  the  orthogonality  of  Equation  (5.9)  is  maintained  and 

j  (k+l)T-T, 

Sk,m  =  l  sk(t)  g*m(t-kT)dt  .  (5.20) 

Assuming  the  guard  interval  length  is  larger  than  xm,  the  maximum  channel  delay,  no  intersymbol 
interference  will  occur.  This  is  illustrated  in  Figure  5.1 .  Note  that  the  message  portion  of  the  maximally  delayed  k* 
symbol  block,  arising  from  multipath,  does  not  overlap  the  message  portion  of  the  direct  path’s  (k+l)*  symbol  block 
since  Tg  >  Tm.  Thus,  there  is  no  intersymbol  interference  between  consecutive  symbol  blocks.  However,  the 
message  portion  of  the  maximally  delayed  k*  symbol  block  does  overlap  the  message  portion  of  the  k‘h  symbol 
block  in  the  direct  path.  This  results  in  fading  which  is  compensated  for  by  channel  equalization. 
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I  II  I 

|  (k-l)lh  Symbol  |  k1*1  Symbol  |  (k+l)lh  Symbol  | 

-  Block  |  Block  |  Block  j 

L.Guard  ^.Message  ^j^Guard  Message  Guard  .^Message 
I  Portion  |  Portion  |  Portion  |  Portion  |  Portion  |  Portion  | _ ^ 

(k-l)T-Ts  (k-l)T  kT  -  Tg  kT  (k+l)T-Tg  (k+l)T  (k+2)T-Tg  ^ 


(a) 


|  (k-1)*  Symbol  [  k*  Symbol  |  (k+l)tb  Symbol 

|  Block  |  Block  |  Block 

L  Guard  g|^  Message  Guard  gjg  Message  Guard  ^ 

I  Portion  |  Portion  |  Portion  |  Portion  |  Portion  | _ ^ 

(k-l)T-Tg  +  Tm  (k-l)T  +  Tm  kT-Tg  +  Xm  kT  +  Tm  (k+l)T-Tg  +  Tm  (k+l)T  +  Xm  t 


(b) 

Figure  5.1 :  The  guard  and  message  intervals  are  shown  for  (a)  the  direct  path  signal  and  (b)  the  maximally  delayed  multipath 

signal. 


5.1.4  Use  of  Fast  Fourier  Transform  in  OFDM  Implementation 

OFDM  is,  in  effect,  a  parallel  system.  A  principal  objection  to  the  use  of  parallel  systems  is  the  complexity 
of  the  equipment  (filters,  modulators,  etc.)  required  to  implement  the  system.  However,  the  modulation  and 
demodulation  processes  involved  with  OFDM  can  be  implemented  in  the  digital  signal  processing  parts  of  the 
transmitter  and  receiver  employing  Fast  Fourier  Transform  techniques.  This  is  discussed  next. 


From  Equations  (5.16)  and  (5.17)  the  kth  OFDM  block  is  given  by 


M-l 

SfcOO  =  BSk,m  e 


j2>rfm(t-kT) 


(5.21) 


Although  the  above  expression  is  valid  for 


kT-Tg  <  t  <  (k  +  l)T  -Tg  ,  (5-22) 

the  message  portion  extends  only  over  the  interval 

kT  <  t  <  (k  +  1)T  -  Tg  ,  (5-23) 

as  illustrated  in  Figure  5.1(a).  An  estimate  of  the  Nyquist  sampling  rate  for  the  kth  OFDM  block  is  obtained  by 
approximating  its  spectrum. 
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The  spectrum  of  the  m*  subcarrier  is 


Gm(f)  =  f  gm(t)  e-j2"ft  dt  =  /  Bej2"(f“-f),dt 


=  B 


ej2*(fn-f)TJ  _  e-j2*(f„-f)Tg 

j2x(fm-f) 


(5.24) 


=  BT  e 


j2jr(fm-f)  (-L-i)  sin[^T(f-fm)] 
*T(f-f m) 


Observe  that  the  first  zeros  in  the  spectral  main  lobe  occur  at  f  =  fm  ±  1/T  where  T  >  Ts.  Therefore,  the  width  of  the 
main  lobe  is  2/T.  As  noted  by  Equation  (5.2),  the  subcarriers  are  spaced  in  frequency  according  to  Af  =  1/TS. 
Consequently,  for  M  »  2,  the  bandwidth  of  each  OFDM  block  is  approximately 

BW  =  MAf  =  —  .  (5.25) 

Ts 

This  corresponds  to  the  bandwidth  of  the  bandpass  transmitted  signal.  The  baseband  bandwidth  is  half  of  that. 
Doubling  the  baseband  bandwidth  to  obtain  the  Nyquist  sampling  rate,  it  follows  that  the  Nyquist  rate  is 

f  =  BW  =  —  .  (5.26) 

Ts 

Let  sk(t)  be  sampled  only  over  the  message  portion  of  the  k*  OFDM  block.  Since  the  time  spacing  between 
the  samples  is 


1  T 

Ats=  7-  =TT  , 

f.  M 


(5.27) 


M  samples  are  required  to  cover  the  message  portion  of  the  k*  block.  The  sampling  instants  are  given  by 


tj  =kT  +  iAts  =kT  +  — ;  i  =0, 1, (M-l)  . 

M 

With  reference  to  Equation  (5.21),  the  i*  sample  of  sk(t)  becomes 

M-l  j2;rf„  |  kT  +  -  kT] 

o  \  =  c  _rtc  o  v  m  J 


Sk(^i)  —  Sk,i  —  ®  Sk.m  e 

ra=0 

M_1 

—  R  Y  <2  o 


(5.29) 


~  B  Z  Sk  e 

m=0 


M-l  j27T  — i 

=  B  Z  e  M 


It  follows  that  the  M  time  samples  (i=0,  1,  ...,  (M-l))  can  be  generated  by  performing  the  inverse  discrete  Fourier 
transform  (IDFT)  of  the  M  modulation  symbols  of  the  k*  OFDM  block,  Skm,  where  m=0,  1,  ...,  (M-l).  Of  course, 
fast  Fourier  transform  techniques  can  be  used  to  implement  the  IDFT.  This  provides  significant  simplification  for 
the  transmitter.  Given  the  M  time  samples,  sk(t)  is  readily  generated  using  the  usual  reconstruction  techniques  based 
upon  the  sampling  theorem.  sk(t)  is  then  used  to  modulate  a  carrier  waveform  at  the  desired  carrier  frequency. 
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At  the  receiver  the  modulated  carrier  is  converted  to  baseband  in  the  usual  way  to  produce  sk(t).  sk(t)  is 
then  sampled  to  generate  the  M  time  samples  skji  =  sk(tj)  where  i=0,  1,  (M-l).  The  message  symbols  are  then 

recovered  by  utilizing  the  discrete  Fourier  transform  (DFT)  such  that 


Jk,m 


1  M-l 

-  E 

B  >=° 


’k,i 


- \2x — m 
M 


(5.30) 


As  with  the  transmitter,  fast  Fourier  transform  techniques  can  be  used  to  implement  the  DFT.  This  considerably 
simplifies  implementation  of  the  receiver. 


When  the  subcarrier  spacing,  Af,  is  designed  to  be  much  smaller  than  the  channel  coherence  bandwidth  and 
the  symbol  duration  T  is  much  smaller  than  the  channel  coherence  time,  the  channel  transfer  function,  H(f,t),  can  be 
approximated  by  a  complex  constant  within  the  bandwidth  of  each  subcarrier  signal  and  the  duration  of  each 
modulation  symbol.  In  such  a  situation,  the  effect  of  the  radio  channel  is  a  simple  multiplication  of  each  subcarrier 
signal,  S^n,  gm(t-kT),  by  the  complex  transfer  constant,  =  H(mAf,  kT). 


As  a  result,  the  k,h  OFDM  block  results  in  the  received  signal 

M-l 

rk(t)  =  E  Skm  Hk  gm(t-kT)  +  nk(t)  (5.31) 

m=0  * 

where  nk(t)  denotes  the  additive  noise.  As  illustrated  in  Figure  5.1 ,  intersymbol  interference  is  avoided  provided  the 
guard  interval  length,  Tg,  is  larger  than  the  maximum  multipath  delay,  xm.  Therefore,  correlation  can  be  used  to 
recover  the  separate  symbols  associated  with  each  subcarrier.  Specifically, 

1  (k+l)T-Tg 

l  ri(<)  g'm(t-kT)dt.  (5.32) 

Alternatively,  the  DFT  could  be  used.  During  the  k*  block  time  interval, 

kT  -  Tg  <  t  <  (k  +  1)T  -  Tg  ,  (5.22) 

the  received  waveform  is 


M-l 

r,(t)  =  X  St  Hk  B  +  nk(t)  .  (5.33) 

m=0 

Sampling  rk(t)  at  the  M  sampling  instants  specified  by  Equation  (5.28),  yields  the  M  time  samples 

M-l  j2)T-i 

rk(t  j)  =  B  E  Sk  Hk  e  M  +  nk(t ,)  .  (5.34) 

m=0 


Application  of  the  DFT  to  these  M  time  samples  generates 


1  M-l  -j2)r— m 

Rlm  =  -  Z  rk(t ,)  e  M 
b  >=o 

=  ^k,m  ^k,m  +  ^k,m 


(5.35) 


where 
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(5.36) 


1  M-l  -'fin — m 

Nk!m  =  —  .2  nk(t  j)  e  M 
B  1=0 


denotes  the  scaled  additive  noise  Fourier  coefficients. 


The  above  detection  schemes  imply  perfect  synchronization  at  the  receiver.  Time  and  frequency 
synchronization  errors  disturb  the  orthogonality  of  the  subcarriers  which  can  considerably  impair  performance. 

A  general  block  diagram  of  the  OFDM  communication  system  discussed  above  is  shown  in  Figure  5.2. 
The  S/P  and  P/S  blocks  indicate  serial-to-parallel  and  parallel-to-serial  encoders,  respectively.  IFFT  and  FFT  refer 
to  the  inverse  fast  Fourier  and  fast  Fourier  transforms,  respectively.  The  bit  interleaver  is  used  to  scramble  the  bits 
so  as  to  mitigate  error  bursts.  In  addition,  convolutional  encoding,  which  is  readily  decoded  by  the  Viterbi 
algorithm,  is  performed  on  the  data. 
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Figure  5.2:  An  ODFM  communication  system  typically  employs  fast  Fourier  transform  techniques,  coding,  and  channel 

estimation. 


5.1.5  Modulation  and  Demodulation  Without  Differential  Encoding 

As  mentioned  in  Section  5.1.2,  each  subcarrier  can  be  modulated  independently  with  the  complex 
modulation  symbol,  Skim.  Generally,  this  complex  value  is  generated  using  multilevel  signaling  whereby  a  sequence 
of  p  binary  bits  is  mapped  into  one  of  2P  possible  symbols. 

If  no  differential  encoding  is  applied,  the  data  bits  to  be  transmitted  are  directly  mapped  to  the  modulation 
symbols,  Sk,m.  M-ary  phase-shift  keying  (MPSK)  and  M-aiy  quadrature  amplitude  modulation  (MQAM)  are 
examples  of  this  modulation  technique. 

Coherent  demodulation  has  to  be  applied  when  nondifferential  modulation  is  used  in  the  transmitter. 
Typically,  this  means  that  the  mixer  in  the  receiver  is  synchronized  in  both  frequency  and  phase  with  the  carrier 
frequency  of  the  received  signal.  However,  in  the  case  of  an  OFDM  transmission  system,  each  subcarrier  has  to  be 
synchronized  in  both  frequency  and  phase.  Also,  if  any  kind  of  amplitude  modulation,  such  as  MQAM,  is  used  in 
the  transmitter,  the  attenuation  of  each  subcarrier  must  be  known. 

To  generate  this  information  in  the  receiver,  channel  estimation  needs  to  be  performed  which  provides 
estimates,  Hkm,  of  the  channel  transfer  constants.  Given 


Rk.m  —  Sk,m  Hk  m  +  Nk>m 


(5.35) 
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the  decision  is  based  on  the  statistic 


r  R 

Dc  =  — 


k.m 


H 


k,m 


H 


k,m 


H 


Sk  + 

k.m 


N 


k.m 


k,m 


H 


k,m 


(5.37) 


Sk  »  an  estimate  of  the  transmitted  signal,  is  determined  by  noting  the  region  of  the  decision  space  containing 
Dck  m .  This  decision  procedure  is  indicated  by 

St.  =dec{DcM,}  .  (5.38) 


The  estimate,  Hk  m ,  is  obtained  by  using  known  pilot  symbols  which  are  inserted  into  the  OFDM  signal 

before  transmission.  The  first  step  for  the  receiver  is  to  extract  the  transfer  constants  at  those  times  and  frequencies 
at  which  the  pilot  symbols  were  transmitted.  Interpolation  is  then  used  to  estimate  the  transfer  constants 
corresponding  to  particular  values  of  k  and  m. 


5.1.6  Modulation  and  Demodulation  With  Differential  Encoding 

Differential  encoding  for  narrowband  channels  means  that  the  transmitted  information  is  contained  in  the 
quotient  of  two  successive  modulation  symbols.  Therefore,  the  data  bits  to  be  transmitted  are  not  mapped  directly 
to  the  modulation  symbols,  Skm.  Instead  the  quotient  of  two  successively  transmitted  modulation  symbols  is 
designed  to  contain  the  data  bit  information. 

This  technique  can  be  applied  either  with  respect  to  time  or  frequency.  In  the  time  direction  the  operation 
of  differential  encoding  is  described  analytically  by 

S|c,m  =Sk_,,m  x  Ck  m  (5.39) 


where  the  p  binary  bits  are  first  used  to  determine  Ck,m.  St™  is  then  determined  by  carrying  out  the  product  between 
Ckm  and  Sk.i  m.  Alternatively,  differential  encoding  in  the  frequency  direction  is  summarized  by  the  product 


Sk.m  —  Sk  m_i  x  Ck 


(5.40) 


As  before,  the  p  binary  bits  are  first  used  to  determine  Ckm.  Skm  is  then  calculated  from  the  product. 

In  both  cases  the  coherence  time  and  coherence  bandwidth  of  the  channel  need  to  be  large  compared  to  the 
symbol  duration  and  subcarrier  spacing,  respectively,  to  ensure  that  the  corresponding  channel  transfer  constants, 
Hkm,  are  adequate  measures  of  the  channel  behavior.  Cases  are  encountered  in  communication  systems  where 
differential  encoding  in  the  frequency  direction  has  some  advantages  with  respect  to  the  system  overhead  needed  to 
monitor  the  channel. 


The  topic  of  demodulation  is  now  discussed.  For  this  purpose,  only  differential  encoding  in  the  time 
direction  is  considered.  When  differential  encoding  is  employed  in  the  transmitter,  demodulation  at  the  receiver  can 
be  carried  out  noncoherently.  From  Equation  (5.35)  Rk,ra  and  Rk_i>m  are  given  by 


Rk,m  =  Sk,m  Hkm  +  Nkm  , 
Rk-l,m  =  Sk-l,m  Hk_,  m  +  Nk_,  m 
Employing  Equation  (5.39)  in  the  expression  for  R^,  there  results 

Rk,m  =^k-l,m  ^k,m  ^k,ni  +  ^k,m 


(5.41) 


(5.42) 


The  quotient  of  R^m  and  Rk_i,m  yields  the  statistic 
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An  estimate  of  the  message  symbol,  Ck,m,  is  obtained  by  noting  the  region  of  the  decision  space  containing  D™m . 
This  decision  procedure  is  indicated  by 

c„,m  =<JecRcJ  ■  (5.47) 

Note  that  channel  estimation  does  not  have  to  be  performed  with  the  noncoherent  demodulation  scheme 
discussed  above.  Thus,  the  computational  complexity  of  the  receiver  is  relatively  low.  Unfortunately,  the 
noncoherent  statistic  of  Equation  (5.46)  is  affected  by  twice  as  much  noise  power  as  the  coherent  statistic  of 
Equation  (5.37).  This  results  in  a  higher  bit  error  rate  than  does  the  coherent  demodulation  approach  with  good 
channel  state  information. 

Another  approach  to  detection,  called  quasi-coherent  demodulation,  involves  both  differential  encoding  and 
coherent  demodulation.  From  Equation  (5.39) 

Ck,m=^^  •  (5.48) 

Coherent  demodulation,  as  specified  in  Equation  (5.38),  is  used  to  obtain  estimates  of  Skm  and  Sk_,,  ra.  The  quasi- 
coherent  detection  statistic  is  then  given  by 
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(5.49) 


dec 


Dkm  = 

k,m 


Jk,m 


^-k,m  ^ 

Hk.m  ‘ 


^k-l,n 


dec 


^  k-l,n 


Hfc-l.n 


In  this  way,  the  channel  influence  is  removed  before  differential  decoding  takes  place.  As  with  the  other  detection 
procedures,  the  estimate  of  the  message  symbol,  Ck,m,  is  determined  by  noting  the  region  of  the  decision  space 

containing  D£cm  and  is  indicated  by  writing 

fy,  =dec{D'k;}  .  (5.50) 

Apart  from  differential  decoding,  the  processing  in  a  quasi-coherent  receiver  is  similar  to  that  in  a  coherent 
one.  However,  because  of  differential  encoding,  an  incorrect  decision  influences  two  successive  symbols.  Provided 
that  the  error  rate  is  relatively  small  (i.e.,  less  than  10%),  it  has  been  shown  that  the  error  probability  with  quasi- 
coherent  demodulation  is  approximately  twice  that  of  coherent  detection  with  perfect  channel  estimation. 


5.1.7  M-ary  Differential  Amplitude  and  Phase  Shift  Keying 

M-ary  differential  amplitude  and  phase  shift  keying  (M-DAPSK)  has  been  promoted  as  an  effective 
modulation  technique  for  achieving  higher  bandwidth  efficiency.  In  fact,  M-DAPSK  has  been  shown  to  have  a 
substantial  performance  improvement  over  M-ary  differential  phase  shift  keying  (M-DPSK)  for  M  >  16. 

Before  delving  into  M-DAPSK,  M-ary  amplitude  and  phase  shift  keying  (M-APSK)  is  discussed  first. 
Assume  p  input  bits  are  used  to  determine  each  symbol.  However,  let  pa  of  these  bits  determine  the  amplitude  of 
each  symbol  while  pp  of  these  bits  are  used  to  specify  the  symbol  phase.  Note  that 


P  =Pa  +Pp  > 

(5.51) 

the  number  of  possible  amplitudes  is 

Na  =  2Pa  , 

(5.52) 

the  number  of  possible  phases  is 

Np  =  2P>  , 

(5.53) 

and  the  total  number  of  symbols  (or  signal  states)  is 

Ms  =2P  =  2P°+Pp  =  Na  Np  . 

(5.54) 

As  an  example,  let  p  =  6,  pa  =  2,  and  pp  =  4.  This  results  in  4  amplitudes,  16  phases,  and  a  total  of  64 
symbols.  A  64-APSK  signal  space  is  shown  in  Figure  5.3.  The  signal  constellation  consists  of  four  different 
amplitude  rings,  each  with  16  phase  states.  The  amplitude  values  of  the  rings  are  1,  a,  a2,  and  a3.  Hence,  successive 
rings  have  amplitudes  which  differ  by  the  constant  factor,  a.  In  general,  the  signal  space  diagram  for  M-APSK  can 
be  described  by  the  signal  set 
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(5.55) 


w  = 


where  A  and  P  are  integers  selected  from  the  sets 

A  ={0,1,...,  (N.-l} 
P  =  {0,1,..„  (Np  — l} 

such  that  Na  and  Np  satisfy  Equation  (5.54). 


(5.56) 


M-ary  DAPSK  is  a  logical  extension  of  M-ary  APSK.  The  p  data  bits  for  M-DAPSK  first  determine  the 
message  symbol,  Ctm.  Let  A(Ck,m)  and  P(Ck)m)  denote  integers  selected  from  the  sets  A  and  P  of  Equation  (5.56). 
The  message  symbol  is  then  given  by 


f'  _  o  A(Ck,m) 
'-'k.m  —  a 


(5.57) 


Figure  5.3:  The  signal  space  contains  64  possible  symbols  when  there  are  4  possible  amplitudes  and  16  possible  phases. 


Differential  encoding  generates  the  transmitted  signal,  Sk,m,  according  to  Equation  (5.39).  Note  that  Sk.1;m  is  of  the 
form 


_  a  A(Sk_,,m) 

—  a 


P(Sk.liD,) 

.  N. 


(5.58) 


where  A(Sk.i>ra)  and  P(Sk.1>m)  are,  also,  integers  selected  from  the  sets  A  and  P.  When  Sk.i,m  and  Ck>m  are  multiplied, 
the  amplitude  exponents  are  added  as  are  the  phase  exponents.  However,  to  assure  that  S^n,  is  contained  in  the  set, 
*F,  it  is  necessary  to  add  A(Sk_iiin)  and  A(Ck,m)  modulo  Na.  Because  of  the  inherent  periodicity  in  the  phase,  it  is  not 
strictly  necessary  to  add  P(Sk_i,ra)  and  P(Ckim)  modulo  Np.  As  a  result,  the  transmitted  symbol  is 


C  _  „  [A(Sk-,,,n )  +  A(Ckim )  ]  mod  Na 

ak,m  —  a 


[p(sk-l,m)  +  P(Ck.m)] 

x  e  Np 


(5.59) 


By  way  of  example,  consider  the  signal  space  diagram  of  Figure  5.3  for  which  p  =  6,  pa  =  2,  and  pp  =  4. 
Assume  the  first  two  message  bits  are  used  to  determine  the  amplitude  of  Q,,,,  in  accordance  with 


00  — ►  1,  01  — ^  a,  11  — ^  a2,  10  — ^  a2  . 


(5.60) 
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The  remaining  four  bits  determine  the  phase  of  C^n,.  Table  5.1  shows  the  magnitude  of  Skm  as  a  function  of  the 
magnitude  of  Sk.iim  and  the  transformation  indicated  by  Equation  (5.60).  Observe  that  |Sk>m|  equals  |Ckm|  when 
|Sk-i,m|  is  unity. 

TABLE  5.1 


DIFFERENTIAL  AMPLITUDE  MODULATION  FOR  pa  =  2 


Sk,m 

1 

1 

a 

a2 

a3 

a 

a 

a2 

a3 

1 

- ? 

a2 

a3 

1 

a 

- 7 

a3 

1 

a 

a2 

AMPLITUDE 

BITS 

00 

01 

11 

10 

Given  a  value  of  p,  the  optimal  values  for  a  pa  and  pp  so  as  to  minimize  the  bit  error  rate  when  employing 
differential  encoding  are  given  in  Table  5.2.  When  p  <  4  (or,  equivalently,  Ms  <  16),  observe  that  differential  phase 
shift  keying  is  preferred  where  all  symbols  have  the  same  amplitude. 


TABLE  5.2 

OPTIMAL  MODULATION  PARAMETERS  FOR  DIFFERENTIAL  ENCODING 


p/Ms 

Pa^Na 

Pt/Np 

Noncoherent 

Demodulation 

a 

Quasi-Coherent 

Demodulation 

a 

<4/<l  6 

0/1 

p/Ms 

1.0 

1.0 

4/16 

1/2 

3/8 

2.0 

1.8 

5/32 

1/2 

4/16 

1.6 

1.45 

6/  64 

2/4 

4/16 

1.4 

1.38 

7/128 

2/4 

5/32 

1.3 

1.21 

It  is  possible  to  eliminate  the  modulo  operation  on  the  sum  of  the  amplitude  exponents  by,  in  effect, 
integrating  it  into  the  amplitude  mapping.  For  this  purpose,  let  the  transmitted  symbol  be  given  by 


:,k,m 


=  s 


k-l,m 


■/k,m 


(5.61) 


Define  the  extended  set  of  integers  given  by 
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A'  =  {(-N,+l),  (-N.+2),  ....  -1,  0,  1,  ....  (N.-l)}  . 


(5.62) 


Let  AXC'k.m)  be  an  integer  selected  from  A'  and,  as  before,  let  P(C'kini)  be  an  integer  selected  from  the  set,  P.  C'^n,  is 
then  defined  by 


r'  -  aA'(c'*.n.>  eJNPP<C|t,m) 

L'k,m  -  a  e 


(5.63) 


The  integer  selected  for  A^C'^)  depends  both  upon  the  pa  data  bits  and  |Sk.i>m|  .  It  is  chosen  in  such  a 
manner  that 


|c  _  [A(St_i,„)  +  A'(C't>m) 

Pk,m  “ 


(5.64) 


is  either  1,  a,  a2,  ...,  a(Na'2),  a(Na'n.  By  way  of  example,  consider,  once  again,  the  signal  space  diagram  of  Figure  5.3. 
Table  5.3  shows  the  magnitudes  of  C'^n,  and  as  a  function  of  the  magnitude  of  Sk.]im  and  the  amplitude  data 
bits.  Observe  that 


r'  = 

''•'k.m 


K„, 

Nk-l,i 


(5.65) 


as  was  the  case  with  modulo  arithmetic.  However  the  amplitude  bits  are  now  decoded  by  recognizing  that 


(5.66) 


lcU = 1 

-y 

00 

|C  ,,m  |  =  a 

or 

j/3 

/a 

->  01 

|c;,„|  =  a! 

or 

y2 

/a 

— >  1 1 

|C'k,„|  =  a3 

or 

1/ 

/  a 

->10 

TABLE  5.3 

DIFFERENTIAL  AMPLITUDE  MODULATION  FOR  pa  =  2 


lsk-l,m 

|Ck,m  | 

1 

1 

a 

a2 

a2 

1 

a 

a2 

a2 

a 

1 

a 

a2 

1/a 

a 

a2 

a2 

1 

? 

1 

a 

1/a2 

1/a 

a2 

a3 

1 

a 

a2 

1 

l/aJ 

1/a2 

1/a 

a2 

1 

a 

a2 

AMPLITUDE 

BITS 

00 

01 

11 

10 

00 

01 

11 

10 
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Table  5.4  summarizes  the  mapping  specified  in  Equation  (5.66)  and  gives  thresholds  that  would  be 
employed  in  the  detection  process.  Due  to  the  division  of  two  successively  received  noisy  symbols,  the  exact 
thresholds  for  noncoherent  demodulation  are  hard  to  calculate.  However,  a  good  approximation  is  given  by  the 
geometric  average  of  two  adjacent  valid  amplitudes.  As  a  result,  the  thresholds  in  Table  5.4  are  of  the  form 

Tanc  =  aj  VI;  j=(— Na  +  1),  0,  (Na  -2)  .  (5.67) 


The  thresholds  used  to  obtain  the  phase  information  bits  are  given  by  half  the  distance  between  two  adjacent  valid 
phases.  Hence, 


rpllC 


-  +  1 


\  7jr 

i  =  0, 1, (N 

NP 


-1) 


J 


(5.68) 


TABLE  5.4 

THRESHOLDING  SCHEME  FOR  EVALUATION  OF  AMPLITUDE 
INFORMATION  BITS 


c'  1 

k,m  | 

I5 

“I5- 

T1 

1 

a 

AMPLITUDE 

BITS 

01 

11 

10 

00 

01 

11 

10 

a'25 

a” 

a0i 

ao.s 

a15 

a2i 

THRESHOLDS 

a25 

a'5 

a05 

a05 

a'5 

a25 

In  the  case  of  quasi-coherent  demodulation  with  independent  identically  distributed  symbols,  the  optimum 
thresholds  are  identical  to  those  used  with  coherent  demodulation.  In  particular,  they  are  half  the  distance  between 
two  adjacent  valid  phases  or  amplitudes.  Thus, 

Tq°  =TC  =aj  0±^)  ;  j=0jl>  _  (Na  -  2) 

(\  \  ?7T  (5'69) 

Tpqc=TpC=Tpnc  =  ^  +  iJ  —  ;  i  =  0, 1,  (Np  —  1)  . 

As  indicated  in  Figure  5.2,  channel  coding  should  be  implemented  with  OFDM.  The  radio  channel 
attenuates  each  OFDM  subcarrier  by  the  complex  transfer  constant,  Hk>m.  If  the  channel  is  a  multipath  channel  with 
many  propagation  paths  and  without  a  line-of-sight  path,  the  amplitude  of  the  transfer  constants  is  Rayleigh.  This 
means  that,  even  at  a  very  large  average  signal-to-noise  ratio,  there  are  always  some  subcarriers  that  are  strongly 
attenuated  and  have  a  rather  low  signal-to-noise  ratio.  Hence,  many  bit  errors  would  occur  on  these  subcarriers.  For 
this  reason,  channel  coding  is  a  very  important  component  of  OFDM  systems.  Soft  decision  decoding  has  been 
shown  to  be  superior  to  hard  decision  decoding.  A  soft  decision  algorithm  first  decides  the  result  based  upon  the 
statistic  being  above  or  below  a  decision  threshold  and  then  gives  a  “confidence”  number  that  specifies  how  close 
the  text  statistic  is  to  the  threshold  value.  In  hard  decisions  only  the  decision  output  is  known.  The  Viterbi  decoding 
algorithm  examines  the  possible  paths  through  a  code  tree  and  selects  the  best  one  based  on  some  conditional 
probabilities. 
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5.2  Ambiguity  Function  of  the  OFDM  Signal 


The  delay-Doppler  ambiguity  function  of  the  interferometer  waveform  is  defined  in  Equation  (2.61)  as 

T(t,v)  =  ]  f,;M  (<J)  f[FM(£T-f)  ej2™CT  d(7  .  (2.61) 

— oo 

Because  of  the  odd  symmetry  for  an  ambiguity  function,  indicated  by  Equation  (4.23),  it  is  necessary  to  consider 
only  nonnegative  values  of  z  and  v. 

Having  chosen  the  interferometer  signal  to  be  the  OFDM  signal  specified  by  Equation  (5.13)  and  (5.3),  the 
ambiguity  function  becomes 

1  “  fK-l  M-l  1 

KMTsB2(Ts2  -  lk=°  m=°  J 

{K-l  M-l  1 

XZ  S,n  gn  (cr  -T-  ^Ts)|  ej2-CT  dtr 

Interchanging  the  order  of  summation  and  integration  results  in 


1  K-l  M-l  K-l  M-l 

r(r,D)= - - — j — r-  X  X  X  X  skm  sln 

KMT  B2<t2  k=0  m=0  *=°  n=0  ' 

s  s  (5.71 

J  gl  (o’  -kTs)  gn  (<r  -T-  £TS )  ei2>rvcr  da 

— oo 

The  four-fold  summation  in  Equation  (5.71)  can  be  grouped  into  terms  for  four  cases: 

1)  k  =  £,  m  =  n,  2)  k  =  £,  m  *  n,  3)  k  *  £ ,  m  =  n,  and  4)  k  *  C,  m  *  n.  This  yields 

K-l  M- 1  .,  “ 

r(T,l»)  KMT.B'o]  =  X  X  S,,„  I  gm(o-kT,)g„  (o -X -  W e>M°  da 

k=0  m= 0  1  1  -oo 


K-l 

M-l 

M-l 

+ 

X 

k=0 

X 

m= 0 

X 

n-  0 

s;,m 

m±n 

K-l 

M-l 

K-l 

+ 

X 

X 

X 

c  < 

k= 0 

m= 0 

e=o 

k*t 

K-l 

M-l 

K-l 

M-l 

+ 

X 

X 

X 

X 

s; 

k= 0 

m- 0 

(=0 

n= 0 

(5.72) 
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Recall  that  the  modulation  symbols,  Sk,ra  are  selected  randomly  depending  upon  the  data  to  be  transmitted. 
Consequently,  T(t  ,  v)  is  a  random  function  of  the  modulation  symbols.  Assuming  statistical  independence  from 
one  symbol  to  another  and  stationarity,  the  statistical  average  of  Equation  (5.72)  involves  expectations  of  the  form 

E[|Si.mf]  •  ElSi.m]  •  a”d  E[Sk,J  •  From  Equation  (5. 11) 

Et|®k,mf]  =<TS2  •  «.ll) 

In  addition,  assuming  circular  symmetry  for  the  symbol  constellation,  such  as  shown  in  Figure  5.3, 

e[s;J=  b[sJ  =  0  .  (5.73) 

It  follows  that 

E[s;,„  Sj=  e[s'Je[sJ=  0;m*n 

E[s;,»  S  J  =  E[s;.J  e[s,  J  =  0  ;  k*  t  (5.74) 

e[s;,  s,  J  =  E[s;,m  ]  e[s(i„  ]  =  0  ;  k  *  l,  m  *n  . 

Consequently,  taking  the  expected  value  of  Equation  (5.72)  results  in 


E[r(r,u)]  = 


1 


M-l 


X  x  f  g*m(tT-kTs)gm(fT-T-kTs)ej2^  d<7  .  (5.75) 


KMTSB  k=°  m=o 

The  expression  in  Equation  5.75  is  further  simplified  by  making  the  change  of  variables 

P  =  a  -  kTs  .  (5.76) 


Then 


E[r  (t,v)]  = 


i 


K-l  M-l  ” 


S  X  i  gKp)g„0?-T)e'2"l',*‘T'>d^ 


KMTSB  k=0  m=0  — 

1  K-l  M-l 

X  X  ej2™kT»  y/mm  (r ,  v) 


(5.77) 


KMB2  k=0  m=0 

where  the  autocorrelation  function  of  gm(t)  was  defined  to  be 

1  g’m(p)gm(p-T)ei2,rvp  dp  . 

lc 


(5.78) 


Note  that  Equations  (5.3)  and  (5.78)  are  identical  in  form  to  Equations  (4.17)  and  (4.29),  respectively.  It  follows 
from  Equation  (4.35)  that 


^mm^)  =  j 


B2 


(T,-H)  sin[;ri>(T, -|r|)] 


v(T,+r)-2^f„r] 


XV(T-\t\) 

0 


m  <  t 


T  >  Tu 


(5.79) 
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Observe  that 

(^„(0,0)=B!  .  (5.80) 

Therefore,  as  would  be  expected  due  to  our  normalization, 

E[r(0,0)]  =  1  .  (5.81) 

Although  the  terms  corresponding  to  the  expectations  in  Equation  (5.74)  average  to  zero,  they  contribute  to 
significant  fluctuations  of  the  ambiguity  function  over  the  x-v  plane.  A  measure  of  these  fluctuations  is  the  variance 
of  r(x ,  v).  By  definition,  this  is  given  by 


Var{r(T,i>)}=  E[|r(r,u|2j-  |E[r(r,w]| 2  . 


E[F(x  ,  v)]  is  found  from  Equation  (5.77).  Therefore,  it  remains  to  determine  E[|  T(x  ,v)| 2],  For  this  purpose,  it  is 
convenient  to  return  to  Equation  (5.71). 


Making  the  change  of  variables 


in  Equation  (5.71),  it  follows  that 


P  =  o-  kTs 


I  g„  (<r  -  kT.)  g,  (<7  -  r  -  <T.)  ia 

— oo 

=  e«”kT’  f  g.(p)e.  IP  -r  ~(<-k)T,]  e'2"'  dp 


=  T.  e 


j2flrukT. 


Wmn  +(^~k)Ts  ,  V  ] 


where  the  cross-correlation  function  between  gn,(t)  and  gn(t)  is  defined  to  be 

Wmn  (T’V)  =  I  gl  (P)gn(P~T)  Q*2KVp  dp  .  (5.84) 

Ts  — 

Note  that  Equations  (5.3)  and  (5.84)  are  identical  in  form  to  Equations  (4.17)  and  (4.31),  respectively.  It  follows 
from  Equation  (4.39)  that 


t(T,-M)  ~M)] 

Ts  namn(l%-\r\) 


,-j[«,«mn(T,  +  r)+  2^fnT] 


r  <T„ 


A  -  Ts 


where  otmn  is  given  by 


Or  =  f  -  f  -v 

mn  m  n 
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Consequently,  Equation  (5.71)  can  be  expressed  as 

r(T,u)  = 


1  K-l  M-l  K-l  M-l  _ 

1  XXX 


KMB2  6 2  k=0  m=0  £=0  n=0 

It  follows  that 

\T(t,v)\2  =T(f,V)  r*( T,v) 


(5.86) 


f 


1 


\2 


\KMB2  g\  j 


K-l  M-l  K-l  M-l 


III  I  S;jmS^ej2"ukT^mn^+(^-k)Ts,u] 

l,_A  C  —  A  r\ — a  *  ’ 


k=0  m=0  e=0  n=0 


K-l  M-l  K-l  M-l  .  ,.T  , 

X  X  X  X  X  Sk...S  .„.  e-,1**kT^;v[T+(<'-k')T1,u] 

k  =0  m  =0  1=0  n  =0  ’  ’ 


(5.87) 


(  i  V 


VKMB  2alj 


K-l  M-l  K-l  M-l  K-l  M-l  K-l  M-! 


i  i  i  i  i  s  i  i  su  s<n  sk.  m.  s;v 

_ ft  A  e A  — A  U  A  »V.  — A  C  — A  n  — A  ’  ’ 


k=0  m=0  t=0  n=0  k'=0  m=0  1=0  n=0 


x  eJ!«k-k'»T'  *r„[T+«-k)T„»]^;...[f+(<'-k,)T.>l/l 

The  expected  value  of  I  E(x  ,v)  1 2  is,  therefore, 

/ 

1 


e[|T(t,u)|j]  = 


V 


KMB2  o]  j 


K-l  M-l  K-l  M-l  K-l  M-l  K-l  M-l 

ZSZZ  SIZE 

k=0  m=0  1=0  n=0  k'=0  m'=0  ('=0  n'=0 


•pfc*  C  <2  O*  1  p  jitfulk-k  )TS 

^  L^k.m  ^k'.m'  'V,n'  J  e 


(5.88) 


^m„^+(^-k)Ts,v]^v[T+(r-k')Ts,u]  . 

Recall  that  the  modulation  symbols,  Skjm,  are  assumed  to  be  stationary  and  statistically  independent  with 
zero  means.  Consequently,  when  E  [s’k  m  n  Sk-  m-  S*-  n- j  can  be  decomposed  into  a  product  for  which  one  or 

more  of  the  factors  is  the  expectation  of  a  single  symbol  or  its  conjugate,  the  expectation  in  Equation  (5.88)  is  zero 
and  the  corresponding  terms  can  be  excluded  from  the  eight-fold  summation. 


In  addition,  for  some  choices  of  the  indices  in  Equation  (5.88),  E|Sk  m  S,  n  Sk<  j  can  be 

decomposed  into  a  product  involving  factors  of  the  form  E[Sk  m]  and  E[(Sk  m)  j  .  For  circularly  symmetric 
signal  sets  these  expectations  are,  also,  zero  as  explained  next. 
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(5.89) 


In  general,  is  a  complex  quantity  which  can  be  expressed  in  polar  form  as 


Sk,m  —  Akm  e 


Consequently, 


SL  =  .  (5.90) 

Assuming  pp  data  bits  are  used  to  specify  the  symbol,  the  number  of  possible  phases  is  given  by  Equation  (5.53). 
When  the  symbol  constellation  has  circular  symmetry,  such  as  that  in  Figure  5.3,  the  phases  are  distributed  by  equal 
angular  increments  where  the  increment  is  given  by 


2k 

A</>  =  — —  . 

2Pp 

Beginning  with  symbols  placed  on  the  real  positive  axis,  the  symbol  angles  are  given  by 

In  Y 

kA0  =  — —  i  k=0,l, ... ,  (N  -1)  . 

2Pp 

A  total  increment  of  n  corresponds  to 

i  it  2nk 

kA0  = -  =  n 

2Pp 


(5.91) 


This  requires  that 


r  i  Nd 

k  =  2  p  =  — 2- 


which  is  an  integer  in  the  set  for  k  since  Np  is  guaranteed  to  be  an  even  integer.  Therefore,  for  each  symbol  in  the 
set,  there  is  a  corresponding  symbol  with  the  same  amplitude  but  displaced  in  phase  by  n  radians.  This  justifies  the 
conclusion  that 


E[ej<?k  m  J  =  0  . 


(5.95) 


Now  consider  Sk  m  for  which  the  angle  is  given  by  20km  .  The  symbols  corresponding  to  Sk  m  are 
distributed  at  equal  angular  increments  given  by 


2  A<p  =  2 


2  K  _  2  K  _  2  K 

2pp  2(Pp-1)  2Pp 


Pp  =Pp“1 
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In  effect,  Sk  m  is  distributed  in  the  same  manner  as  Sk,m  except  that  the  number  of  data  bits  used  to  specify  the  phase 
of  S2k  m  is  reduced  by  unity.  Therefore,  the  number  of  distinct  phases  for  Sk  m  is  half  of  those  for  S^.  Based  on 
these  observations,  it  is  concluded  that 


E[ej20kra  J  =  0  . 

Because  the  amplitude  data  bits  are  statistically  independent  of  the  phase  data  bits,  it  follows  that 

Efe,  ]  =  Ek,„  e12  ““  J  =  e[a;,„  ]  E^1'- 1=0  . 

Also,  since  the  complex  conjugate  of  zero  is  zero, 

eI(s;J2|=o  . 


(5.98) 


(5.100) 


Excluding  those  terms  for  which  E[s*k  m  S,  n  Sk>.  S’.  n-]  is  zero,  the  only  remaining  terms  in  the  eight¬ 
fold  summation  of  Equation  (5.88)  are  those  having  the  form 

i  |4  I  I  |2  |  |2 

E  S,  I  and  E  IS,  IS.  n  I  where  k*t  and  m?m.  These  arise  in  the  following  ways: 

L  K,m  I  J  L|  j  |  t,n  | 

1)  Let  k  =  l  -  k'  =  €'  and  m  =  n  =  m'  =  n'. 


The  corresponding  terms  in  Equation  (5.88)  are  given  by 


kmbV.2 


K-l  M-l  [j  .4]  .  . 


(5.101) 


2)  Let  k  =  £,  m  =  n  and  k'  =  C',  m'  =  n'  but  exclude  k  =  l  =  k'  -  l'  and  m  -  n  -  m'  -  n'.  The 
corresponding  terms  in  Equation  (5.88)  are  given  by 


K-l  M-l  K-l  M-l  ,2 1  C  .2 

EES  Z  E  Skm  E  Sk.J 

lr=(l  tm=0  k'iO  m  =f)  ’  1  ’  * 


KMB  2<J2  J  k=0  m=0  k'=0  m’=0 

k^k'  &  m±m' 


(5.102) 


x  ej2^(k-k')is  ¥mm[T,v]wmW[r,v\ 


3)  Let  k  =  k',  m  =  m’  and  £  =  €',  n  =  n'  but  exclude  k  =  £  =  k'  =  £'  and  m  =  n  =  m'  =  n'.  The 
corresponding  terms  in  Equation  (5.88)  are  given  by 
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( 


V 


vKMB2<72y 


K.-1 

M-l  K-l 

M-l 

s 

z  z 

Z  E 

k„i  J 

E 

|s,.|  J 

k=0 

m=0  ^=0 

o 

II 

c 

k*£  &  m*n 


Letting 


^4  =E 


Jk,m 


and  making  use  of  Equations  (5.1 1),  (5.101),  (5.102),  and  (5.103),  Equation  (5.88)  becomes 


£[|r(r,u)|!]  = 


'  1  V 


KMB1  a)  j 


+ 


1  V  K- 1  M-l  K- 1  M-l  „  „ 

1  X  £  £  £  e>MU-t<T-ylmmV,v\YMliM 

k-0  m- 0  *  =0  m  =0 

k±k'& 


KMB 


1  A2  AT-l  M-l  K- 1  M-l  , 

-  £.5  ,5  £  lr-fr+«-*)r.,»i| 

'  k*l  &  m*n 


KMB 


With  reference  to  Equation  (5.82),  it  is,  also,  necessary  to  evaluate  I  E[  T(x  ,v)]  1 2. 
follows  that 


\e[y{t,v)]2  =£[rMMrMF 


AT-l  M-l  K- 1  M-l 


Z  Z  Z  Z  e“^i% nC„. n , u ] 

KMB2) 

r  12 


'  1 

KMB  J  k=0  m= 0 
2 


/  i  V  *-i  M-i  k-i  M-i  , 

kmb  ;  *-•  — 


Combining  the  results  for  E[|  f(x  ,v)| 2]  and  I  E[T(x  ,v)]| 2  in  Equation  (5.82)  yields 


(5.103) 


(5.104) 


(5.105) 


From  Equation  (5.77)  it 


(5.106) 
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(5.107) 


Var{T(T,  v)}  = 


1  1 


K  MB  cr 


/  V  M-I 


K-l  M-l  K-l  M-l 


22  k_[»+«-k)T..»]| 


1C  MR 2  /  k=0  m=0  £=0  m=0 

y  &m?tn 


An  obvious  choice  for  Ts  is 


Ts  =  T 


(5.108) 


where  T  is  the  duration  of  a  subpulse  in  the  Costas  single  proposed  for  the  radar  pulse  waveform.  Also,  let 


K  =  M  =  N 


(5.109) 


where  N  is  the  number  of  subpulses  in  the  Costas  single.  Then  the  frequencies  of  the  Costas  and  OFDM  signals 
become  identical  and  the  duration  of  the  OFDM  signal  equals  that  of  the  radar  pulse. 

Note  that  Equation  (5.77)  can  be  expressed  as 


1  K-l  M-i  .  _ 

E[T(t,v)]  =  T'(t,v)  + - 7  Z  Z  ysmm(T,v)  ej2™kTs 

L  v  ’  v  7  KMB  k=°  m=0 


(5.110) 


T,(r,v)  =  ——T  Z  ysmm(T,v)  ej2™mT‘ 
KMB  m=° 


(5.111) 


while  Equation  (4.32)  can  be  expressed  as 


Z  Z  <f>nm[T+(m-n)T,v)  ej2nvnT 

A  fy  »=0  m= 0 


(5.112) 


where 


X\w)  =  i  <pmm  (T,v)  e 
A  N  m=0 


j2>ri)mT 


(5.113) 


As  pointed  out  in  Chapter  4,  the  central  peak  of  the  Costas  signal  ambiguity  function  is  associated  with  %' (f  ,V) . 
Under  the  conditions  of  Equations  (5.108)  and  (5.109) 
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(5.114) 


g2 

¥mm(W)  =  — p  0mm  (t,v) 

Then  T'(z ,v)  becomes  equal  to  %' (T ,V)  within  a  scale  factor.  This,  in  turn,  implies  that  the  central  peaks  of 
the  ambiguity  functions  for  the  OFDM  and  Costas  signals  have  identical  shapes  with  identical  resolutions  in  delay 
and  Doppler.  Because  the  non-cooperative  bistatic  radar  does  not  know  a  priori  in  which  direction  the  main  radar  is 
pointing,  it  must  perform  the  difficult  task  of  pulse  chasing.  However,  it  is  likely  to  be  unsuccessful  in  completing 
this  task  because  it  will  be  unable  to  distinguish  between  reflections  involving  either  the  OFDM  or  Costas  signal.  In 
addition,  it  will  be  difficult  for  the  non-cooperative  bistatic  receiver  to  separate  the  OFDM  and  Costas  signals 
because  they  contain  the  same  frequencies  and,  assuming  time  synchronization,  will  experience  amplitude 
transitions  at  the  same  time  instants.  As  a  result,  performance  of  the  non-cooperative  bistatic  radar  will  be  severely 
degraded.  This  issue  is  investigated  further  in  Chapter  6. 


Ill 


CHAPTER  6 


NON-COOPERATIVE  BISTATIC  RADAR  RECEIVER  PERFORMANCE 


6.1  Cross  Ambiguity  Function  for  the  Costas  and  OFDM  Signals 

As  discussed  in  Section  2.3,  the  non-cooperative  bistatic  radar  receiver  performs  detection  by  utilizing  an 
estimate  of  the  main  radar  transmitted  signal  for  correlation  with  the  received  data.  In  this  report  it  is  assumed  that 
the  estimate  is  obtained  from  the  direct  path  signal  observed  by  the  bistatic  radar.  As  revealed  by  Equation  (2.63), 
the  detection  statistic  used  by  the  non-cooperative  bistatic  radar  receiver  involves  the  cross  ambiguity  function, 
Z(t,u),  between  the  Costas  and  OFDM  signal  in  addition  to  their  delay-Doppler  ambiguity  functions.  An  expression 
for  Z(r,u)  is  derived  in  this  section. 

With  reference  to  Equation  (2.54),  the  cross-ambiguity  function  is  defined  to  be 


Z(T,V)=  !  f,*FM  (cKrad (<*-*)  ej2*u<r  d<J  ■ 


(2.54) 


Substitution  of  Equations  (4.16)  and  (5.13)  into  Equation  (2.54)  yields 
1  ~  (k-i  m-i 


Z(T,V)  = 


AB<7s  ^kmntj  - 


1  X  X  s;,m  g*m(0--kTs) 


k=0  m=0 
N-l 

Z 

n=0 


(6.1) 


Z  pn  (<7  -  t  -  nT)  >  ei2irvcr  da 


Once  again,  assume  the  conditions  of  Equations  (5.108)  and  (5.109).  Interchanging  the  order  of  summation  and 
integration  yields 


Z(T,V)  = 


1 


N-l  N-l  N-I 

z  z  z  s: 


ab<tstnVn  kr°  m“  “°  k,m 

2  g;(a-kT)  pn(cr-T-nT)  ei2*vff  da  . 

Making  the  change  of  variables 

p  =  a  -  kT 


(6.2) 


(6.3) 


in  Equation  (6.2),  it  follows  that 

I  gi>-kT)  p„(<7  -r-nT)  eilnvcr  dcr 

—  OO 

=  ei!'”lT  ?  g;  (P)p.[p-T-(n-k)T]  ip  (6.4) 

=  T  ei2"‘Te„„[5-+(n-k)T,p] 


where  the  cross-correlation  function  between  gn,(t)  and  pn(t)  is  defined  to  be 
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(6.5) 


0m„(^)=7f  I  gm(P)Pn(P-‘T)  ej2?rUP  dp  . 

1  -°° 

Utilizing  the  same  procedure  employed  to  develop  Equation  (4.39),  it  follows  that 


fAD(T-M>  sin[g«„(T-|r|)] 


(*.*>)  = 


0 


)[jianm(7+T)+2x{DT] 


H  -T 
H  >T 


(6.6) 


where  dm,,  is  given  by 


=  f  m  -  f  n  "  ^ 


Consequently,  Equation  (6.2)  can  be  expressed  as 

1 


Z(T,V)  = 


ABcr  nVn  k=°  m=0  n=0 


S  I.  Z  S:mej2jrukT^mn[f+(n-k)T,y]  . 


(6.7) 


(6.8) 


Because  the  modulation  symbols,  Sk,m,  are  selected  randomly  based  upon  the  data  to  be  transmitted,  Z(t,u) 
is  a  random  function  of  the  modulation  symbols.  Note  that 

E[Z(T,t»)]  =  s  X  X  5S;,Jel"’9,[It(n-k)T,»]  =  «  (6.9) 

ABcrsNVN  k=0  m=0  n=0 

where  use  was  made  of  equation  (5.73).  Thus,  as  seen  from  Equation  (2.63),  Z(x,u)  behaves  as  a  zero-mean  additive 
noise  component  of  the  detection  statistic  for  the  non-cooperative  bistatic  radar. 


6.2  Non-Cooperative  Bistatic  Radar  Detection  Statistic  Involving  the  Costas  and  OFDM  Signals 

As  expressed  in  Equation  (2.44)  the  detection  statistic  for  the  non-cooperative  bistatic  receiver  is  the  sum  of 
nine  terms  where  the  first  four  terms  are  given  by  Equations  (2.51),  (2.58),  (2.59),  and  (2.60),  respectively.  The 
remaining  five  terms  involve  receiver  noise  and  are  developed  next  by  carrying  out  the  products  indicated  in 
Equation  (2.43). 

Making  use  of  Equations  (2.3),  (2.38),  (2.39),  and  (2.40),  these  five  terms  can  be  expressed  as 


T  _  i_TP  -j2>r  [fDP  rT?  +  fBr+ur'IP  -  vr] 

A5  —  KRAD  e 

I  Wt-rw)  [nDP(t-rB  +4  e™'""*  dt , 

T  _  1.TP  -j2;r  [fDPrTP  +  fBr+urTP  -  vt] 

*6  ~  KIFM  C 

i  f,F„(t-TTP)[nDf(t-rB+r)r  dt , 

— oo 

T  _  „-j2^[-fBTTP  +fBT+urTT  -VT] 

1?  —  c 

J  nTP(t)  [nDP(t-TB  +  r)]*  dt , 


(6.10) 


(6.11) 


(6.12) 
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(6.13) 


i  _  (.DP  V  p.j2>r[fTPrTP  -  fTPr-urTP  +  vr] 

*8  “  \KRAD)  e 

In^Wf’  (t-T^+T)  e-i2'<'”-*  dt , 


[fleT  -f‘rT-VT,r  +VT] 


T  _  (kDP  V  12*  [ 
x9  VMFM/  c 


ln"(t)f'  (t-f'+f)  dt  . 


(6.14) 


The  statistic  can  be  written  as 


where  the  signal  portion  of  the  statistic  is 


e  ~es  + 


(6.15) 


and  the  interference  component  of  the  statistic  is 

lx  =  i  Ii  .  (6.17) 

1  i=2 

Recall  from  Equation  (2.39)  and  (2.40)  that  t  and  t>  denote  the  errors  experienced  in  estimating  tb  and  fB, 
respectively.  Perfect  estimates  are  obtained  when 

T  =  V  =0  .  (6.18) 

In  order  to  evaluate  the  best  possible  performance  of  the  non-cooperative  bistatic  radar,  it  is  assumed  throughout  the 
remainder  of  this  chapter  that  Equation  (6.18)  is  valid.  The  expressions  for  the  nine  terms  of  the  detection  statistic 
then  simplify  to 


I,  =  k,  ^f(0,0) , 

12  =  k2  Z*(0,0)  , 

13  =  k3  Z(0,0) , 

14  =k4r(0,0), 


(6.21) 

(6.22) 


I5  =kfAD  e~j2*fDP,rTP  J  U-Tff)  [nDP(t-rB)]*  ej2<rfD,‘  dt  , 

— oo 

I,  =  C,  [n”'(t-TB)I  e*""'  dt  , 


I,  =  j  nTP(t)  dt  , 


(6.25) 


Is  =  te)’  I  nTP(t)  e-^fT?‘  dt  , 
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dt 


(6.27) 


I»  =  (CJ  =|!,,'’'i.T’(l)fra(.-t")  e-2"”' 

— OO 

In  general,  the  detection  statistic  is  a  complex  random  variable.  However,  because  the  Costas  signal  is 
deterministic,  x(0,0)  is  a  constant.  Utilizing  Equation  (4.47),  it  follows  that 

E[/.]  =  E[iJ  =  k,  .  (6.28) 

Thus,  the  signal  portion  of  the  detection  statistic  is  a  constant. 

With  respect  to  the  interference  component  of  the  detection  statistic,  note  that 

E[I2]  =  E[I,1  =  0  (6.29) 

due  to  Equation  (6.9).  Also,  with  reference  to  Equation  (5.81) 

e[i4]  =  k4  E[r(0,0)]  =  k4  .  (6.30) 

Finally,  assuming  the  complex  envelopes  of  the  receiver  noises  in  the  direct  and  total  paths  to  be  statistically 
independent  with  zero  means,  it  follows  that 

E[IS]  =  E[lJ  =  E[I,]=  E[I,]=  E[l,]  =  0  .  (6.31) 

Consequently, 

E[/i]  =  k4  .  (6.32) 

Hence,  the  interference  component  of  the  detection  statistic  is  a  random  variable  with  a  nonzero  mean. 

The  variance  of  is  now  derived.  By  definition, 

Var{^1}=  e[|(,|]2  -|E[«,ir  ■  (6-33) 

Note  that 

\L\2  =£.  £\  =  i  £  LI*  .  (6.34) 

1  11  1  1  i=2  j=2  1  1 

Therefore, 

E fki|2]=  £  £  e[lI*]  =  £  E [ll. |2 ]  +  £  £  e[l  1*1  .  (6.35) 

L|  1 1  i=2  j=2  1  i=2  L|  1 1  J  i=2  j=2  1  1 1 

i*j 

Observe  that  the  double  summation  in  Equation  (6.35)  contains  64  expectations  that  need  to  be  evaluated. 

To  evaluate  the  expectations  involving  I2  and  I3,  it  is  necessary  to  obtain  an  expression  for  Z(0,0).  From 
Equation  (6.8) 

*  •~r~~~F7  S  s'  Z  S^»„[(n-k)T,0]  (6,36) 

AB(TsNVN  k=0  m=0  n=0 

where  0ran(r,O)  is  given  by  Equation  (6.6).  Since  0ra„(T,  u)  is  zero  for  |t|  >  T  ,  it  follows  that 

0mn  [(n-k)T,0]  =  0  for  n^k  .  (6.37) 
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As  a  result,  Z(0,0)  becomes 


(6.38) 


However, 


<U»,o)  = 


AB  ,  m=k 
0  ,  m^k 


Thus,  the  expression  for  Z(0,0)  simplifies  to 

^  ^  £  s»  ■  (640) 

E  [  1 1, 1 2  ]  and  e[|I,|  2  ]  are  evaluated  first.  Note  that 

E[|l,|f  =  |k2|J  e[|Z(0,0)|2]  .  (6.41) 

Utilizing  Equation  (6.40),  it  follows  that 

Elinor)  =  ^  £  I  e[s;,s,J 

f  S  (6.42) 

=  —TT  I  I  E  |skk|2  +x  i  E [s*  k ]  e[s,J  • 

/t2N3  k=0  L'  k’k'  J  k=0  (=0  L  K’KJ  e’eJ 

k*i 

V  X 

where  the  statistical  independence  between  different  data  symbols  was  invoked.  However,  the  data  symbols  are 

2 

random  variables  with  zero  mean  and  variance  O'. 


Consequently, 


e[|z(o,o)|2]= A- 


(6.43) 


ferl-l 


(6.44) 
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In  a  similar  manner, 


_  |k3 


N2 


Also, 


E[i!i;]=k!k;  E{[z'(0,0)f} 


(6.45) 


(6.46) 


Use  of  Equations  (6.40)  and  (5.99)  and  the  statistical  independence  of  the  zero  mean  data  symbols  results  in 


p  [r 

T*  1  _  k2k3 

N-l 

y 

V  E[s  s  1 

£|1 

2  3j  (Js2N3 

k=0 

ML 

(t'N3 


jl  I  e[sJ  e[s„] 


k=o  r=o 
k*; 


=  0 


Similarly, 


e[i,i;J  =  0 


Next  observe  that 


E[z(o,o)]  = 


1  N-1 

cr  nVn  k=° 


I  e[s;J  =  o 


due  to  the  zero  mean  value  of  the  data  symbols.  It  is  concluded  that 

e[i2i;J=e[i4i;J=e[i2i;J=e[i5i;J 


=  e[i2 

i;1=e[ 

I,  i'J=e[ 

i2i;1=e[i 

=  e[i, 

i;]=e[ 

i.«;]=e[: 

[,  i;]=e[i 

si;] 

=  e[i, 

i;]=e[ 

i,  i;]=e[: 

i2i;1=e[i 

,i;l 

(6.47) 


(6.48) 


(6.49) 


(6.50) 


Also,  because  the  receiver  noises  are  zero  mean  and  statistically  independent,  it  follows  that 
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E  k  /61 =  E  k  II  ]  =  E  [/,  /;  ] : =  E  [i,  /•  ] 

=  e  [/3  /,*]=£[/«/;]=£  [/,  i;]=e  [i,  i; } 
=  e  [/,  /;]=£■[/,/;]=£  [/,  /;]=£[/«  /;  ] 
=  e  [/4  /;]=£[/,/;]=£  [/,  /;]=£[/,/;  ] 

=  e  [?,  /,*]=£  [/,  /;  ] = e  [i,  i;]=e  [/,  /;  ] 

= 

=  e  [/« /;]=£  [/» /;]=£[/,/;]=£  [/,  i;  ] 

=  £[/,/»•]=£[/,/;]=  0 


Finally,  from  Equation  (5.13),  it  is  noted  that 

l  K-l  K-l 


E  [fffM  0)]  “ 


£  E[st,Jg.(t-kT.)=  0 


(6.52) 


(6.53) 


B<xsA/KMTs  k=°  m=° 

due  to  the  zero  mean  value  of  the  data  symbols.  Thus, 

eIisi;I=e[i6i;]=e[i,i;1=e[i,i;J=  o 

At  this  point  58  of  the  64  expectations  in  Equation  (6.35)  have  been  evaluated.  It  remains  to  evaluate 

E |l 4 1 2 J»  e|i5|2J.  e|i,|2 j,  e|i7|2 j,  e|ib|2J,  and  e|i9|2  . 

Note  from  Equation  (6.22)  that 

Ellul]2  =  |k4|  E[|r(o,o)|2]  .  (6.54) 

An  expression  for  E[|r(T, l>)|  2  ]  is  given  by  Equation  (5.105).  Observing  from  Equation  (5.79)  that 


^mm(0,0)  =  B2 

and  from  Equation  (5.85)  that 

^m„(0>0)  =  0;  m*n, 
the  expression  for  E[|r(0,0)|  2]  simplifies  to 


(6.55) 

(6.56) 
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:[|r(o,o)|!]= 


N-I 

//4n  z  b4 

m=0 


n2b2 


N-l  N— 1  N-I  N-l 

XXX  X  B4 

k=0  k'=0  m=0  m'=0 

k*k'  m*m' 


where  the  conditions  ofEquations  (5.108)  and  (5.109)  have  been  employed. 


(6.57) 


With  respect  to  Equation  (6.57),  observe  that  there  are  N  terms  in  the  single  summation  involving  the  index 
m.  In  addition,  there  are  (N2  -  N)  terms  in  the  double  summation  involving  the  indexes  k  and  k'  where  k  *  k'  and 
there  are  (N2  -  N)  terms  in  the  double  summation  involving  the  indexes  m  and  m’  where  m  *  m’.  Consequently, 
there  are 

(n2  -  n)  (n2  -  n)  =  N2  (N  - 1)2  (6.58) 


terms  in  the  four- fold  summation  involving  the  indexes  k,  k',  m,  m'  where  k  *  k'  and  m  *  m'.  It  follows  that 


E[|r<o,o)|J] 


N2 


^-  +  (N-l)2 


(6.59) 


Thus, 


[ll4f]=^ 

+  (N-1)2 

l<*i  J 

An  expression  for  I5  is  given  by  Equation  (6.23).  Therefore, 

e[|i, n  f >-0 

E{[nDt(t  -tb)T  nDP(//  itiM 


(6.60) 


(6.61) 


Denote  the  autocorrelation  function  of  nDP(t)  by  RnDP  (f)  where 

RnDP(r)  =  E{[nDP(t)]’  nDP(t  +r) }  .  (6.62) 

Consequently, 

E{[nDP(t-TB)]*  nDP(// -rB)}  =  RnDP(//-t)  .  (6.63) 

Assume  that  nDP(t)  can  be  approximated  as  white  noise  with  a  two-sided  power  spectral  density  given  by  N0.  Then 

RnDP(r)  =  N0£(T)  (6.64) 


where  8(t)  is  the  Dirac  delta  function.  Hence, 
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*[w!]=te,f  /™>('  -  T”) 


N0S  (//—/)  e 


jlKfDP(t-H) 


dtdjU 


(6.65) 


-  /fcrp  jy 
~  \krad\  -,vo 


since,  by  Equation  (4.21),  fRAD(t)  has  unit  energy. 


An  expression  for  I6  is  given  by  Equation  (6.24).  As  a  result, 


4w2]=l*£.f  ii 

*  I  I  — oo 

£{[«"VOr 

-tef  n  £k*(r  -!■”•)  /;„(/<  -  ol 

I  I  — oo 


(6.66) 


N0  S(fii-t)  e 


jilt  fDe  (t~H) 


dtdfi 


--KVk  1  E[\u„(t-t^]dt 


-  lit77*  N 

~  \^IFM\  ■'v0 

where  use  was  made  of  Equation  (6.64)  and  the  fact  that  fffM(t)  is  normalized  by  Equation  (5.13)  to  have  unit 
average  energy. 


An  expression  for  I7  is  given  by  Equation  (6.25).  It  follows  that 

J.[|/,|2]= jj  k'oo]-} 

dtdfi 

The  autocorrelation  function  of  n^ft)  is  defined  to  be 

RnTp(T>  =  e{[  nTP  (t)]*  nTP(t  +  f)}  . 

Making  use  of  Equations  (6.63)  and  (6.68),  Equation  (6.67)  becomes 


(6.67) 


(6.68) 
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(6.69) 


E 


/J  RnTP  (t  — //)  Rn„(ji-t)  ej2*f®(//  0  dtd n 


At  this  point  it  would  be  convenient  to  approximate  both  nDP(t)  and  nTP(t)  as  statistically  independent  white 
noise  processes.  However,  this  leads  to  the  integrand  of  Equation  (6.69)  containing  a  squared  impulse  with  the  end 

result  being  that  e[|  I7  |  2  J  is  infinite.  This  unrealistic  result  is  caused  by  the  simplifying  but  idealized  assumption 

of  white  noise. 

In  actuality,  the  non-cooperative  bistatic  radar  will  have  a  receiver  bandwidth  of  2W.  A  realistic 
assumption  is  that  the  power  spectral  density  of  the  receiver  noise  is  constant  over  this  bandwidth.  Therefore,  it  is 
assumed  that  the  power  spectral  density  of  the  noise  complex  envelopes  is  rectangular  with  height,  N0,  and 
extending  from  -W  to  +W  Hz.  The  autocorrelation  function  of  this  bandlimited  noise  is  then  given  by 

,  .  ,  i*  _  y  *  , »  sin  2 tc  IV f 

R  DP(r)  =  R  ^(r)  =  f  N0  eJ  df  =  2N0W  —  — —  .  (6.70) 

"  n  -w  2kWt 


Another  point  that  need  to  be  addressed  is  the  infinite  limits  of  integration  that  appear  in  Equation  (6.69). 
They  arise  from  Equation  (2.42)  where  the  infinite  limits  of  integration  are  permitted  because  of  the  actual  finite 
time  extent  of  the  integrand.  (Recall  that  a  radar  will  employ  time  gates  to  extract  time  segments  of  the  received 
signal  equal  in  duration  to  the  radar  pulse  width.)  In  the  absence  of  estimation  errors,  Equation  (2.42)  would  be 
expressed  as 

i  -  ?  [sDP  (t--rB)]*  STP(t)  e~j2*fB(t~TDP“rB)  dt  .  (6.71) 


Under  the  conditions  of  Equations  (5.108)  and  (5.109)  both  the  main  radar  and  interferometer  signals  contained 
within  sDP(t-rB)  and  sTP(t)  are  zero  outside  the  time  interval 

Ttp  <  t  <  rTP  +  NT  .  (6.72) 

Consequently,  Equation  (6.71)  can  be  written  as 


.TP  +  NT 


T 

l 

rTP 

[sDP(t-rB)]*  sTP(t)  e-j2*fB("rTP)dt 

(6.73) 

where  use  was  made  of  Equation  (2.3).  Introducing  the  change  of  variable 

ii 

1 

H 

(6.74) 

Equation  (6.73)  becomes 

NT 

[sdp(t  +  ttp  -tJ  sjp(t  +  ttp)  e-j2*fBp  d t 

0 

NT 

(6.75) 

=  1 

[sDP(t  +  rDP)r  sTP(t  +  TTP)  e-j2"fBt  dt 

where  use  was,  again,  made  of  Equation  (2.3)  and  the  variable,  t,  was  substituted  for  the  variable,  x.  Focusing  on  the 
noise  components  of  Equations  (2.28)  and  (2.35),  the  component  of  £  corresponding  to  I7  is  given  by 

NT 

I  =  J  [nDP  (t+rDP)]*  nTP(t+rTP)  e-j2,rfBt  dt  .  (6.76) 

0 
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Equation  (6.76)  is  now  used  to  determine  an  expression  for  e[|  I7  |  2  ] . 

Utilizing  Equation  (6.76),  it  follows  that 

£[|/,|2]=y  £{k,«+T0,)r  nDr(u+zBP)} 

E{nTr(t+ tTP)  [nTP(u+fTr)] }  dtdu  (6.77) 

NT 

=  JJ  RDP(u-t)  Rw(t-u)  gV'2*/s(“-o  dtdu  . 

0  "  " 

Observe  that  Equation  (6.77)  is  identical  to  Equation  (6.69)  except  that  the  limits  of  integration  are  now  finite.  Note 
from  Equation  (6.70)  that  the  noise  correlation  function  is  an  even  function  of  t.  Consequently, 

R[  TP  (t-u)  =  Rn7T  (u  -t)  .  (6.78) 

Employing  Equations  (6.70)  and  (6.78),  it  follows  that 

e[||7|2]=  (2N„W)!  fl  dt(Ju  _  (6.79) 

U  71  0  o  l  2ffW(u-t) 

Attention  is  now  devoted  to  evaluation  of  the  double  integral  in  Equation  (6.79).  Introduce  the  change  of 

variable 

r  =  U  - 1  .  (6.80) 

We  wish  to  replace  the  double  integral  involving  u  and  t  by  a  double  integral  involving  x  and  t.  The  regions  of 
integration  in  die  t  vs.  u  and  t  vs.  x  planes  are  shown  in  Figure  6.1.  Observe  that  the  boundaries  for  the  region  of 
integration  in  the  t  vs.  x  plane  are  given  by 


I: 

ii 

2 

H 

X  -  u  - 

II: 

t  =  o  , 

X  =  u 

III: 

u  =  NT  , 

x  =  NT 

IV: 

u  =  0  , 

x  =  -t 

t 

NT 

i 

0  jl  NT  u 


t  i 

i 

NT 

1  rv  n. 

1 _ \ 

in 

\  , 

-NT 

0 

n 

NT 

Figure  6.1 :  E  [  1 1 7  |  2  ]  is  evaluated  by  integrating  over  the  regions  indicated. 


Carrying  out  the  integration  in  the  t  vs.  t  plane, 


ivi 

E  /7 12  =  (2 N0W)2  1 

*  i\ 


NT 

l 

0 


sin2;r  W z 
2kWZ 


f NT-t  I 

eJlKf»x\  f  dt\  dz 


-NT  L 

(2N0  Wf 

NT 

\ 

0 

sin 2k  Wz 

2k  Wz 

(2  N„Wf 

0 

1 

sin  2k  W  z 

-NT 

2  kWz 

(2  N0Wf 

NT 

\ 

-NT 

sin 2k  W z 

2  kWz 

sin  2k  Wz 

2kWz 


-|2 


eJ2nf°T\ T  dt\  dt 


-|2 


~\2 


ejl7Cf*T  (NT -z)  dr 
ej2nflT  (NT +z)  dz 
(7Vr-|r|)  ejlnf°x  dz 


(6.82) 


where  the  strips  of  integration  with  respect  to  t  are  shown  in  Figure  6. 1 .  Unfortunately,  even  though  the  expression 
for  E[|  I7  | 2  ]  has  been  reduced  to  a  single  integration,  a  closed  form  result  does  not  appear  possible. 


To  obtain  an  estimate  of  E[|  I7 1  2  J ,  a  rectangular  approximation  is  made  for  sin  2  n  W  t  /  2  n  W  t.  The 

height  of  the  rectangle  is  chosen  to  be  unity,  which  is  the  maximum  value  of  sin  2  n  W  x  /  2  n  W  x,  and  the  width  of 
the  rectangle  is  chosen  to  be  1/W,  which  is  the  width  of  the  main  lobe  of  sin  2  it  W  t  12  7t  W  x.  The  rectangular 
approximation  is  centered  at  the  origin  because  the  main  lobe  of  sin  2  n  W  x  12  7t  W  x  is,  also,  centered  at  the  origin. 
Hence,  we  introduce  the  approximation 

1  „  ^  1 

1, - <Z< - 

2W  2W  (6-83) 

0,  otherwise 

Consequently, 


sin  2k  Wz 
2kWz 
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-1  2  W 

£[|/7| j  «  {2N0W)2  J  (NT -\r\)  ejlMfst  dr 


=  (2 N0W)2  1  (NT +t)  eilKhz  dr 


z.rr 

+  (2NqW)2  \  ( NT-r )  eJ2rtfsT  dr 


Utilizing  the  integral  relation  that 


J  r  ear  dr  =  ^-(ar  -1) 


and  combining  results,  it  follows  that 


;[|I7|2U  (2N0W)2  - - l—r  f"l  -  cos-^-^-1  +  ( NT - — 

L|  7|  V  0  ’  2 (zr f B ) 2  L  W  J  ^  2^ 


.  rcfu 


2W  j  nU 


From  Equation  (4.65)  and  (5.26)  the  baseband  bandwidth  of  both  the  Costas  and  OFDM  signals  is  approximately 
given  by 

N 

BW  =  —  (6.87) 


where  the  conditions  of  Equations  (5.108)  and  (5.109)  are  assumed  to  apply.  Assuming  the  receiver  bandwidth  of 
the  non-cooperative  bistatic  radar  receiver  is  twice  the  baseband  bandwidth  of  the  signals. 


w-H. 

T 


As  a  result,  Equation  (6.86)  can  be  expressed  as 


:[|i7|2]  - 


oxr  NY  t  f*fBT 

2Nn—  1  -  cos  — — 


2n  fn 


f  T  .  [  ^fBT  ) 

+  NT - sin  - —  >  (6.89) 

2N  N  Jf 


where  fB  is  the  bistatic  Doppler  shift  given  by  Equation  (2.38). 


Finally,  following  the  same  reasoning  that  led  to  the  results  of  Equations  (6.65)  and  (6.66),  it  is  concluded 
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E[|l«|!]=|k^f  N0 


E[|l,|2]  =  |kS,r  N0 


With  reference  to  Equations  (6.47),  (6.48),  (6.50),  (6.51),  and  (6.53),  it  is  noted  that  all  terms  are  zero  in  the 
double  summation  of  Equation  (6.35)  for  which  i  *  j.  Hence, 


e[|<,|!]=  x  e[|i,|2] 


Substitution  of  Equations  (2.45),  (6.44),  (6.45),  (6.60),  (6.65),  (6.66),  (6.89),  (6.90),  and  (6.91)  into  Equation  (6.92) 
yields 


»  I2  _ 

■i  |  ~ 


\kDP  I2  \kTP  I2 

\"~RAD  y^IFM 


,DP  2  iTP  2 

kifm  krad 


DP  I2  I;  TP  I2 


\klFM\  \UA 


+  (N-iy  +  kj^  n0 


+  Km\2  N0  +  I 2  N0  +  k°pM  2  N0 


1-cos 


*f,T 


2nfB 


( XTrr  T\  .  (nfBT 

+  NT - sin  JB 

l  2  N  l  N 


On  the  other  hand,  making  use  of  Equations  (2.45)  and  (6.32), 


eM!=N2  =|k™f 


The  variance  of  th  as  given  by  Equation  (6.33),  then  becomes 
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Var «,}-  \k 


1  -  cos 


One  measure  of  performance  for  a  radar  is  signal-to-noise  ratio.  For  a  monostatic  radar  operating  in  an 
additive  thermal  noise  background  with  variance,  N0,  the  signal-to-noise  ratio  is  given  by 


(SNR)„. 


_  i4<,r 

nos.a.U  Var{Noise}  Nq 


(6.96) 


For  the  non-cooperative  bistatic  radar,  the  signal-to-noise  ratio  is 


,  ,  |  E[es]2 

[SNR  )Non~cooperative  =r  - 

'Bistatic  Var  {l  J  } 


(6.97) 


where  Var{£[}is  given  by  Equation  (6.95).  The  loss  factor  is  defined  to  be 


noise  ratio  of  the  monostatic  radar  is  expected  to  be  larger  than  that  of  the  non-cooperative  bistatic  radar.  The  more 
effective  is  the  interferometer  in  masking  the  transmitted  radar  waveform,  the  larger  will  be  the  value  for  LF. 


6.3  Non-Cooperative  Bistatic  Radar  Receiver  Performance 


6.3.1  False  Alarm  Probability 


Because  Z(0,0),  T(0,0),  and  fiFM(t)  consist  of  a  large  number  of  random  components  due  to  the  randomness 
of  the  data  symbols,  the  central  limit  theorem  can  be  invoked  to  justify  modeling  the  detection  statistic,  t,  as  a 
complex  Gaussian  random  variable.  By  definition,  the  probability  of  false  alarm  is  the  conditional  probability  of 
detecting  a  target  given  that  a  target  is  not  present.  With  reference  to  Equations  (2.35),  (2.36),  and  (2.37),  the 
absence  of  a  target  is  equivalent  to 

=  k-RAD  =  ^ifm  =  0  •  (6.99) 

Then 


sTP(t)  =  nTP(t)  .  (6.100) 

The  mean  value  of  the  signal  component  of  the  detection  statistic  is  given  by  Equation  (6.28).  In  view  of 
Equations  (2.45)  and  (6.99),  it  follows  that 

E[/.]  =  ki  =  0  (6.101) 

in  the  absence  of  a  target.  In  addition,  the  mean  value  of  the  interference  component  of  the  detection  statistic  is 
given  by  Equation  (6.32).  Hence,  in  the  absence  of  a  target,  it  is  concluded  from  Equations  (2.45)  and  (6.99)  that 
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(6.102) 


E[/,]=k4  =  0 

Thus,  given  that  a  target  is  not  present, 

e[/]  =  e[/J+ e|a]  =  o  .  (6.io3) 


Let  <jj0  denote  the  variance  of  the  detection  statistic  in  the  absence  of  a  target.  Because  is  is  a  constant, 


<  =  Var{e,}=  |C„f  K  +  |Crf 


(6.104) 


( 2N0N/T )2 

*fB 


1  -  cos 


N 


2  wfB 


( 

NT  -  — 
\  2N 


dn*f*T 


N 


where  use  was  made  of  Equations  (6.99)  and  (6.95).  The  absence  of  a  target  is  conventionally  referred  to  as 
hypothesis,  H0.  Therefore,  under  H0,  the  detection  statistic,  £,  can  be  modeled  as  a  zero-mean  complex  Gaussian 
random  variable  with  its  variance  given  by  Equation  (6. 104). 


and 


Denote  the  real  and  imaginary  parts  of  the  detection  statistic  by 

ix  =Re{f} 

=  Im{/} 


In  the  absence  of  a  target,  the  probability  density  function  of  £  is 

Pt\H,  (L=Lr+  J  L,\Hq)  =  P(  >tilHo  (Lr,  L,  |tf0) 

_  wl 

1  la1,.  1  2  a. 


(4  +  4) 


2no: 


2no, 


where  (7 >a  is  given  by  Equation  (6. 104). 


(6.105) 

(6.106) 


(6.107) 


Let  Hj  refer  to  the  hypothesis  that  a  target  is  present.  The  statistical  test  employed  for  detection  is  then 
described  by 


|z|  >  r  (6.108) 


where  y  is  a  constant  that  denotes  the  threshold. 


A  statistically  equivalent  test,  which  is  easier  to  analyze,  is 
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L\2  =  L\  +  L)  >  f 


(6.109) 


As  a  result,  the  probability  of  false  alarm  is  given  by 

pF  =  Pr  {el  +e]>  r2|tf„}=  Pr{er,e,  e  r  |h0} 

where  R  is  the  shaded  region  sketched  in  Figure  6.2.  Consequently, 

PF  =  II  Ptr,t,  |h0(4  >A  dLi 


(4  +  4) 


Jf  2^(7, 


e  dLdL , 


(6.110) 


(6.111) 


Figure  6.2:  A  target  is  declared  present  when  |L|  falls  into  the  shaded  region. 


To  evaluate  the  double  integration  in  Equation  (6.1 10),  it  is  convenient  to  introduce  polar  coordinates.  Let 


Lr  -  Z  COS  p 


(6.111) 


L;  =  Z  sin 


(6.112) 


It  follows  that 


z  =  vi/r  +  M 


(6.113) 


P  =  tan  1  — 
Lr 


(6.114) 


Equation  (6.1 10)  can  then  be  expressed  as 
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271(7]  -* 


I  d/?  J  Z  e  2<  dZ 


where  use  was  made  of  the  differential  relationship 

dLr  dLj  =  ZdZ  dyS 

Evaluation  of  the  integrals  results  in 


1  f  j  o  2  2  a? 

PF  = - r  J  dP  e  ' 

F  2;r<r2  -* 


(6.115) 


(6.116) 


(6.117) 


=  JL  e  2<  f  dp  =  e  la'° 

In 


From  Equation  (6. 1 1 7)  it  follows  that 


7  =  (7 1  J 2  In  (1/PF) 


(6.118) 


Consequently,  the  threshold  is  readily  determined  from  a  knowledge  of  <j\0  and  specification  of  the  false  alarm 
probability.  Note  that  the  smaller  is  PF  and  the  larger  is  <7(0  the  larger  is  the  required  value  of  the  threshold. 

6.3.2  Detection  Probability 

By  definition,  the  probability  of  detection  is  the  conditional  probability  of  detecting  a  target  given  that  a 
target  is  present.  The  presence  of  a  target  is  typically  designated  as  hypothesis,  H] .  Under  H]  the  detection  statistic, 
t,  is  modeled  as  a  complex  Gaussian  random  variable  with  mean 


mti  =  +  j  mti 

and  variance  c]l  .  From  Equations  (6.28),  (6.32),  and  (2.45)  it  follows  that 

mti  =  Eb.frhBk'fri-t&y  kfLAD  +  te)‘  K 

<J2n  is  given  by  Equation  (6.95). 


(6.119) 


(6.120) 


In  the  presence  of  a  target,  the  probability  density  function  of  l  is 


2  no 


(6.121) 


However, 
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(6.122) 


L~mt!\  =  [(4  ~me)  +  j  (4  ~mi, )]  [(4  ~mir  )  -7(4  -«/,)] 


=  K  +  4r  +  4,  -20»*,  4  +  4) 


As  a  result,  Equation  (6.121)  becomes 


P,|„,(L|Hi) 


2kg) 


m<r  +  L2  +  L2  -  2(m<f  L,  +  m,.  L,) 

e"  e"  2£7'2' 


(6.123) 


Employing  the  test  specified  in  Equation  (6.109)  the  probability  of  detection  is  given  by 

4>  =  Pr{4+4  >r2  Pr{4, 4.  £  /?!//,} 


(6.124) 


where  R,  once  again,  is  the  shaded  region  sketched  in  Figure  6.2.  Consequently, 

4  =  y  j#, (4  >4|#i)  <*4  ^4 


4+ I2  -2 (m/fZ.r  +  m,,4) 


(6.125) 


2;r<72 


<7  4  <7  L, 


Evaluation  of  the  double  integral  in  Equation  (6.125)  is,  as  before,  facilitated  through  introduction  of  the 
polar  coordinates  defined  by  Equations  (6. 1 1 1)  and  (6. 1 12).  It  follows  that 


m2  rr,2 

m/r  + 


po  =  — ^  e 

2k 


/  Z  e  "1  J 


Z(m<r  cos/?  +  m,.  sin /?) 


d/3  dZ  .  (6.126) 


From  a  table  of  integrals  it  is  known  that 

1  71 

1  f  _  r  co 


—  I  eims/,  +  ssinfi  d/3  =  10  [V r2  +  s2 

2k  -*■ 


(6.127) 


where  I0(-)  is  a  zeroth-order  modified  Bessel  function  of  the  first  kind.  Therefore,  the  detection  probability  becomes 


2  ,  2 

mtr  +  mtj  _  *  r 

1  „  2£7'i  7  ^  „  2<7'i  r  ^  F3  nr 


4,  =  -V  *  2"'  J  Z  e  I0 


.2  \’"tr  '  ”‘t, 


mt  +  m(  dZ  .  (6.128) 


With  reference  to  Equation  (6. 119),  observe  that 


=  m2  +  m2 


(6.129) 


Thus,  the  expression  for  PD  can  be  simplified  to 
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PD  =-L  e"2"1'  J  Ze"“>  I0  -^r  Z  dZ  . 


(6.130) 


Having  integrated  out  p,  it  has  been  determined  that  the  probability  density  function  on  z  is  given  by 


(6.131) 


/>.(*)  H/r* 


—  e  2<T'-  Z  e  2<  /0  Z  ,  Z  >  0 

O’,2  01, 


,  Z  <  0 


Consequently, 


1  Pz(z)  dz  =  1 


(6.132) 


and  the  probability  of  detection  can  be  expressed  as 


Pd='  ~  -r  e 

CF  p 


k  z  r 

72  v  -rrr  m,  . 

2ff'‘  \  Ze  2°1'  10  — y  Z  dZ 

0  O, 


(6.133) 


The  above  result  can  be  further  simplified  by  introducing  the  incomplete  Toronto  function  [13,14]  which  is 
defined  as 


TB(m,n,r)  =  2r“'  e"'‘  f  t”  e'z’  I„[2rt]  dt  . 


(6.134) 


V2  <7 1 


(6.135) 


V2  <7,  * 

Substitution  of  Equations  (6.135)  and  (6.  136)  into  Equation  (6.133)  yields 


(6.136) 


PD=  \-2e~r2  !  t  e~‘2  I0  [lrt]dt 


a 


1  T  j  o 

G ,  r\  ’  '42  G„  • 


(6.137) 
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Plots  of  (1,0,  ->Jq  )  are  presented  in  Figures  13  and  14  of  the  Marcum  paper  [13]  for  various  choices  of  u 
and  q.  Although  the  Toronto  function  needs  to  be  evaluated  numerically,  Equations  (6.130)  and  (6.133)  clearly 
reveal  that  the  detection  probability  is  decreased  by  increasing  y  and  <7^  ,  This  observation  is  useful  in  selecting 
system  parameters  to  effectively  degrade  the  operation  of  the  non-cooperative  bistatic  radar. 
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CHAPTER  7 


CONCLUSIONS 


7.1  Summary 

Radar  receivers  typically  utilize  a  replica  of  the  transmitted  waveform  for  processing  the  radar  return.  In  an 
attempt  to  estimate  this  transmitted  waveform,  non-cooperative  bistatic  radars  normally  intercept  a  direct  path  signal 
emitted  through  the  side  lobes  of  the  host  radar. 

This  report  investigates  the  use  of  an  interferometer  along  with  the  host  radar  to  deny  coherent  reference 
signals  to  unwanted  observers.  The  interferometer  signal  is  intended  to  be 

1)  effective  in  masking  the  host  radar  signal  emitted  through  the  side  lobes, 

2)  orthogonal  to  the  signal  transmitted  through  the  main  lobe  of  the  host  radar  so  as  not  to  degrade  radar 

performance, 

3)  useful  for  data/voice  communications, 

4)  controllable  such  that  separate  communication  signals  can  be  directed  to  different  receivers  located  in 

different  directions. 

In  addition,  the  host  radar  signal  is  intended  to  be  flexible  to  allow  for  multi-mission  operations. 

The  time  domain  cannot  be  used  to  achieve  orthogonality  between  the  interferometer  and  radar  waveforms 
because,  assuming  perfect  synchronization  and  temporal  orthogonality,  the  non-cooperative  radar  could  correlate  the 
direct  path  signal  with  its  radar  return  without  experiencing  any  loss  due  to  presence  of  the  interferometer  signal.  In 
fact,  for  effective  masking,  the  radar  and  interferometer  signals  should  be  highly  correlated  in  the  time  domain. 

Consequently,  the  desired  orthogonality  must  be  achieved  in  the  spatial  domain.  Analysis  of  the  two- 
element  interferometer  array  factor  reveals  that  a  relatively  broad  spatial  null  exists  in  its  end  fire  position  when  the 
element  spacing  is  an  odd  integer  number  of  half  wavelengths  of  the  transmitted  interferometer  signal.  This  null 
cannot  be  steered  electronically.  Therefore,  when  the  host  radar  main  lobe  points  broadside  to  the  x-axis,  the 
interferometer  elements  must  be  placed  along  the  y-axis  in  order  to  take  advantage  of  the  broad  null.  Should  the 
main  lobe  be  steered  to  different  directions,  it  would  be  necessary  to  mechanically  rotate  the  interferometer  elements 
in  order  to  maintain  the  spatial  orthogonality  between  the  interferometer  and  radar  waveforms. 

In  order  to  communicate  separate  signals  in  different  directions,  it  will  be  necessary  to  steer  the 
interferometer  without  destroying  the  broad  null  in  the  end  fire  position.  This  can  be  accomplished  by  frequency 
hopping  the  transmitted  carrier  such  that  the  odd  integer  relationship  between  the  element  spacing  and  wavelength  is 
maintained  while  the  carrier  changes  from  one  frequency  to  another. 

The  flexibility  of  frequency  hopping  makes  it  attractive  for  design  of  the  host  radar  waveform.  Depending 
upon  the  particular  mission,  the  frequency  hops  can  be  arranged  to  synthesize  a  discrete  frequency  quantized 
approximation  to  the  frequency  modulated  signal  that  produces  a  desired  delay-Doppler  ambiguity  function. 

Because  the  ideal  shape  of  the  ambiguity  function  is  an  impulse  function,  the  Costas  signal  is  proposed  for 
the  host  radar  waveform.  It  yields  a  thumb  tack  shaped  ambiguity  function  with  a  relatively  low  pedestal.  For  a 
fixed  number  of  frequency  hops  within  a  radar  pulse,  there  are  many  different  hopping  patterns  that  result  in 
essentially  the  same  thumb  tack  shaped  ambiguity  function.  Hence,  different  frequency  hopping  patterns  can  be 
utilized  to  further  complicate  the  coherent  reference  estimation  task  of  the  non-cooperative  radar. 

Having  selected  the  Costas  signal  for  the  transmitted  radar  waveform,  orthogonal  frequency  division 
multiplexing  (OFDM)  becomes  an  attractive  communication  technique  for  the  interferometer.  The  hopping 
frequencies  of  the  Costas  signal  are  required  to  be  orthogonal  over  a  hopping  interval  of  duration,  T.  The  subcarrier 
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frequencies  of  an  OFDM  signal  are  required  to  be  orthogonal  over  a  symbol  duration,  Ts.  By  selecting  T  =  Ts,  the 
two  sets  of  frequencies  become  identical.  In  addition,  the  use  of  M-ary  amplitude  and  phase  shift  keying  to 
modulate  the  subcarriers  makes  it  very  difficult  to  separate  the  Costas  and  OFDM  signals  when  their  amplitude 
transitions  are  synchronized  in  time. 

Utilization  of  OFDM  for  communications  has  many  important  advantages.  Both  modulation  and 
demodulation  can  be  implemented  using  fast  Fourier  transforms.  The  difficult  tasks  of  channel  estimation  and 
equalization  are  simplified  significantly  when  communicating  over  fading  and  time-variant  broadband  channels. 
Also,  the  parallel  nature  of  OFDM  makes  it  an  efficient  scheme  for  the  transmission  of  data. 

Performance  of  the  non-cooperative  bistatic  radar  receiver  was  analyzed.  In  addition  to  noise  terms,  the 
detection  statistic  was  shown  to  involve  the  delay-Doppler  ambiguity  functions  of  the  host  radar  and  interferometer 
waveforms  as  well  as  their  cross  ambiguity  function.  Analytical  expressions  for  these  functions  were  obtained  under 
the  assumption  that  the  host  radar  and  interferometer  employ  Costas  and  OFDM  signals,  respectively.  Taking 
advantage  of  the  result  that  the  detection  statistic  consists  of  a  random  sum  of  many  terms,  none  of  which  are 
dominant,  the  detection  statistic  was  approximated  as  a  complex  Gaussian  random  variable.  Analytical  expressions 
were  then  derived  for  the  false  alarm  and  detection  probabilities  of  the  non-cooperative  bistatic  radar. 


7.2  Suggestions  for  Future  Work 

7.2.1  Performance  Evaluation  of  the  Complete  System  Using  Realistic  Scenarios 

Although  analytical  expressions  needed  for  performance  of  the  non-cooperative  bistatic  radar  were 
obtained,  time  did  not  permit  for  assessment  of  the  complete  system.  It  is  desirable  to 

1)  Construct  plots  of  the  OFDM  ambiguity  function  and  the  cross  ambiguity  function  between  OFDM 
and  Costas  signals  for  both  individual  sample  functions  and  the  entire  ensemble  of  the  OFDM  signal, 

2)  Obtain  analytical  expressions  for  the  false  alarm  and  detection  probabilities  of  the  host  radar, 

3)  Postulate  realistic  scenarios  along  with  the  corresponding  numerical  values  for  the  system  parameters, 

4)  Evaluate  complete  system  performance  for  the  realistic  scenarios  postulated  in  (3)  and  adjust  the 
system  parameters  where  necessary  to  improve  performance, 

5)  Perform  computer  simulations  to  verify  the  predicted  theoretical  results. 

7.2.2  Extension  of  the  Interferometer  to  More  than  Two  Elements 

The  two-element  interferometer,  under  a  specified  constraint,  produced  a  broad  null  in  the  end  fire  position 
which  was  useful  in  obtaining  spatial  orthogonality  between  the  host  radar  and  interferometer  waveforms.  However, 
this  null  could  not  be  steered  electronically.  It  is  conjectured  that  introduction  of  a  third  interferometer  element  will 
enable  steering  of  the  broad  null.  Should  this  prove  to  be  possible,  the  advantages  to  be  gained  by  adding  additional 
elements  would  be  investigated. 
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The  concept  of  the  tri-element  interferometer  is  illustrated  in  Figure  7.1.  Elements  1  and  2  can  be  replaced 
by  an  equivalent  element,  denoted  by  E.  Elements  E  and  3  then  comprise  an  equivalent  two-element  interferometer. 
By  adjusting  the  amplitudes  and  phases  of  the  sinusoidal  signals  applied  to  elements  1  and  2,  the  position  of  element 
E  can  be  moved  horizontally  along  the  line  connecting  elements  1  and  2.  Selection  of  the  frequency  of  excitation 
such  that  the  distance,  d,  is  an  odd  integer  of  half  wave  lengths  would  then  produce  a  broad  null  in  the  end  fire 
position  of  the  equivalent  two-element  interferometer.  In  this  manner  the  broad  null  could  be  electronically  steered 
by  appropriately  selecting  the  amplitudes  and  frequencies  of  the  sinusoidal  signals  applied  to  elements  1, 2,  and  3. 

1  E  2 


d 


3 

Figure  7.1:  The  tri-element  interferometer  can  be  viewed  as  an  equivalent  two-element  interferometer. 


7.2.3  Generalization  of  the  Ambiguity  Function  to  Include  Direction  of  Arrival 

Because  both  spatial  and  temporal  processing  is  performed  in  a  radar,  an  ambiguity  function  generalized  to 
include  the  direction  of  arrival  would  be  extremely  useful.  As  a  first  step,  it  is  proposed  to  carry  out  this 
generalization  by  incorporating  the  antenna  radiation  patterns  as  a  function  of  angle. 

More  generally,  the  transformation  of  operators  approach,  introduced  by  Cohen  [15],  may  be  useful  in 
developing  a  combined  delay-Doppler,  spatial  displacement-spatial  frequency  ambiguity  function.  This  would  be 
applicable  to  pulse  propagation  where  there  is  a  strong  correlation  between  temporal  and  spatial  frequencies. 
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